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TRANSLATORS’ PREFACE 


The exciton concept was introduced into solid-state physics by Frenkel 
in 1931, in connection with the transformation of excitation energy into 
heat in rare-gas solids. Not until Davydov’s now-famous 1948 paper 
on the splitting of electronic terms in the naphthalene crystal was the 
concept applied to geometry-determined molecular problems. Since 
then the treatment of an impressive variety of problems has been pub- 
lished by Davydov and his colleagues and by numerous other researchers. 
The exciton model has now been used in treating general formulations 
of energy transfer, in the crystal splitting of electronic states which are 
nondegenerate in the isolated molecule, in the splitting of infrared bands 
and Raman lines in molecular crystals corresponding to nondegenerate 
vibrational modes in the isolated molecule, in the calculation of excited 
states of molecules consisting of two or more isolated or not-strongly- 
conjugated chromophoric units, in optical rotatory power of helical 
polymers containing nonconjugated chromophoriec units, etc. Alto- 
gether, the range of problems covered suggests that the role of the exciton 
model may be as great in the fields of molecular physics and quantum 
chemistry as that of the molecular orbital model, the valence-bond model, 
or the charge-transfer model. Davydov’s book summarizes many of 
the problems, and its perceptive treatment suggests many starting points 
for further extensions. 

The translators are pleased to make this unique book available in 
English. In spite of widespread current interest in the topics covered, 
copies of the book in Russian are something of a rarity, and even its 
existence is unknown to many researchers. 

We wish to give a few notes on our translation. Throughout the book 
we have translated: 


“Combinational scattering spectra” as combinational scattering (Raman) 
spectra 

“Partial vibrations” and ‘normal vibrations” as normal vibrations 

“Long and intermediate axes” as long and intermediate axes 


In the first case, the Russian terminology was preserved, but the 
more familiar terminology was introduced. In the last case, the author’s 
Vv 
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usage was again preserved in preference to the more common “‘long and 
short axes.”’ 

We have not attempted to revise the form of the mathematical expres- 
sions except for the replacement of some half-dozen Cyrillic characters by 
Roman characters. A number of typographical errors among symbols 
in the original text have been corrected. A few fundamental corrections 
supplied by the author have been included. 

This book was published originally under the title “Theory of the 
Absorption of Light in Molecular Crystals” by A. 8. Davydov, Works of 
the Institute of Physics, Ukrainian Academy of Sciences, Kiev, Ukrainian 
§.8.R., 1951. We feel that the title we have adopted stresses the princi- 
pal theme of the book—the resonance interaction of excited states of 
weakly coupled molecular systems. 

To bring this edition more up to date, we have included two appendixes 
which cover principal literature on related topics to December, 1961. 

We acknowledge with thanks the careful reading of the final proof and 
suggestions by Dr. Robin M. Hochstrasser. 


Michael Kasha 
Maz Oppenheimer, Jr. 
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INTRODUCTION 


1. Molecular Crystals 


For the theoretical description of the physical properties of a solid 
body, it is convenient to classify crystalline substances according to the 
magnitude of the ratio of the forces of interaction of the particles (atoms 
or molecules) of which the crystal is formed to the forces of interaction 
between the composite parts (electrons and atoms, respective'y) of these 
particles. If the forces of interaction between the particles are signifi- 
cantly smaller than the forces of interaction between the composite parts 
of one particle, then when forming a solid body these particles preserve 
their individuality in the first approximation. We shall define such solid 
bodies as molecular crystals. Usually molecular crystals are composed 
of molecules between which van der Waals forces of attraction act within 
the crystal. If the energy of intermolecular interaction is close in order 
of magnitude to chemical energies (the bond energies of atoms in a mole- 
cule or of electrons in an atom), then the molecules lose their individu- 
ality, and the crystals can now be identified as of the valence, metallic, 
or ionic type. 

Molecular crystals are formed from atoms of inert gases (Ne, A, Kr, 
and others) and from molecules with saturated bonds (H2, CHy, and 
others).f Molecular crystals have an especially wide occurrence among 
organic substances. 

Molecular crystals are characterized by low melting and boiling points. 
Evaporation takes place by whole molecules. The heat of sublimation 
of typical molecular crystals amounts to several kcal/mole, i.e., is 10 to 
20 times smaller than the dissociation energy of isolated molecules or 
than the dissociation energy of an ionic crystal lattice (Born, Ref. 1, 
p. 309). 

Molecular crystals are electrical nonconductors. Almost all molecular 
crystals are diamagnetic (except, for example, crystals of O2 and NO), 
inasmuch as they consist of closed-shell molecules having no resultant 
spin. 

In this work we shall investigate the optical behavior of molecular 

{ Translators’ note: Crystals formed of molecules containing unsaturated bonds 


(Ox, benzene, naphthalene) are also included in this category. 
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crystals and, more specifically, the problems of light absorption in molecu- 
lar crystals. 


2. Method and Aims of the Investigation 


The theoretical examination of the question of light absorption by a 
solid body (as well as by other aggregates) leads to the solution of two 
problems: (1) the determination of the stationary states of the crystal and 
(2) the calculation of the transition probability, under the action of a light 
wave, from one state to another, or (a less complete but simpler problem) 
the determination of the selection rules for the optical transitions. 

In general, the solution of these problems for any solid body cannot be 
completed in practice, inasmuch as the problem narrows down to solving 
the equations of the quantum mechanics of stationary states for a system 
with an enormous number of degrees of freedom. The solution of this 
type of problem can be achieved only by means of approximation meth- 
ods. In this connection, the approximation method utilized depends 
upon (1) the region of wavelengths in which it is necessary to investigate 
the absorption and (2) the character of the solid body (metals or molecu- 
lar, valence, or ionic crystals). 

When investigating absorption in the infrared region of the spectrum, 
one can, with sufficiently good approximation, examine only the move- 
ment of the atomic nuclei in the lattice, utilizing classical mechanics. 
When investigating the absorption spectra in the visible and ultraviolet 
regions, coupled with vibrations of the nuclei, one must investigate the 
change in electronic states of the crystal, but in this case, it is necessary 
to utilize quantum mechanics. The methods used for the calculation of 
molecular crystals will be essentially different from the methods for solv- 
ng the same problems for other types of crystals. 

When calculating the optical properties of molecular crystals, it is 
natural to start with the premise that these properties are determined 
basically by the optical behavior of isolated molecules and to consider 
the influence of their interaction as a perturbation. Such an assumption 
is useful from two points of view. First, it points out the method of 
solving the problem, namely, by utilizing quantum-mechanical pertur- 
bation theory. Second, it has great practical significance. If it is possi- 
ble to determine how the optical properties of molecules change under 
the influence of intermolecular forces when molecules are united into 
crystals, then the investigation of the absorption spectra of crystals 
enables one to decipher the absorption spectrum of the free molecule, the 
structure of the crystal itself, and the nature of the molecular interaction. 

It is known that the study of the character of the interaction of matter 
with radiation is one of the methods of investigating the structure of 
matter. It is sufficient merely to remember that the structure of atoms, 


INTRODUCTION 3 


and also, to a certain degree, of molecules, became known by virtue of 
the development of methods of spectroscopy for the gaseous state. The 
large amount of experimental information assembled served as a basis 
for creating quantum concepts and, later, quantum mechanics. In the 
spring of 1928, L. I. Mandelshtam and G. S. Landsberg? (for crystalline 
quartz) and C. V. Raman (for liquid benzene) discovered the comb:- 
national scatteringt of light. This discovery soon took a prominent place 
as a& means of investigating molecular structure. The combinational 
scattering (Raman scattering) enables one to obtain the vibrational fre- 
quencies of atoms in molecules and the frequencies of intermolecular 
vibrations in solid and liquid bodies. Vibrations of the latter type were 
first discovered by the Soviet physicists E. F. Gross and M. F. Vuks and 
named by them low-frequency spectra.* 

One can expect that the study of absorption spectra and of light scatter- 
ing of crystals also will be one of the methods for studying their structure 
and the forces of intermolecular interaction as well. However, in the 
visible and ultraviolet regions of the absorption spectrum, the majority 
of solid bodies at room temperature exhibit wide bands without evident 
structure. Salts of the rare-earth elements, salts of chromium and man- 
ganese, aS well as a few others, form an exception: the absorption spectra 
of these substances show a large number of narrow bands. The anoma- 
lous behavior of the salts of rare-earth elements is conditioned by the 
following: When absorbing light, there is a change in the energy states 
of the 4f inner-shell electrons, which do not participate in valence bonds 
and are therefore protected by valence electrons from the effect of other 
atoms. The large number of absorption lines is explained by the split- 
ting of the degenerate energy levels by the electric field (created by 
neighboring ions); the symmetry of this field is determined by the sym- 
metry of the crystal. Splittings of this type are similar to the Stark effect 
and were qualitatively investigated by H. Bethe. 

J. Becquerel® first studied the spectra of the salts of rare-earth ele- 
ments at the temperature of liquid air and showed that a significant 
narrowing of the absorption bands occurred when the temperature was 
lowered. At the same time, V. Kurbatov® showed that the color of a 
very large number of crystals changed upon immersion in liquid air. 
This demonstrates the change of the absorption spectra of crystals when 
the temperature is lowered. 

Very significant are J. V. Obreimov’s works,’—® in which it was shown 
that, when the temperature is lowered, the absorption bands of many 
molecular crystals become discrete. 

These works supplied a basis for the systematic investigation of absorp- 
tion and fluorescence spectra of molecular crystals. 


t Translators’ note: Commonly called Raman scattering. 
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The lowering of crystal temperature changes the absorption spectrum 
essentially. Eight wide absorption bands of benzene, superposed at room 
temperature, transform?*!! at 78°K, into a series of separate absorption 
maxima each of which exhibits a complicated structure at 20.4°K. 

As was shown by I. V. Obreimov and A. F. Prikhotko,’? 3 wide and 
diffuse absorption bands of naphthalene at room temperature are trans- 
formed into 40 absorption bands at a temperature of 78°K and into more 
than 300 absorption bands!3 at a temperature of 20.4°K. 

At low temperatures a discrete spectrum is observed even in those 
compounds which have a continuous spectrum in the gaseous state, mak- 
ing their analysis in this state impossible. As I. V. Obreimov showed, 
the diffuse absorption spectrum of the vapor of azobenzene is broken up 
into a multitude of very sharp absorption bands in the crystal at the 
temperature of liquid hydrogen and helium. The continuous absorption 
spectrum of the vapors of iodobenzene, according to the data of E. M. 
Bronstein, acquires a distinct structure in the crystal at 20.4°K. 

It is difficult to analyze the absorption spectra of complex molecules in 
the gaseous state because of the great vibrational and rotational struc- 
ture. When studying the absorption spectra of the same substances in a 
crystalline state at low temperatures, it is apparently possible to decipher 
the absorption spectrum of the molecules. Owing to the small inter- 
action forces between the molecules, one could expect only insignificant 
spectral changes in molecular crystals upon variation of the aggregate 
state. These changes can be limited to the displacement of the entire 
spectrum, to the disappearance of the rotational band structure, to the 
appearance of small lattice vibrations (Gross frequencies), and to the 
splitting of the absorption bands in the crystal. The last phenomenon 
has always been linked with the removal of degeneracy from degenerate 
states of free molecules. However, if the molecule did not have degener- 
ate states, then presumably the splitting of terms could not take place. 

It is roughly these considerations which furnished I. V. Obreimov and 
A. F. Prikhotko the basis for their numerous investigations of absorp- 
tion and fluorescence spectra of crystals.12~1® However, these assump- 
tions, which at first sight seemed obvious, required an experimental 
and theoretical basis. The author expresses his sincere gratitude to 
A. F. Prikhotko and I. V. Obreimov, who called his attention to the 
indispensability of analyzing this problem from the theoretical point of 
view. 

Experience has fundamentally confirmed the correspondence between 
the spectra of the gas and crystal. However, until recently, a detailed 
verification could not be made, since there were insufficient experimental 
data available (especially concerning vapor spectra). It is true that the 
great pleochroism of molecular crystals has been well known for a 
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comparatively long time. Thus, I. V. Obreimov’® showed, taking azo- 
benzene as an example, that the absorption spectrum for light in this 
crystal is markedly different for two mutually perpendicular directions 
of the electric vector of the impinging light wave. Furthermore, the 
pleochroism in the absorption spectra of crystals of anthracene,!” phe- 
nanthrene,!® and a series of others was observed. If we consider that 
almost all molecular crystals contain somewhat differently oriented mole- 
cules in one unit cell (most frequently, two or four molecules per cell), 
then it becomes clear that one cannot explain the pleochroism of the 
crystal without the assumption of changes in the absorption spectrum 
of the molecule. 

There are very few theoretical works devoted to the investigation of 
absorption spectra of light in molecular crystals. The first investigations 
along this line were completed by J. I. Frenkel!*-2! and somewhat later 
by Peierls.2? However, all these investigators examined only the cubic 
crystals of rare monatomic gases (neon, krypton, and others); therefore, 
a series of interesting physical results was lost in their work. In spite of 
this, the works of J. I. Frenkel are of great theoretical interest. Firstly, 
these works were the first to introduce the concept of the exciton, whose 
movement is equivalent to the migration of the excitation through the 
crystal from one atom to another. This concept turned out to be very 
fruitful and was utilized subsequently in a series of other works. Sec- 
ondly, it was shown in these works that, in the ideal crystal, light can 
excite only the excitons with wave numbers equal to the wave number 
of the light wave. The violation of this condition is linked with the 
lattice vibrations. Thirdly, the works of J. I. Frenkel have posed and 
partially solved the problem of the transformation of absorbed light 
energy into heat. 

The method of obtaining the energy states of a crystal from the states 
of isolated atoms was also utilized by D. I. Blokhintsev**-?* in some of his 
works devoted to establishing a theory of light absorption in heteropolar 
crystals. 

8. I. Vavilov?-2? very successfully enlisted the phenomenon of induc- 
tive resonance, which leads to the migration of excitation energy from one 
molecule to another, to explain the changes in polarization and duration 
of the excited state of molecules in solutions. It was possible to prove 
experimentally the presence of energy migration in luminescent sub- 
stances and in crystals which contain foreign atoms. In this connection, 
it was shown that the migration can take place over very significant dis- 
tances, exceeding by a million times the distances between two neighbor- 
ing molecules or atoms.#*31 

In the present work, we shall show that the resonance interaction of 
molecules, which leads to the migration of energy in the crystal, plays an 
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extremely essential role in explaining the peculiarities of absorption and 
luminescence spectra in molecular crystals. 

We shall apply the exciton concepts to molecular crystals of any 
desired symmetry, which are composed of complicated polyatomic mole- 
cules. In this connection, we shall show that it is thereby possible to 
explain a series of absorption and fluorescence spectral peculiarities in 
molecular crystals. We shall see that taking into account the resonance 
interaction of molecules in the crystal leads to often basic, albeit small, 
changes in the spectrum of the isolated molecule. In particular, with 
the help of these concepts one can successfully explain the large pleo- 
chroism of molecular crystals. The resonance interaction between 
molecules sometimes leads to the appearance in the crystal of crystal 
excited states, which are distinct from the corresponding states in the 
gaseous phase. The appearance of special crystal excited states, which 
as a rule are strongly polarized, is experimentally confirmed by work per- 
formed in the laboratory under the direction of A. F. Prikhotko, Corre- 
sponding Member of the Academy of Sciences of the Ukrainian 8.8.R. 

We shall see that the intramolecular vibrations also vary under the 
action of intermolecular resonance forces. Therefore, the infrared spec- 
tra and the combinational scattering (Raman) spectra must also reflect 
the crystalline structure of a solid body. In the case of combinational 
scattering (Raman scattering), noticeable changes need be expected 
only for the vibrations which are simultaneously active in the infrared 
spectrum. 


3. Scope of This Work 


The following summary of the contents of the various chapters will 
facilitate reading this book. 

Chapter I is of an auxiliary nature. Therein are stated the problems 
linked with the vibrations of atoms and molecules in a molecular lattice. 

Chapters II and III contain the classical and quantal theory of absorp- 
tion spectra in the infrared, visible, and ultraviolet regions, as well as 
several consequences related to the combinational scattering (Raman) 
spectra. We shall show that in crystals there occurs a splitting of energy 
terms into a number of components equal to the number of molecules in 
the unit cell. With the aid of group theory, the selection and polari- 
zation rules are established for certain crystals. 

Starting with Chapter IV, we shall give the absorption theory, which 
considers the possibility of the movement of molecules in the lattice. 
Chapter IV formulates the problem. It states the circumstances under 
which the adiabatic method is applicable to systems with many degrees 
of freedom. Therein are shown the conditions under which either free or 
“localized” excitons appear in crystals at the instant light is absorbed. 


INTRODUCTION 7 


Equations are derived which determine the movement of molecules in a 
crystal, when excited by free excitons. 

Chapter V investigates the solutions of equations obtained in the 
zeroth-order approximation. It is shown that the zeroth-order approxi- 
mation coincides with the case of rigidly fixed molecules. The limits 
within which the results obtained in Chapters II and III are applicable 
are thereby shown. 

Chapter VI investigates the solutions of the same equations in the 
first approximation. This enables us to describe the processes of exci- 
tation of free excitons simultaneously with the excitation of lattice vibra- 
tions. Also examined are the problems of the structure of an absorption 
band and its dependence on temperature. 

Chapter VII investigates the process of light absorption when “‘local- 
ized” excitons areformed. A comparison is made of the excitation energy 
upon formation of “localized’’ and free excitons. 

Chapter VIII investigates the behavior of free and “localized” excitons 
subsequent to their formation. The necessary conditions are established 
for the appearance of fluorescence in solid bodies that do not contain 
admixtures. 

Chapter IX investigates the application of the theory to several experi- 
mental data in the field of combinational scattering (Raman) spectra and 
infrared spectra, as well as absorption and fluorescence spectra in the 
visible and ultraviolet regions. 

Chapter X applies exciton concepts to the calculation of excited energy 
states of certain individual polyatomic molecules. 


CHAPTER I 


TRANSLATIONAL AND ROTATIONAL VIBRATIONS OF 
MOLECULES IN A MOLECULAR LATTICE 


4. Atomic Lattices 


When investigating the movement of molecules in a molecular lattice, 
it is convenient in the beginning to consider the lattice as an atomic one 
and, subsequently, to take into account the fact that in a molecular 
lattice the atoms group themselves into molecules. The interaction 
forces between these molecules are considerably smaller than the inter- 
action forces between the atoms which enter into the make-up of one 
molecule. The possibility of considering a lattice as such when investi- 
gating the movements of molecules is based on the fact that the displace- 
ments of atomic nuclei are essential for such movements. Considering 
the lattice in this manner is at times inadmissible when investigating the 
optical properties of molecular crystals in the visible and ultraviolet 
regions, the magnetic properties, and other properties. 

Let us consider a crystal whose unit cell is determined by three non- 
coplanar vectors, ai, a2, and a3, and contains ¢ molecules, each of which 
has patoms. Let us assume that the crystal has the form of an oblique- 
angled parallelepiped with axes Nya;, N2a2, and N3as, so that the general 
number of unit cells in the crystal is equal to N = NiN2N; (N assumed 
very large). In order to avoid the difficulties of considering the super- 
ficial layers of the crystal, we shall take as limiting conditions a cyclic 
system with large periods that correspond to the infinite repetition of 
the crystal in all directions. 

Let us designate by R2* the z component of the displacement of the 
atom which occupies the ath position in the mth cell; herea = 1,2, .. 
prandn = mn = nya, + Node + 33 is the vector of the direct lattice. 

In the decomposition of energy in powers of displacement from the 
equilibrium positions, if one limits oneself to quadratic terms, then the 
total energy of small vibrations about the equilibrium positions will have 
the form! 


Sy 


E=4Y (m\Re|? + Y xgege’ RerRze’). (1.1) 
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In this circumstance, the equations of motion become 


mals + YS nznzie’ Ree’ = 0. (1.2) 
nats’ 

Here m, is the mass of the atom, occupying the ath site in the cell; the 
coefficients 2:7,’ for the ideal lattice depend only on the difference 
n—n’, not onnand n’ separately. These coefficients form a difference 
matrix in relation to the indices n and n’. Each element of this matrix 
is, in turn, in ultimate arrangement the matrix relative to the indices x, a. 
Owing to the translational symmetry of the crystal, the equations (1.2) 

have the solutions 
Rz = e*#(q)eian-et) (1.8) 


We wrote (1.3) in the complex form. However, the displacements of 
atoms in the lattice are physically real magnitudes; therefore, one must 
take the real part of (1.3). The expressions (1.3) represent a running 
wave with a wave vector q. 

In order to satisfy the cyclical conditions, the wave vector q must have 


the discrete values 
3 
1=| 
Qa 
a= DRG on (1.4) 


where 7: are whole numbers which satisfy the inequalities 


N; N; 


—y <n Rp t= 1, 2, 3, 


and (a,)—! are the vectors of the reciprocal lattice 


[a,ax] 


(a:)-! = a;[acas| * 


(1.5) 


The vectors e*(q) characterize the vibrational directions in the running 
wave; they are defined as the solutions of the equations obtained after 
substitution of (1.3) into (1.2): 


—wme + Yo pzargia’eia(aaer'a’ = 0, (1.6) 
n'a! 
Let us introduce the designation 
Lrae'e'(q) =) azege’etataa), (1.7) 
o 


L7«2'e'(q) is a continuous function of the wave number q; this function 
does not depend on n, the indices of the cells. 


10 DAVYDOV: THEORY OF MOLECULAR EXCITONS 


The matrix \2*:’ is symmetrical in the indices n and xa; that is, 


Nisaa ena (1.8) 
Furthermore, it is also symmetrical in the indices 2; that is, 
Neeser ze 2 (1.8a) 
From the symmetry properties [(1.8)] it follows that 
Lroa'a'(q) = Lha'e'9(q), (1.9) 
[t4,2'a! (q) = D*ta.z'a'( gq), (1.10) 


For those crystals in which each atom is the center of symmetry, the 
coefficients L**’e’(q) will be real. 

With the help of (1.7), the equations of motion [(1.6)] may be written 
in the form 


> [702 (q) ere’ aed w?m,e7% = 0. (1.11) 


The squares of the eigenfrequencies w? are the roots of the equations of 
Sopth order 


[Le*#’"(q) — w(q)madza,z/a’| = 0. (1.12) 
Owing to the Hermitian nature [(1.9)] of the coefficients of this equa- 
tion, the roots of the latter will be real and positive (they cannot be 


negative since the equilibrium positions correspond to an energy mini- 
mum) continuous functions of the wave number q: 

w? = w3(q), j=1,2,... , 3c. (1.18) 
As is known, equation (1.12) has three roots, equal to zero when q = 0; 
these solutions correspond to the acoustical class of vibrations. The 
remaining 3(p7 — 1) roots of (1.12) are different from zero when q = 0; 
they correspond to the optical class of vibrations. 

Each root [(1.13)] of equation (1.12) will correspond to its own sum 
total of solutions of the system of equations (1.11), which we shall dis- 
tinguish by means of the index j; the latter runs through the values 
1,2, ... ,8p0: 

eF*(q). (1.14) 


The vectors [(1.14)] form orthogonal systems, which may be normalized 
in the following manner: 


x MaeF er" = opdy;. (1.15) 


za 


The arbitrary displacement of a desired atom of the lattice may be 
decomposed into a series according to solutions obtained: 


Rat = (Noe)—* Y) eF*(q) (aie + aye). (1.16) 
jg 
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The coefficients a,; = aj,e characterize the amplitude of the jth 
class, which has a wave vector q and a vibrational direction that is 
determined by the index 7. For R2* to become physically real it is 
necessary to fulfill the conditions 


a, = ag. (1.17) 


The summation over q in (1.16) is carried out only for gz > 0. Substi- 
tuting (1.16) into (1.1), and taking into consideration that 


» effa-a2 = NG qa’, 
n 


we obtain the total energy for lattice vibrations in the following form: 
E=2) wta,a%.. (1.18) 
a 
From the complex variables a,; and a%,, one can go over to the real 
canonical variables x,; and P,; with the aid of the transformation 


tq = — (2 + ize) (1.19) 


Then the total vibrational energy will be reduced to a sum of quadratic 
terms: 


EB = '6) (P+ 03,03)). (1.20) 
a7 


The transition from classical to quantum mechanics is accomplished in 
the usual manner. It is necessary to consider the canonical variables of 
each vibration as noncommuting quantities which satisfy the commuta- 
tion relationships 

Pa j%qijt — Lqt7tP gs = — th’ 8 qq biz". (1.21) 

The quantum-mechanical equation which describes the movement of 
atoms in the lattice is reduced to the equations for a certain number 
(83Npc) of independent oscillators, which are distinguished from each 
other by the quantities g and 7. The wave function is decomposed into 
the product 

Il W(2q;), 

q 
where W(2,;) satisfies the equation 

BY (x 2M 4; M jw? 2? 
As ) a ie (z = . gi ) we) = 0. 


The solution of this equation has the well-known form 


e—l2?/ 2299] H x 
Vy, ;(2) = A/ tn Ngj (2 ) 
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where N,; represents whole numbers and where 


x = Rae 
~T M 30295 


Ex; = hegj(Nq; + 34). 


The energy levels are 


Consequently, the total energy of the crystal is 
E = J (Ng + ¥4)eash, (1.22) 
QI 


and the wave function for the corresponding state is 


IT Wn; (Lei)- 
qi 
The whole numbers N,; determine the number of phonons which have 
wave number g and polarization 7. 
If we consider the canonical variables x,; as independent, we can repre- 
sent them by means of Hermitian matrices: 


(Ng; fle 1) e 


= 7* = 
UNoiNgitl SN t1.Ng; J 209i 


TNqiNyi’ = 0, if NY ¥ Nai + 1. 


tig ;t 
, 


Consequently, z,; has matrix elements different from zero only for such 
transitions in which N,; decreases or increases by unity. 

However, because they are independent variables, it is more convenient 
to consider the complex quantities a,;, which can be represented by non- 
Hermitian matrices: 


h(N 4; + 1) ; 


ONy.Ngtl = Da tage (1.23) 
ANG, 
OF arya = af Be eet (1.24) 


These matrices are proportional to one and the same exponential factor 
when the matrix elements xz,; contain various time-dependent factors. 
The quantities a,; satisfy the commutation relations 


h 


a (1.25) 


* Ky 
AgiAg; —~ Aq Aq; = 


§. Molecular Lattices. Internal and External (Lattice) Vibrations 


A distinguishing characteristic of molecular crystals is that in such 
crystals the molecules preserve their individuality in the first approxi- 
mation. This is because the interaction forces between molecules are 
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small in comparison with the forces acting between the atoms of a single 
molecule. Therefore, all possible atomic vibrations in molecular crystals 
may be divided into two groups: 

1. Internal vibrations of the atoms of a molecule relative to one another. 
In these vibrations, the center of gravity of the molecule is not displaced 
and there is no rotation of the molecule as a whole. As a result of the 
strong bonds between atoms in the molecule, these vibrations possess 
relatively large frequencies (800 to 3000 cm—'). Temperature and other 
external agents exert little influence on these eigenfrequencies; therefore, 
these same frequencies (with small variations which will be investigated 
later) are preserved in various states of aggregation. In view of the 
small connection between the vibrations occurring in neighboring mole- 
cules, the frequencies of these vibrations are but little dependent on the 
wave number g. 

2. External vibrations, or lattice vibrations, which appear because of 
rotational and translational degrees of freedom of the molecule. For the 
majority of molecular crystals, the frequencies of these vibrations rarely 
exceed 100 cm-!. These vibrations are not encountered in vapors and 
solutions, and therefore, in the future, we shall simply call them lattice 
vibrations. The frequencies of these vibrations depend essentially on 
the wave vector q. 

In reality both types of vibration are lattice vibrations and cannot, 
strictly speaking, be separated from each other. However, in molecular 
lattices the separation described above is entirely possible in the first 
approximation and results in a significant simplification when investi- 
gating vibrations. 

At times, the influence of the crystalline lattice on the interna] vibra- 
tions of the molecule leads to a change of selection rules and of polari- 
zation, as well as to a small displacement and splitting of the frequencies 
of internal vibrations. All these changes are relatively small (especially 
for vibrations with no change in dipole moment) and will be investigated 
in Chapter II. In the meantime, we shall consider that the internal 
vibrations of the molecules in the crystal are known. 

As we saw in the preceding paragraph, there are 3p classes of vibra- 
tions, that is, 30p types of vibrations for each wave vector, in a molecular 
crystal which contains in the unit cell «» molecules, each having p atoms. 
Three classes of vibrations are acoustical, the remainder optical. 

Out of all 3cp classes of vibrations of the molecular crystals, 6c classes 
belong to the lattice vibrations; out of these, 6 — 3 are optical classes 
and 3 are acoustical. Out of 6¢ — 3 optical classes of molecular vibra- 
tions 3(¢ — 1) classes are translational vibrations and 3c classes are 
rotational. 

All remaining (8p — 6)o classes will be determined by the internal 
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vibrations of the atoms in the « molecules of the unit cell; they all belong 
to optical classes. 
Figure 1 shows the dependence on the wave vector g of the various 
classes of vibrations in the molecular crystal. 
Hereafter the term lattice vibration will mean only external vibrations. 
In other words, out of all classes of vibrations considered in Section 4, 
we shall concern ourselves only with the first 


60 classes (including the 3 acoustical ones). 
We can utilize all of the results obtained in 
Section 4, considering that the index j, which 


determines the type of vibration, runs only 
through the values 


o(3p—6) 


[-~-----—- j=1,2,..., 60. (1.26) 


Internal vibrations 


6. Classification of Lattice Vibrations 
According to Irreducible Representations 
of the Symmetry Group of a Crystal 


3(ap-1) 


| 
H 
l 
{ 
{ 
| 
! 
| 
\ 
Optical branches 


At the present time it is hardly possible to 
calculate the external vibrations of a com- 
plicated polyatomic lattice. Therefore, we 
must utilize experimental data. 

Itis comparatively simple to classify lattice 
vibrations according to the irreducible repre- 
sentations of the symmetry group of the 
crystals and to establish the laws of behavior 
Fic. 1. Dependence of vibra- of these frequencies in various types of spec- 
ene aedaney one wave tra [infrared and combinational scatterin 
vector for various classes of & 
vibrations. (Raman) spectra, etc.] and their polarizations 

as well. Furthermore, the lattice vibrations 
can be divided into rotational and translational ones. The vibrations of 
the translational type usually have small frequencies (especially in crys- 
tals of low melting point, i.e., those with weak binding forces between 
molecules). 

Such a qualitative examination can be carried out by the methods of 
group theory and is extremely useful when analyzing experimental data. 

By way of an example, let us consider molecular crystals having two 
molecules in a unit cell. Such crystals usually crystallize in the mono- 
clinic system and belong to the space symmetry group C,. 

The symmetry group C$, has four symmetry elements: 

1. The identity element EF. 

2. Rotation by 180° about the monoclinic axis C,. In this operation, 
the molecules in the unit cell change places. 


Lattice vibrations 
60 
—— 
Yev 
Acoustical 
branches 


q--- 
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3. Inversion 7. Both molecules of the unit cell remain in their places. 

4. Reflection o in a plane perpendicular to the monoclinic axis. In 
this operation, both molecules of the unit cell exchange places. 

The characters of the irreducible representations of the group C%,, 
showing the transformation properties of the radius vector components 
Te, To, aud r,, are given in Table 1. The table also shows the characters 
of the reducible representations: x(T), the characters of all translational 
motions of the two molecules in the unit cell; x(T.), the characters of 
translational motions corresponding to the acoustical classes; x(R), the 
characters of rotational motions of molecules in the lattice. 


TaBLE 1. CHARACTER TaBLE FoR THE C3, Space Group; JRREDUCIBLE AND 
REDUCIBLE REPRESENTATIONS 


Cc E C t o 


5 
2h 2 
A, 1 1 1 
rr, Au 1 1 ~1 -1 
B, 1 -1 1 -1 
Ta, Te, Bu 1 -—l -l 1 
x(T) 6 0 —6 0 
x(T a) 3 -1 —3 l 
x(R) 6 0 6 0 


Decomposing the characters x(T) of the reducible representation I(T) 
for translation into the irreducible representations of the group C3,, we 
obtain 

I(T) = 3A, + 3B, (1.27) 


In order to obtain only the relative motions of both molecules belong- 
ing to the optical class, it is often necessary to separate from all possible 
translations the translations of the lattice asa whole. (The lattice trans- 
lations correspond to the extreme acoustical vibrations.) If the reduci- 
ble representation of the translations of the lattice as a whole is desig- 
nated by I'(T.), then 

T(T.) = Au + 2By,. (1.28) 


Subtracting (1.28) from (1.27), we obtain the representation T'(7',) of the 
optical class of frequencies corresponding to the progressive shifts of the 
molecules relative to each other: 


I(T) = 24, + Bu. (1.29) 


Utilizing the characters of the reducible representations T(R) of the 
molecular rotations in the lattice, we obtain the decomposition of this 
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representation into irreducible components: 
T(R) = 8A, + 3B,. (1.30) 


In this manner, the six possible rotational vibrations in the lattice are 
divided into two groups: three vibrations belong to the totally symmetri- 
cal representation A, and three vibrations belong to the irreducible repre- 
sentation B,. Asa result of the small intermolecular forces in molecular 
erystals, the frequencies corresponding to the vibrations of symmetry A, 
and B, are little distinguished from each other. 

The presence of lattice vibrations in the combinational scattering 
(Raman) spectra was first discovered by E. I’. Gross and M. F. Vuks.? 
There are now many works covering the investigation of these vibrations. 
Among these we may mention the works of Gross and collaborators, ***4 
of Vuks,°* of Raskin,®* of Nedungadi,*” and of Rousset.*® The frequen- 
cies of crystal vibrations in single crystals of naphthalene were especially 
thoroughly studied in reference to combinational scattering (Raman 
scattering). Six vibrations, distributed in three pairs of similar frequen- 
cies (38, 47; 72, 78; and 107, 121 cm™), were found. All these frequen- 
cies belong to rotational vibrations of the molecules. The presence of 
approximately the same pairs of rotational-vibrational frequencies (54, 
46; 76, 74; and 109, 127 cm—!) was confirmed in the more recent work of 
Rousset.# 


CHAPTER II 


THE INFLUENCE OF THE INTERACTION FORCES 
BETWEEN MOLECULES OF A MOLECULAR 
LATTICE ON THE INTRAMOLECULAR 
VIBRATIONS + 


7. Free Vibrations 


In this chapter we shall examine the approximation method of account- 
ing for the influence of the crystal lattice on the internal vibrations of 
molecules. This problem is of great interest (1) because its rigorous 
solution by the method set forth in Chapter I is infeasible in practice 
even where the interaction forces between the atoms making up a mole- 
cule and between the molecules are known and (2) because very often 
investigators are specifically concerned with the problem of the change in 
molecular vibrational frequencies (calculated or experimentally deter- 
mined when investigating gases at low pressures) when molecules enter 
into a crystal lattice. 

The theory of such changes must relate them to the interaction forces 
of molecules in the lattice and with the structure of the crystal. This, 
on the basis of the study of spectra [in particular combinational scatter- 
ing (Raman) and infrared], will enable us to deduce certain features of 
the crystal structure. It is evident that spectral methods of investigating 
the crystal structure are extremely useful for verifying and supplement- 
ing results obtained by other means. 

It is known that, when investigating the vibrations of atoms relative 
to each other, we can, with great approximation, utilize the laws of classi- 
cal mechanics, which are simpler and more obvious. In the following 
chapters we shall give the more rigorous quantum-mechanical theory and 
show the limits within which the approximation results obtained here are 
applicable. 

Thus we shall proceed from the assumption that we know the frequen- 
cies and symmetry properties of the normal vibrations of each molecule 
taken separately. Let us consider the changes of the normal vibrations 
of ¢N molecules when they are united into a molecular crystal containing 


t A short summary of the principles in this chapter was published by the author.*? 
17 


18 DAVYDOV: THEORY OF MOLECULAR EXCITONS 


o molecules in the unit cell and when they are rigidly fixed in equilibrium 
positions corresponding to the equilibrium state of the crystal. 

If each molecule possesses p atoms, the molecule will have 3p — 6 inter- 
nal normal vibrations. A crystal containing oN molecules will have 
oN (8p — 6) ‘internal’ vibrations (the molecules being fixed). In the 
zeroth approximation, such a crystal may be considered the aggregate 
of completely noninteracting molecules. In this case, all oN(8p — 6) 
vibrations of the crystal break down into 3p — 6 different oN multiply 
degenerate vibrations, if, that is, the normal vibrations of the separated 
molecule do not include degenerate vibrations. In the converse case, 
the number of different vibrations of the crystal in the zeroth approxi- 
mation decreases, and their multiplicity increases proportionately. 

In reality, however, the molecules of a crystal interact with each other, 
and consequently, the vibrations of separated molecules cannot be con- 
sidered as independent. As we have already shown, the problem is 
reduced to determining the changes introduced by the interaction of 
normal molecular vibrations. 

By taking into account the specific resonance character of vibrational 
interaction, the solution of this problem can be significantly simplified. 
By way of example, let us consider the interaction of two normal vibra- 
tions having frequencies wio and wa and an interaction energy of the 
dipole type yaia2, where a, and a» are coordinates of the corresponding 
normal vibrations. The equations of motion of such a system will have 
the form 


Gy + wea, + 7 an = 0, 
(2.1) 
do + we ote + way = 0. 


Solving this system of equations by substituting 
a; = azet, 
we obtain two values for the frequency w: 
2 
ola = [ute tate + afloat +4 (2) | 22) 
If 19 = Wea, then 


(2.3) 


Consequently, if the normal frequencies coincide (resonance), the correc- 
tion to the square of the frequency will be of the first order relative to 
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the bond coefficient y/m. Therefore, if 


|wty — w3o| > a (2.4) 
then 
2 2 (y/m)? 2 2 (y/m)* Z 9 
oF “ Of) — 50’ oe oe Win — Wo ( =) 


Consequently, if the normal frequencies do not coincide, the correc- 
tions to the frequencies will be quadratic relative to the bond coefficient 
y/m. Therefore, if the inequality [(2.4)] is fulfilled, it is possible not to 
take into account the changes of frequencies caused by the interaction of 
vibrations. 

We shall assume that the following inequality is fulfilled: 


Jo3 — w?] > 2% i=1,2,...,3p-6 (2.6) 


and consider oN vibrations of the crystal, corresponding to the frequency 
wo in the zeroth approximation. Let us designate by an. the normal 
coordinate of the molecular vibration in the molecule occupying the ath 
site of the nth cell. Then in the isolated molecule, this coordinate will 
satisfy the equation 

Gna + WOna = 0 (2.7) 


where wo is the normal vibrational frequency (vibrational frequency of 
the free molecule). 

Initially we shall consider the case in which a change of the electric 
dipole moment of the molecule corresponds to the vibration with fre- 
quency wo. In this case the interaction between molecules will have a 
dipole-dipole character. The total energy of oN vibrations will be 


E = ca) {), @. + v3 ) ake 
+ YR mplPraPme — 3(Pmitname)(Profnaims)|}, (2.8) 


na,mp 


where Phe = @ Vf Cadna (2.84) 


is the effective dipole moment of the molecule na; ea is the unit vector 
of the direction of the dipole moment of the molecule occupying the site 
a in the unit cell and is determined by the type of vibrations in the 
molecule and by the orientation of the molecule relative to the crystallo- 
graphic axis of the crystal (the vector, here, assumed to be known); ¢ is 
the electronic charge; froma is the unit vector in the direction of the line 
connecting the molecules na and m8; and f is the oscillator strength, 
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which here, in accord with Ladenburg,*° is introduced as a certain factor, 
allowing for the correction of the quantum behavior of the oscillator. 

The equations of motion of the system of oscillators with energy given 
by (2.8) will be 


Gna + wena + > Yna,mbAmél (Cas) = 3 (€p¥na,mB) (Cal na,ms)] = 0, (2.9) 
Hacend 
fe? 

where Yuams = 553,-—" 2.9a 
2 Rae ms ( ) 
The solutions of this system of equations may be looked for in the form 

of running waves: 
Gna = Age ea—ot) (2.10) 


where @z is the amplitude of the vibrations and k is the wave vector run- 
ning through N discrete values as a result of the cyclic conditions (with 
large period). We assume that the crystal, in the form of an oblique- 
angled parallelepiped with axes Niai, Noae, and N2as, repeats itself end- 
lessly in all directions (V = NiN.2N; and is very large). Then 


k = y az ni, (2.11) 


where 7: is a whole number satisfying the inequality 


N: N: 3 F 
—-a<ugy t= 1,23, (2.12) 
and the az! are the vectors of the reciprocal lattice. 
Substituting (2.10) in (2.9), we obtain 
(w3 — w?)aa + Y Tes(k)ag = 0, (2.13) 
a= 


where 
Tae(k) = bY Yne,mpe'*™—| (@e@p) me 3 (CaF na,mp) (€pTna,ms)]. (2.14) 


If a ¥ 8, then the summation over m includes m = n; if a = 8, then the 
summation includes only m ¥ n. 
The matrix [(2.14)] is Hermitian: 


Tag(k) = 3. (k). 


However, for those crystals in which each molecule is a center of sym- 
metry, this matrix is real and symmetric relative to the indices a and 8. 
The system of homogeneous equations [(2.13)] has nontrivial solutions if 
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the determinant composed of the coefficients of these equations is equal 
to zero: 
A(w?) = |Tas(k) + (wi — w?) daa] = 0. (2.15) 


The preceding equation determines ¢ roots of w?, as functions of the 
wave vector k: 
wz = wi(k), w=1,2,...,¢. (2.16) 


In this manner, as a result of the interaction between molecules, instead 
of one degenerate frequency wo of the vibrations, one obtains, generally 
speaking, » different (some of them may coincide) quasi-continuous bands 
of frequencies (corresponding to the number of molecules in the unit cell) 
or classes of frequencies. Each of these bands of frequencies will con- 
sist of N closely spaced frequencies distinguished by various values of the 
wave number k; these values satisfy the conditions given by (2.11). 

For the sake of greater concreteness, let us examine, by way of example, 
crystals having two molecules per unit cell. In this case, the equation 
(2.15) takes the form 


2 2 
ws — w? + Tii(k) Ti2(k) : 
T21(k) we — w? + Ti1(k) Fy, 0. (2.17) 
Its solution is 
w(K) = w8 + Dinde) + Tir(k), nae 


w? (kk) —_ ws + T(k) = Ty0(k). 


We obtain, with the aid of (2.13), two systems of coefficients a. corre- 
sponding to the two roots of (2.18): 


He ay a, 

(2.184) 
1 1 
2 1 —1 


These determine the solutions [(2.10)]. 


8. Forced Vibrations. Absorption in the Infrared Region Corresponding 
to the Excitation of Intramolecular Vibrations in Crystals 


Let us now consider the forced vibrations arising in a crystal under 
the influence of a plane electromagnetic wave: 


E = Eyei(or-00), (2.19) 


where Q is the wave vector of the electromagnetic wave and w is its 
frequency. 
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To clarify the optical behavior of a crystal, it is necessary to deter- 
mine the electric moment of a crystal: 


P= Pye, 


the electric moment induced by the electromagnetic wave [(2.19)]. This 
moment is expressed in terms of the dipole moment P.. of separate mole- 
cules in the following way (all molecules being identical): 


P, =e Vif» a,ene", (2.20) 


where aj, satisfies a system of nonhomogeneous equations 


ov 


(w? — w)al + bY Taa(k)ah = Mz, (2.21) 
f=1 
a 74 (Q—k)n 
where M, = eV} Eves > : aa (2.214) 


n 


Here, owing to the fact that the wavelength of light is great in com- 
parison with the distance between molecules, we disregard the difference 
of the phases of the light wave’s action on the molecules occurring in a 
unit cell. 

Particular solutions of the equations (2.21) may be represented in the 
form 


>) AseMs 
SS ) (2 = 
A(w?) ° 


a= 


where A.s is the minor of the determinant A(w*), corresponding to the 
element 
Task) + (wp — w?) dag. 
Since 
5 eka = 5(Q — 10, 
it becomes immediately apparent that, under the influence of the electro- 
magnetic wave, only the eigenfrequencies 


wy = w4(Q) (2.23) 


may be excited. These frequencies correspond to the value of the wave 
vector k = Q, since only for these values of the squares of the frequen- 
cies is the numerator of (2.22) different from zero and the denominator 
equalto zero. Thus the values of the frequencies are the resonance points 
of our system. At the resonance points, the amplitude [(2.22)] is infinite. 
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This result is obtained because the model we are considering is too 
idealized. We do not take into account the transformation of energy 
to the lattice vibrations (see Section 11). This transition of energy is 
equivalent to the damping of the vibrations we have considered. The 
damping of vibrations causes the amplitude to remain finite and within 
the region of resonance. 

Because the wavelength of infrared radiation is very great, the wave 
vector Q is small, and we can set Q = 0. Thus we come to the con- 
clusion that under the influence of an impinging electromagnetic wave, 
only limited wave frequencies may be excited in each band of frequencies 
[(2.16)]. Substituting in (2.20) the solution obtained in (2.22), we may 
determine the electric moment arising in the system: 


P, = Posem*e", (2.24) 
where 
¥ AasMs 
ease Bian teak 
Py, =evfN yy ee ae) (2.25) 
w=1,2,...,¢ 


In particular, for crystals containing two molecules per unit cell, we obtain 


2 
P, = Poe", Poi = GINT (e1 + 2), (2.26) 
A(a}) 
ee 2 MINT ae 2.2 
Py = Pore, Poo = Alo’) (e: — @2). (2.27) 


Consequently, for the frequency w; determined by (2.18) for k = 0, 
the electric moment of the system is proportional to the sum of the elec- 
tric moments of the two molecules occurring in the unit cell; but for the 
frequency we, the electric moment is proportional to their difference. 


our 


| 
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Fic. 2. Projection of electric moments of two molecules (left) and resultant moments 
of entire system (right) onto ab plane of monoclinic crystal, with two molecules per 
unit cell. 


Figure 2 shows the possible position of the projection onto the plane ab 
of a monoclinic crystal of the electric moments of the two molecules e; 
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and @, of the unit cell of the crystal, as well as the projection onto the 
same plane of the direction of the vectors of the resultant electric moment 
of the entire system. Thus, because of the resonance interaction of the 
intramolecular vibrations, the infrared absorption and emission bands 
must double in comparison with the corresponding bands for the gas, 
and they must be strongly polarized. The ratio of the intensities of the 
absorption bands will be 


I(w1) = (e1 + €2)? 
I(we) = (@1 — @2)? 


If such an interaction did not exist, only one frequency would be absorbed, 
and the strong polarization would be absent for e; not parallel to es. 
It is interesting to note that if the dipole moments of the normal molecu- 
lar vibrations of two molecules in a unit cell are parallel to each other, 
then the electric moment of the system will differ from zero only for 
(2.26). Here the polarization of the absorption bands of the crystal will 
coincide with the polarization of the absorption bands of the free mole- 
cule, and the position of the absorption bands of the crystal will be dis- 
placed relative to their position in the spectrum of the free molecule. 

Investigations of the infrared absorption of molecular crystals in polar- 
ized light are extremely complex and comparatively crude. In order to 
detect the effect described above, one must conduct experiments at the 
lowest temperatures possible and with the smallest possible thicknesses 
of the absorbing substance. Nevertheless, in the experiments of Elliott 
and Ambrose,‘! there was observed a noticeable pleochroism (at room 
temperature) of the infrared absorption of naphthalene crystals when the 
wave fell perpendicular to the plane ab of the crystal (for the cases in 
which the direction of the electric vector of the wave is parallel to the 
axes b and a of the crystal) in the region 2 to 4 uw, corresponding to the 
intramolecular vibrations of the crystal. 


(2.28) 


9. Intramolecular Crystal Vibrations in Which the 2/-pole Moment of 
the Molecule Changes 


In Sections 7 and 8 we considered the interaction between molecules 
that occurs upon excitation of intramolecular vibrations of the dipole 
type. We shall now consider the general case in which, to the normal 
frequency wo, there corresponds the formation of a 2/-pole moment of a 
molecule (for 2 = 1, a dipole moment; for ] = 2, a quadrupole, etc.). 
Let us designate the corresponding moment of the molecule occupying 
the ath place in the nth cell by 


P2 


ra 


20 


= HOna€yray 


where e?4 is a unit tensor of the corresponding rank characterizing the 
transformation properties of a given multipole of the molecule, and x is 
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the coefficient of proportionality between the multipole moment and the 
value of the normal coordinate ane of the intramolecular vibration. Since 
we consider the molecules rigidly fixed in the lattice, the axes of the tensor 
will have a constant orientation relative to the crystallographic axes of 
the crystal, the only change being that of the normal coordinate ana, 
of the intramolecular vibration, characterizing the magnitude of the 
corresponding multipole. 

Let D2, .@ determine the interaction between two 2/-pole moments of 
molecules na and m@. It is essential that D%., be inversely propor- 


na,mp 


tional to the (22 + 1) power of the distance between molecules: 


D2, yg ~~ ROO) (2.29) 


na,mp 


For the case in which | = 1 (dipole-dipole interaction), 
DE cis = Yna,mbl (Calg) a 3 (Cal na,mp) (CT nama). (2.29a) 


The magnitudes entering into (2.29a) are determined in Section 8. 

We must, however, keep in view that besides the interaction between 
molecules [(2.29)] (occurring at the moment of excitation of the mole- 
cules), there is an interaction even before excitation. It is precisely this 
interaction (van der Waals forces and exchange repulsion forces) which 
determines the configuration of molecules in the crystal in its normal 
state. Upon excitation of intramolecular vibrations, this interaction also 
changes; therefore, in a more rigorous theory, it must be taken into 
account. Let us designate this type of interaction energy by the letter V. 
We assume that this energy depends on, besides the distances and the 
mutual orientations of the molecules (here fixed), the vibrational state, 
ie., on the normal intramolecular coordinates azo. When the normal 
coordinate Qn. changes, the force 


aV 


Pra ees dana 


will act upon the molecule ne. If the displacements from the equilibrium 
positions are small, this force may be decomposed into a power series of 
normal coordinates: 


aVv aeV 
Fao = — (20) (er) et. 280 
mp 


The first, constant, term of the force [(2.30)] will result in a small dis- 
placement of the equilibrium position for the intramolecular vibrations. 
The second term in (2.30) is, generally speaking, also small and falls off 
with distance not slower than R-*. It is therefore expedient to take 
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(2.30) into account only when the interaction between excited molecules 
is not lower than quadrupole-quadrupole. 

When taking (2.30) and the 21 — 2i-pole character of the interaction 
between excited molecules into account, the equations of the vibrations 
of No molecules may be written in the form of the following system of 
nonhomogeneous equations: 


Gna + 2200 + ) M4, gdm = — (=) (2.31) 
mp Tee, 
mprtna 
where 
ov 
02 = we + (<r), (2.31a) 
a7V 
2 = pr _\. 2 
Me. mB D2 mB Gi =) (2.31b) 


The solution of the system of equations [(2.31)] can be searched for in the 
form 
One = Oe, + aac hao), (2.32) 


where the wave vector k runs through the same N discrete values as in 
(2.11). 

The magnitudes az, determine the constant displacements; they are 
calculated with the aid of a system of nonhomogeneous equations: 


Ma, + ) MY, ppg = — (2). (2.33) 
0 
mp 


OGna 


These magnitudes are, generally speaking, very small and will not 
interest us. 
The amplitudes a, now satisfy the system of homogeneous equations 


edz + Y V25(k)as = 0, (2.34) 
ga 
where 
e = Q2 — o?, (2.344) 
Ti (k) =) Oe ee (2.340) 


Here the summation over m has the same meaning as in (2.14). 
From the conditions of solubility of the system of equations [(2.34)], 


“A*(e) = [Ta5(K) — €da8| = 0, (2.35) 


we shall obtain the algebraic equation of the oth degree relative to «. 
The roots of this equation will give o functions, 


& = e(k), 
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and consequently, also « quasi-continuous bands of frequencies: 


2 
ot = oh + (F4) Ee: GeBee so ke 1036) 
na f 0 
Thus, in the crystal, splitting and displacement of the frequencies of 
normal molecular vibrations take place. In the particular case of crystals 
with two molecules per unit cell, 


2 = 2 3° 

w} (k) = Wh + aa2 , + T1,(k) + T y2(k), 
2 ans 2 3? — 

we(k) = w? + daz.) + Ti1(k) T12(k). 


Thus the splitting in the crystal of the molecular frequency wo into 
o quasi-continuous bands of frequencies [(2.36)] applies generally to all 
types of vibrations of the molecule. However, the distance between 
bands and the width of the bands of frequencies, which is determined 
by the magnitudes of 1'24(k), will decrease with the increase in the vibra- 
tional multipole index. This was to be expected, since the resonance 
interaction between molecules falls off rapidly with an increase of the 
multipole index of the vibrations. 

The investigation which has been conducted permits us to draw a 
rather interesting conclusion. In a crystal, splitting of the molecular 
frequency wo into « quasi-continuous bands of frequencies occurs, and 
the polarization character changes essentially. The magnitudes of the 
splittings and the widths of the bands of vibrational frequencies are sig- 
nificantly larger when accompanied by a change in the molecular dipole 
moment (in which case they appear in the infrared spectrum) than are 
the corresponding magnitudes for those vibrations which do not result in 
changes in the molecular electric moment. 


10. Influence of the Interaction Forces between Molecules on 
Combinational Scattering (Raman Scattering) 


As is well known, spectra of combinational scattering (Raman scatter- 
ing) of light play an important role in the spectral study of the structure 
of matter. The value of this method consists in the possibility of investi- 
gating intramolecular vibrations in the visible and ultraviolet regions of 
the spectrum but not in the experimentally difficult infrared region. 
Recently, thanks to the discovery by E. Gross and M. Vuks? of the low- 
frequency spectrum, one can also utilize combinational scattering (Raman 
scattering) to elucidate the structure of molecular crystals and liquids. 

Naturally, it is of great interest to investigate the change in the combi- 
national scattering (Raman) spectrum upon transformation of the mole- 
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cules from a gaseous to a liquid or solid state. Repeated efforts were 
made to take into account these changes by means of examining the influ- 
ence on the spectrum of the field created by the van der Waals forces. 
One of the first theoretical studies in this direction was that of Buch- 
heim.42, This study was not without basic defects and has already 
attracted criticism.** 

We must accept as significantly more successful the effort undertaken 
in this connection by M. V. Volkenstein.‘? However, even Volkenstein’s 
calculation is not without a series of defects, even within the framework 
of a classical model. 

When investigating intermolecular interaction, Volkenstein considered 
the interaction of only one pair of molecules. Since, generally speaking, 
the changes of intramolecular vibrations resulting from interaction with 
surrounding molecules are not additive, such consideration was insuf- 
ficient. Furthermore, Volkenstein did not consider the forces of repul- 
sion between molecules and took into account only the van der Waals 
attraction for the energy, for which he used the expression 


1/2, 
— (ser Diet iieed diae $4Taxas)) 


where uw: and yw, are the molecular dipole moments, a; and ae are the 
polarizabilities, J is the ionization potential, and T is the temperature. 

The first term of this formula has meaning only for gas molecules, 
since it is obtained by averaging the interaction of the dipoles over all 
possible orientations. This term is not applicable to liquids and is even 
less applicable to solids. Its use may be justified only at sufficiently high 
temperatures, but then the term is small; as a result of heat fluctuations 
the change of frequencies is mainly determined by the change of inter- 
molecular distances, which, in the calculation, are assumed to be fixed. 
Since the main change of combinational scattering (Raman) frequencies 
is, according to Volkenstein, determined by precisely the orientational 
interaction, one must be wary of the results obtained. 

A much more essential observation is that, when calculating the influ- 
ence of the interaction forces between molecules on the intramolecular 
vibrations, one must take into account the interaction forces between 
excited vibrations. Here the so-called resonance forces of interaction 
between molecules occur. When the molecules are excited with a con- 
comitant change in dipole moment (i.e., vibrations active in the infrared 
spectra), these resonance forces vary only as the inverse fourth power 
(and not the seventh, as do the van der Waals forces) and, therefore, 
surpass by a hundredfold the van der Waals forces. 

If vibrations which are also simultaneously active in the infrared spec- 
trum participate in the combinational scattering (Raman scattering), it is 
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inadmissible to ignore the forces of the resonance type. In molecules 
not having a center of symmetry, many frequencies are simultaneously 
active in both spectra. The forces of the resonance type (quadrupole- 
quadrupole) will be comparable with the van der Waals forces only in 
the study of the influence of intermolecular forces on the combinational 
scattering (Raman) spectrum of molecules having a center of symmetry. 
However, in this case, the interaction effect will be very small, as was 
shown in the preceding paragraph. In particular, the magnitude of the 
line splittings must be small in combinational scattering (Raman scatter- 
ing) of crystals containing few molecules per unit cell. 

Nedungadi“‘ carried out a thorough investigation of the intramolecular 
vibrations of the naphthalene monocrystal. When studying the polari- 
zation of the lines of combinational scattering (Raman scattering) it 
became evident that in the crystal there occurs a splitting of the lines 
512 and 1576 cm-!. Here the difference of frequencies between the split 
components of both lines equals 4 to 5 cm-!. A similar splitting was also 
observed for the lines 1384 and 3058 cm7!. Nedungadi considers that 
the splitting effect is a general property of internal vibrations; however, 
the splitting of components is in many cases so small that it cannot be 
observed. These deductions find a theoretical confirmation in our work. 

In crystals which consist of molecules not having a center of symmetry, 
one must expect great changes for the vibrations simultaneously active 
(in the isolated molecule) in the infrared and combinational scattering 
(Raman) spectra. Owing to the quasi-crystalline structure of liquids, 
the splitting of frequencies must be smeared out, and as a result, the lines 
will be wide. Since molecules with a dipole moment do not have a 
center of symmetry, the displacement and splitting of frequencies must 
be greater for such molecules than for molecules possessing a center of 
symmetry. The presence of a permanent moment is not essential for a 
widening of the lines connected with the splitting of vibrations of free 
molecules. The only indispensable requirement is that an induced 
moment occur with the vibration. We must, however, keep in mind that 
another essential cause for the widening of lines also exists and shall be 
described in the following section. 


11. Several Generalizations 


In the three preceding sections, we suggested a few simplifications in 
order to facilitate calculation. Let us consider how the results obtained 
there will vary if we discard these simplifications and which of the 
simplifications are of little importance. 

1. We proposed that the frequency of normal molecular vibrations wo 
under study satisfies the inequality [(2.6)]. This allowed us to consider 
only the interaction of molecular vibrations of frequency wo. 


30 DAVYDOV: THEORY OF MOLECULAR EXCITONS 


Let us now assume that there are two normal frequencies, wo; and wo», 
located sufficiently close together so that the inequality [(2.6)] is not ful- 
filled. It is now necessary to consider the interaction of each of the 
normal frequencies (wo1 and wos) of the molecule, together with vibra- 
tions of types wo1 and wo: of other molecules. In place of the system of 
oN equations, we shall obtain the system of 20N connected equations 
(here, for simplification, we assume V = Q): 


Gat white + Df De, mall) ang + Diemp(12)anal = 0, 
ree ae 
G+ Wh2Gea + a D2 mp(22) Ging + D7, ma(21)ane] = 0. 
mes ae 
Solving this system of equations by substituting 


a}, = asetmw,, ak, = afela=un, 


we obtain the system of 2c equations 
(wor — w?)ag + » [T2p(11)a5 + T2h(12)as] = 0, 
(wor — w)a2 + or [124,(22)az + 124(21)a}] = 0 
This system will have nontrivial solutions if the determinant, composed 


of the coefficients of the unknowns, is equal to zero. We shall obtain 
an algebraic equation of the 2oth degree with respect to w: 


A(w*) = 0. 
The solution of the latter yields 2c bands of frequencies: 
w= wk), pp =1,2,..., 20. 


Thus, in this case also, to each of the normal molecular frequencies 
will formally correspond its own o bands of frequencies. However, the 
width and distance between them will depend not only on the vibrations 
connected with one of the normal frequencies, but also on the vibrations 
corresponding to another normal frequency. In particular, one can also 
expect noticeable displacements and splittings of the combinational 
scattering (Raman) frequencies in molecules having a center of sym- 
metry, if the molecule possesses a frequency in the infrared spectrum 
which happens to coincide with or is close to the Raman frequency (see 
Section 35). 

2. A second, more essential point was our assumption that the mole- 
cules were rigidly fixed in the lattice points (i.e., that the mass of the 
molecule was infinitely large). This allowed us to exclude from con- 
sideration the connection between the intramolecular vibrations and the 
external vibrations (rotations and translations) of the molecules. This 
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assumption is sometimes very crude and results in the loss of several 
physically interesting possibilities. These we shall briefly (and incom- 
pletely) indicate here, leaving a more detailed examination for later sec- 
tions, in which we shall develop the general quantum-mechanical theory 
of light absorption in the infrared, visible, and ultraviolet regions. 

Upon excitation of the intramolecular vibrations, the forces of inter- 
action with the surrounding molecules change. The change in interaction 
forces results in a local lattice deformation in the region of the vibrating 
molecule and, consequently, in an excitation of external vibrations (lattice 
vibrations). Thus, when taking into account the possibility of the dis- 
placement of molecules from their equilibrium positions, it is impossible to 
consider intramolecular vibrations as independent from lattice vibrations. 

The transfer of energy from intramolecular vibrations to lattice vibra- 
tions is equivalent to the damping (vibrations with friction) of intra- 
molecular vibrations. Asa result, in the resonance points w, = w,(0) of 
the expression (2.22), the amplitudes will now be finite, and the resonance 
curve will have a finite width. Consequently, the absorption bands will 
be wide. The widening of absorption bands will be simultaneously 
accompanied by a decrease in the width of the bands and a decrease in 
the distances between the bands of frequencies caused by the splitting 
of normal frequencies. Therefore, the above-mentioned polarization 
effects will be less noticeable. 

The width of the absorption bands will be greater as energy is more 
rapidly transferred to the lattice vibrations. This last process depends 
on the correlation of two magnitudes: (a) the time of transfer 7; of intra- 
molecular vibrations from molecule to molecule, which is inversely pro- 
portional to the group velocity of running vibrational waves [(2.23)]; 
(b) the time of displacement rz of molecules into new equilibrium posi- 
tions, which depends on the masses of the molecules and on the change 
of the interaction forces when intramolecular vibrations occur. 

If 7: > ra, a local lattice deformation occurs. This results in a transfer 
of the energy of intramolecular vibrations to the lattice vibrations, and 
vice versa, at high temperatures. 

If rs < ra, the intramolecular vibrational state is transmitted from 
molecule to molecule at such speed that the displacements of the mole- 
cules from their equilibrium positions cannot attain significant magnitude. 
The probability of transfer of energy of intramolecular vibrations to 
lattice vibrations will be small. 

3. The results obtained in this chapter may also be extended to spectra 
in the visible and ultraviolet regions, since it is well known that a rough 
picture of the interaction of molecules (atoms) with light may be obtained 
by having a certain virtual oscillator correspond to each quantized tran- 
sition in the molecule. The interaction of molecules upon the absorption 
of light will then be reduced to the interaction of a system of oscillators. 


CHAPTER III 


QUANTUM THEORY OF ABSORPTION AND 
COMBINATIONAL SCATTERING (RAMAN) 
SPECTRA IN MOLECULAR CRYSTALS 
(ZEROTH APPROXIMATION) { 


12. Energy Terms of a Molecular Crystal with Fixed Molecules 


We shall begin the investigation of absorption spectra of molecular 
crystals by means of quantum-mechanical methods, assuming that the 
molecules are rigidly fixed in equilibrium positions (temperature of the 
crystal near zero, and mass of molecules very great). In Chapter IV the 
limits of applicability of such an approximation will be shown. 

Let us examine, as we did in Chapter IT, a molecular crystal contain- 
ing o identical molecules in a unit cell. We further consider the station- 
ary energy states of the isolated molecule, i.e., the wave functions g,, and 
the energy levels Ef,,, to be known. The symbol f determines the char- 
acter and number of the excited state. This may be one of the vibra- 
tional excited states, a purely electronic excitation, or an electronic- 
vibrational excitation. 

In order not to complicate calculations, we at first assume that the 
molecules do not have degenerate states. For many molecules of a low 
symmetry class (which will chiefly interest us), this is actually the case. 
It is not very difficult to extend the generalization of the theory to the 
case of degenerate excited states. If there are degenerate levels in a 
molecule, the degeneracy may be removed when the molecule is situated 
in a crystal of suitable symmetry. This results in a splitting of molecu- 
lar terms, similar to the splitting of atomic terms already studied by 
H. Bethe.4 Splitting of this type is analogous to the splitting of atomic 
terms in an electric field. It is essential in ionic crystals but cannot play 
a significant part in molecular crystals consisting of neutral, nondipolar 
molecules. Study of this type of splitting does not come within the scope 
of our work. Therefore, we shall limit ourselves to examining the non- 
degenerate energy states in an isolated molecule, and we shall show that 

+ The results of the first two sections of this chapter have been published by the 


author in a somewhat less general form.*® 
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even in this case, as a result of the presence in the unit cell of several 
differently oriented molecules, a characteristic splitting of the bands of 
light absorption by crystals appears. 

The wave functions o£, of the molecule will depend on the relative 
coordinates of the nuclei of the molecule and on the coordinates and spins 
of electrons that make up the molecule. For the sake of brevity, all these 
coordinates will henceforth be called internal coordinates of the molecule 
and will be designated by the letter ry. As is known for molecules, one 
can always separate the group of electrons whose energy states determine 
the light absorption in the long-wave spectral region and whose wave 
functions overlap weakly with the wave functions of the remaining elec- 
trons. We shall call these electrons optical. In many organic com- 
pounds these electrons are the so-called z electrons. Let us assume that 
each molecule possesses S optical electrons. 

We shall designate by H,. the operator for the energy of the molecule 
occupying the ath place in the nth cell of the crystal and by Viema the 
operator for the energy of interaction between the molecules na and m8. 
The wave functions and the energy levels — of the stationary states 
of the entire crystal are then determined from the quantum-mechanical 
equation: 

(S Ane +34 SY’ Vnams — E) & = 0. (3.1) 

na na,mBp 
Here the summation is carried out over all oN molecules that make up 
the crystal. The prime of the second sum shows that the term for which 
simultaneously m = n and a = Sisabsent in the sum. The wave func- 
tion @ depends on the internal coordinates r of all the molecules of the 
crystal and of the coordinates determining the spatial locations of the 
molecules; we consider the latter coordinates fixed. 

If the interaction between molecules is small, it may be neglected in 
the zeroth approximation. Then the wave function of the normal state 
of the crystal may be presented in the form of antisymmetrized products 
of the wave functions of g2, normal states of isolated molecules. 


©y = [(SeN)I]-* ) (—-1)"P. Mo, (3.2) 
where VY = [| 2. (3.3) 


?,, is the nonantisymmetrized wave function for a separate molecule, 
and P, is one of (ScN)! possible permutations of electrons. The permu- 
tations are numbered arbitrarily but in such a manner that each suc- 
ceeding one is obtained from the preceding by permutation of one pair of 
electrons; the summation is carried out over all possible permutations 
of the electrons. 


34 DAVYDOV: THEORY OF MOLECULAR EXCITONS 


With the aid of (3.2) and (3.1), the usual methods [substitution of 
(3.2) in (3.1), multiplication by the complex conjugate of function (38.3), 
and integration over all internal coordinates of the molecules] enable us 
to obtain the energy of the normal state of the crystal in the first 
approximation: 


Bo = { Be, + 6 Y [f lea) |?Vna.mal@ma (11) |? dr 
mp 


— Yi f BED) a2. (IL) Vine ma @eh(I) G2 9(1) dr ]}» (3.40) 


where @no(I) and gme(II) are the wave functions of the molecules na and 
m@ when the electrons are in a certain arbitrary position in the molecules, 
and $,2(II) and @ng(I) are the same functions upon permutation between 
the molecules na and m8 of any pair of electrons. The sum over yu covers 
all possible permutations between the molecules na and mf of pairs of 
electrons. { 

The expression (3.4a) is derived by assuming that the wave functions 
@°, of isolated molecules do not change when a crystal is formed. In 
reality, of course, these functions do change as a result of the mutual 
polarization of the molecules. As is well known, this change is the cause 
of van der Waals attraction between molecules. 

Therefore, to allow for this effect, somewhat different wave functions 
should be substituted for the wave functions 2, in (8.3): 


Pra — Pra = C° Pew + Une (3.3a) 


To find a correction for the wave function, the variational method of 
quantum mechanics could be utilized. For this purpose, one must resolve 
Una into a series over the full system of orthogonal functions 94, for an 
isolated molecule 


Una = By oD (3.3b) 
£0 


alld use instead of (3.3) the wave functious of (3.3a); then instead of 
(3.4a), we obtain the new expression 


eae by fea on % » [ / lpee(L)|?V nemal Pe (I1)|? dr 
—¥ f 8D eza(I) Vana moons) ena(1) dr}}» (3-4) 
where Bea = (0°)*Bta +S (0) Ele. (3.40) 
F#0 


+ If each molecule has S electrons, the number of terms in this sum will be equal 
to 8%. 


ABSORPTION AND RAMAN SPECTRA IN CRYSTALS 35 


The unknown coefficients ef for molecular crystals must be small magni- 
tudes, and c° + 1. These coefficients are determined from the minimum 
condition for (3.4), relative to the variation of these coefficients. 

However, inasmuch as the wave functions ¢0, for the isolated mole- 
cules are usually unknown, the solution of the problem formulated by us 
does not have practical meaning. It is more convenient to consider that 
the energies [(3.4b)] and the functions [(3.3a)], corrected in this manner, 
enter into the integral [(3.4)] determining the energy of the system’s 
normal state. These energies and functions take into account the mutual 
polarization of the molecules and, consequently, are functions not only of 
the internal coordinates of the molecules, but also of their mutual dis- 
tances and orientations. In other words, instead of (3.4a), we will use 
(8.4), considering that the values #°, and ¢°, entering into (3.4) are 
known. 

If one of the molecules of the crystal (for instance na) passes to the fth 
excited state, and if all molecules are sufficiently distant from each other, 
such an excited state of the crystal will be described by the wave function 


Xba = [(SoN) ]-* ¥ (—1) "pe, (3.5) 
where wi, = of. [] eee. (3.6) 
ae 


When the molecules are drawn together to the distances corresponding 
to the normal location of the molecules in the crystal, the excitation will 
no longer be localized in the molecule na but will be transmitted from one 
molecule to another. The transmission of the excitation from molecule 
to molecule is analogous to the spreading over the crystal of excitation 
waves, or excitons according to the terminology first introduced by 
J. I. Frenkel?°:?1 and now universally recognized. 

The crystal states, differing only in the site of the molecular excitation, 
belong to one and the same energy of the crystal. In other words, the 
wave functions [(8.5)], which differ from each other by their indices na, 
form a degenerate system of functions in the zeroth approximation. 
Therefore, the wave function of the excited state of the crystal (in the 
zeroth approximation) may be written in the form of a superposition of 
states [(3.5)]. 


® = (oN)-# Y dnaxde, (3.7) 


where |G@na|2/oN determines the probability of the molecule na being in 
an excited state for a given location of the molecules in a crystal. In the 


t Here also it is assumed that the functions gf, and ¢°, take into account the 
effect of mutual polarization of the molecules. 
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ideal crystal lattice formed from identical molecules, the probability that 
precisely molecule na will turn out to be excited does not depend on the 
indices na; that is, 

|Gna|? = const. (3.8) 


Substituting (3.7) into (3.1), multiplying by x4,, and integrating over 
internal molecular variables, gives a system of algebraic equations for 
determining @n. and the energy E/ of a system which is in an excited state: 


YY Moa mpna — Omp = 0. (3.9) 
Here the summation covers all molecules except the mth; the matrix 
element 
Mh me = [ PSL) e8a(L) Vnemayet (IL) efq (IL) dr 
= > f ph(1) one (IT) Vnameyne(Il) y%,(I) dr (3.10) 
Lb 


determines the exchange of excitation between the nath and mth 
molecules. 


Bl = + Bhat Y {Eoa t+ f lehe(D)[?V¥ nama|@re(II)|? ar 
— ) J 85D e2s(ID Vana maesi (IL) ofa (L) dr 
tine YL leSe(D IV anna) a 
— > f oP) e2p(ID Vnemayet(II)¢2,(1) dr ]}- (8-11) 


The prime in the sum indicates that the summation is carried out over 
all molecules except the mth. 

By subtracting the energy of the normal state [(3.4)] from (3.11), we 
obtain the excitation energy of the crystal: 


AE! = AES, + Dig + e, (3.12) 
where AES, = EL, — Ex, 
is the excitation energy of one molecule,t and the expression 


Dia = YT f lee)? Vna.ma|e2a (1) |? dr 
— f le8eC)|?Vnemalyea(II)|* dr — exchange terms] (3.13) 


+ Generally speaking, this energy differs from the excitation energy of an isolated 
molecule because, as a result of mutual polarization of molecules in a crystal, their 
energy levels change. However, in molecular crystals, as we have previously shown, 
this change must occur in the second approximation of the theory and will be 
insignificant. 
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determines the difference in the interaction energy between the excited 
and normal mf molecule and all the remaining molecules. 

Thus the excitation energy [(3.12)] of a crystal differs from the exci- 
tation energy of a molecule by the two terms Dyg and e. 

To calculate «, one must solve the system of equations (8.9), which are 
homogeneous relative to daa. By virtue of the translational symmetry 
of the crystal, the system of equations (8.9) can be reduced to a simpler 
system with the aid of the substitution 


= Baek, (3.14) 


where k is the wave vector characterizing the excitation wave in the 
erystal and running through N discrete values. 


Here »; are whole numbers satisfying the inequality 


NE -N; oes 
=e SWS og? 2 = 1, 2, 3, 
and az} are vectors of the reciprocal lattice. 


The coefficients B. [(3.14)]} satisfy the system of equations 


> Lig(k)Ba — &B, = 0, (3.15) 


where Lite) = y Mine (8.16) 


The prime in the sum [(8.16)] indicates that if a # 8 the summation is 
carried out over all the values of , including n = m, but that if a = 6 
then the summation is over all n ¥ m. 

The conditions under which the system of equations (3.15) has non- 
trivial solutions leads to the equation of the oth degree relative to e: 


A(€) = |ELg(k) — e845] = 0. (3.17) 


After solving the (3.17), we obtain o values of e, as functions of the wave 
number k: 
ef = ef(kr), Bm=1,2,...,¢. (3.18) 


The energy distances between the values e/, corresponding to various u, 
are determined by the matrix elements Li f(Ic), which, with the aid of 
(3.16), are expressed through the matrix elements for the transfer of exci- 
tation between two molecules. In molecular crystals, these matrix ele- 
ments are small, even for neighboring molecules (~-0.01 ev and less). 
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For a rough calculation of the matrix elements of (8.10), we can utilize 
the decomposition in a power series of inverse distances between mole- 
cules Rzing of the classical interaction energy Vnems. If we limit our- 


selves merely to the first term of the series (i.e., dipole-dipole interaction), 
then 


e? ee aig Be ms fe 
Vnamé = 3 > (224 2p a Craving as Yra¥ng)s (3.19) 
na,mB ; 
where Ler Ynay ena and Ling, Ying eB 


are the coordinates of electrons 7 and j in two molecules, relative to the 
systems of coordinates located in the centers of the molecules, so that 
the line linking the centers of the molecules coincides with the 2 axes of 
both systems. 

If we substitute (3.19) in (8.10), we are immediately convinced that, 
for dipole transitions in the molecules, only the first integral, which deter- 
mines the resonance interaction between two molecules, will have practi- 
cal meaning. With this approximation, we obtain the matrix element 
[(3.10)] in the following form: 


tee, te 


— cos 07,4 cos 8%, — cos O%, cos O¥,), (3.20) 


2 
ME ccimB =~ ps 


na,mp 


2 
(2 cos 02,4 cos 0%, 


where efyftre?, dr is the matrix element of the electric moment vector 
determining the quantum transition within the molecule, and cos 0%, and 
cos 62, are the cosines, formed by the electrical moment vectors in the 
molecules m§ and na, respectively, with the axes of the coordinates. 

The square of the matrix element for the transition may be expressed 
through the oscillator strength Fy for the corresponding transition and 
through the cyclic frequency of the transition w with the aid of the known 


formula 
2 
| i. e*re dr 


where yu is the mass of the electron. Then the expression (3.20) assumes 
the form 


_ Fogh 
~~ Que 


d 


e* DF oy 


f eee re LL 
BL 2noks, nb 


na,mgp — 


(2 cos 07,4 cos 02, 
— cos 07, cos 07, ~ cos 04%, cos O4,). (8.20) 


For dipole forbidden transitions, fef*ren.dr = 0. In this case, an 
interaction of a higher order (quadrupole-quadrupole, etc.) must be taken 
into account in the decomposition of (3.19). It is then not possible to 
discard the second term in (3.10), since in certain cases it will have a 
value identical with the first term in order of magnitude. 
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Thus, to the molecular excitation energy E4, in the crystal, » quasi- 
continuous bands of excited states will correspond: 


AES(k) = AEL, + Dig + e{(k). (3.21) 


By substituting the » values of e from (3.18) in the equations (3.15), 
we obtain uw systems of coefficients B+, which determine u wave functions 
for excited states of the crystal, corresponding to the fth excitation of the 
molecule 

$/(k) = («N)-* > Bt xl ek, (3.22) 


To determine the selection and polarization rules for quantum tran- 
sition from a normal state ®° to an excited state ®{(k) of the crystal under 
the action of light, we must calculate the matrix element 


(uf|rlo) = Jfel*(k)r&° dr. (3.23) 


Here we have disregarded the wave vector Q of light for the same reasons 
as in Section 8.f 

After substituting the values for the wave functions into (8.23), we 
obtain 


(uflrlo) = (oN) et YY Be f oftres, dr. (3.24) 


The integrals entering into (3.24) do not, depend on n; therefore, the 
matrix element [(3.24)] will differ from zero only when 


k = 0. (8.25) 


Thus the transitions may occur only at the points of quasi-continuous 
zones of energy ¢,(k), which correspond to the valuesk = 0. This agrees 
with the classical statement: Only the bounding (limiting) frequencies are 
excited by light. 

Thus the matrix element [(3.24)] must have the form 


N-H(uflrlo) = o-# Y Be f pftrg®, dr. (3.26) 


For a further investigation of the matrix element [(3.26)], we need to 
know the values of the system of coefficients B# for each u. 


13. Application of Group Theory for Determining the Selection and 
Polarization Rules of Split Components upon Absorption and 
Combinational Scattering (Raman Scattering) of Light in Crystals 
Which Contain Two Molecules per Unit Cell 


The solution of equation (3.17) to determine ¢/(k) and the coefficients 
Bt is comparatively simple in the case of a small number of molecules in 


tA more rigorous investigation, which considers Q ~ 0, will be carried out in 
Sections 21 to 23. 
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a unit cell. To establish a qualitative picture, it is sufficient to obtain 
the values of the coefficients Be in (3.22) and to determine the symmetry 
properties of the wave functions [(3.22)] by utilizing group theoretical 
methods. After determining the symmetry of the wave functions for the 
energy bands distinguished by the designation », we can establish the 
selection rules and polarization of the corresponding transitions under 
the action of light. 

In lieu of the immediate solution of equations (8.15) and (3.17), we 
may use the following method to find the coefficients Bé. 

Let us write the wave functions [(3.22)] for the case of two molecules 
per unit cell and k = 0: 


B= ON) Y OBE + x6,BD) 


Ne 2 


= (2N)- y (xf, B2 + xf,B3) — 


The functions [(3.27)] must be orthogonal to each other and normalized; 
i.e., they must satisfy the conditions 


JOLOL dr = by. (3.28) 
Furthermore, by virtue of (8.28), 
[Be] = | By). (3.29) 


Utilizing (3.27), we shall, with the aid of (3.28), obtain the following 
equations: 
|Bt|? + [Bil? = 2 
|B3? + [B32 = 2 (3.30) 
BiB? + BiBz = 0 


Equations (3.30) and (3.29) satisfy the system of coefficients 


Me BY By 
1 1 1 (3.31) 
2 1 —-1 
Thus the wave functions [(3.27)] will have the form 
é{ = (2N)-* » Og y tay 
(3.32) 


- a Y (xa) 


Inasmuch as the wave functions ¢°, for the normal state of the mole- 
cules belong to the totally symmetrical representation, the symmetry 
properties of the wave functions x4, entering into (3.32) will be deter- 
mined by the symmetry properties for the wave functions ¢/, of the fth 
excited state of the molecule. To determine the irreducible representa- 
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tions according to which the functions [(3.32)] transform, we must know 
the irreducible representations of the wave functions for the excited states 
of the molecules entering into the crystal. 

a. Crystals of the Naphthalene Type. To illustrate the use of group 
theory, let us examine a monoclinic crystal belonging to the space group 
C8, and containing two molecules in a unit cell. The space symmetry 
group C,—the fifth group of the monoclinic-prismatic class—is very often 
encountered among organic crystals. This group has four symmetry ele- 
ments: screw axes of the second order Cy, parallel to the b axis of the 
crystal; centers of symmetry 7 situated in the center, the corners, and 
the centers of the sides of the faces ab and of the plane parallel to the 


Fic. 3. Coordinate axes for the naphthalene molecule. 


face ab and c/2 distant from it; glide planes o perpendicular to the b axis; 
and identity element &. The characters of the irreducible representa- 
tions of this group are given in Table 1 (see Section 6, Chapter I). 

The minimum number of two molecules per unit cell in group C8, is 
realized when the molecules possess a center of symmetry. Then, when 
the centers of symmetry of a crystal and a molecule coincide, the sym- 
metry operations of the crystal will transpose one molecule to the place 
of another or will leave them in place. 

Crystals of this type include naphthalene CioHs, anthracene Ci.H io, 
biphenyl (CsHs)2, ethylene C2H,, paradichlorobenzene CsH,Cl2, and 
others. 

By way of example, let us examine a crystal of naphthalene. All 
qualitative results obtained in this case will apply equally to other 
crystals of this type. 

The naphthalene molecule has a planar structure which is schematically 
depicted in Figure 3. The carbon atoms are indicated by white circles 
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and the hydrogen atoms by black ones. The molecule belongs to the 
point symmetry group D»,. The characters of the irreducible representa- 
tions of this group, together with the designation of the transformation 
properties for the components of the radius vector (2, y, 2), are given in 
Table 2. 


TaBLE 2. CHARACTER TaBLE FOR THE D2, Point Group, IncLupING NUMBER OF 
NorMalL VipraTIONAL Mopgs aNp THEIR SELECTION RULES FOR NAPHTHALENE 


Combina- 
Dy» |E| cz | cy} cz] i | ot | ov | ot | N| Infrared nie 
(Raman) 
Aig 1 1 1 1 I 1 1 1 | 9 | Forbidden | Permitted 
Bi, 1] -1] -1 1 1; -1/ -l 1 | 4 | Forbidden | Permitted 
Aw 1 1/ -1] -1 1 1 | —1 | ~—1/ 8 | Forbidden | Permitted 
Bog 1} ~1 1] -1 1] -1 1} —1 | 3 | Forbidden | Permitted 
Aiu 1 1 1 1/ —-1} -1) —1 | —1 | 4 | Forbidden | Forbidden 
z, Bru 1; ~1/] -1 1} -1 1 1 | —1 | 7 | Permitted | Forbidden 
z,Au | 1 1/ —-1/ ~-1] -1) -1 1 1 | 4 | Permitted | Forbidden 
y, Bu | 1] —1 1 ee 1; -1 1} -1l 1 | 8 | Permitted | Forbidden 


The number of intramolecular vibrations N belonging to each irreduci- 
ble representation are also noted in this table. The last two columns 
show whether the corresponding vibrations are forbidden or permitted in 
the infrared and the combinational scattering (Raman scattering) spectra. 

As the author has shown,‘** the excited states for z electrons in the 
naphthalene molecule may belong to four irreducible representations, 
Ay, Ag, Biz, and Bo, of the symmetry group Dy,. The normal states 
of these molecules belong to the totally symmetrical irreducible repre- 
sentation Ay. Optical transitions (of the dipole type) from a normal 
state are permitted only to the excited states By, and B., (with polari- 
zation along the long and intermediate molecular axes, respectively). 
The transitions to the states Ay, and Ae, are forbidden. 

These results are obtained by the methods of group theory and do not 
depend on the approximation method of calculating the molecular energy 
terms themselves. 

Knowing the symmetry of the wave functions of the electronic tran- 
sitions in the molecule, one can determine, with the aid of group theory, 
the symmetry properties of the corresponding wave functions of the 
erystal and, consequently, can determine the selection rules and the 
polarization of the spectral lines of the crystal. 

The wave function of the normal state of the crystal will transform as 
the totally symmetrical representation A, of group C§, (Table 1): 


bw A,. (3.33) 
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Naturally, the character of the excited states of the crystal will depend 
on the type of the corresponding excited molecular state. Therefore, 
for the upper index f of the wave functions of the crystal ®/ we substitute 
the symbol for that irreducible representation of the molecular symmetry 
group to which the corresponding molecular transition pertains. 

Let us take into consideration the fact that the identity operation # 
and the inversion 7 coincide in the molecule and the crystal. The trans- 
position of the molecules with indices 1 and 2 and the subsequent turn 
by 180° of the molecule around its intermediate axis correspond to the 
crystal symmetry operations C, and o (i.e., the operation C¥ of the mole- 
cule, and correspondingly, the reflection in the plane, perpendicular to 
the same axis is the operation o” of the molecule). Then, utilizing the 
table of characters of the irreducible representations for the symmetry 
group of the molecule, we obtain the following transformation properties 
of the wave function of the crystal 6?™: 


Eps = Bu Cob P™ = — Pew 
IDB = — Bu ob = paw, 


If we compare the results obtained with the character table, we are 
satisfied that the wave function of the crystal 6?" transforms like the 
irreducible representation B, of the crystal symmetry group. This may 
be written in abbreviated form: 


Ef ~ By, (3.34) 


We similarly establish the transformation properties of the remaining 
functions, which determine the excited electronic states of the crystal: 


bi", df" ~ A, 
4", BA" ~ B, 
bP, be ~ Au 
BPs, SBM WB, 


(8.35) 


Using group theory, we may at once show the selection rules and 
polarization of the corresponding optical dipole transitions in the crystal. 
All transitions from a normal state to excited states of a crystal are for- 
bidden with wave functions &¢", 64%, 64%, and é4%. The transitions to 
all other levels are permitted. Thus the transitions in the crystal are 
forbidden if they correspond to forbidden transitions in the molecule, 
and they are permitted if they correspond to permitted transitions in the 
molecule, for the case of two molecules per unit cell. 

Two permitted transitions in the crystal correspond to each per- 
mitted transition in the molecule. In other words, the energy term of 
the molecule in a crystal is split into two terms, which differ from each 
other in polarization. One term, pertaining to the irreducible represen- 
tation A, of the crystal, is polarized in such a fashion that the vibrations 
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of the electric vector have a direction parallel to the 6 axis of the crystal. 
The other term is in the plane perpendicular to the b axis. The splitting 
of molecular terms is caused by the nonidentical spatial orientation of 
anisotropic molecules in a crystal and by the resonance interaction 
between molecules. 

Similarly, results pertaining to the infrared spectra may be obtained. 
The absorption is permitted for molecular vibrations pertaining to the 
irreducible representations A2,, Biz, and B», of the molecule. To each 
type of molecular vibration in the crystal, two vibrations with a different 
polarization will correspond. Table 3 shows the irreducible represen- 
tations to which the split intramolecular crystal vibrations belong and 
their polarizations. 


TasBLe 3. CorRRELATION TABLE FOR SYMMETRY SPECIES OF INFRARED ACTIVE 
NAPHTHALENE VIBRATIONS IN THE Mo.LecuLe (Do,) AND THE CrystaL (C%,) 


re 
Molecular | Wave function for| Irreducible | Direction of vibration 
vibration | crystal vibration | representation in crystal 
a of Bu 1b 
ou Ca Ay i b 
B peu By 1b 
i ay Au I| 6 
B pu Ay || b 
ae ee Bu Lob 


It must be noted that, in a crystal, molecular vibrations in the isolated 
molecule belonging to an irreducible representation A, will also occur in 
the infrared absorption spectrum. In the naphthalene molecule there 
are four such vibrations, and they are forbidden not only in the infrared 
spectrum, but also in the combinational scattering (Raman) spectrum. 
To each of these vibrations, in a crystal there will correspond two vibra- 
tions with wave functions 


iu ~w A, vibration parallel to b axis, 
4% ~ B, vibration perpendicular to b axis. 


However, inasmuch as the separations between split components are pro- 
portional in the first approximation to the squares of the dipole moments 
for the corresponding molecular vibrations (in the given case they are 
equal to zero), these separations must be very small. Therefore, in the 
absorption spectrum, the split components corresponding to the molecu- 
lar vibrations of the Ai, type apparently cannot be permitted. 

The vibrations occurring in combinational scattering (Raman scatter- 
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ing) of light may also be analyzed in the same manner. Table 4 shows 
the irreducible representations to which the split components of vibra- 
tions in an isolated molecule belong. 


TaBLe 4. CoRRELATION TABLE FOR SyMMETRY SPECIES OF COMBINATIONAL 
ScATTERING (RAMAN AcTIVE) NAPHTHALENE VIBRATIONS IN THE 
Mo ecuLe (Dx) AND THE CrysTAL (C3,) 


Molecular vibration | Crystal wave function | Irreducible representation 
As of x 
Bu a iy 
An oe i 
Be oe B 


From this table we may infer that each type of vibration occurring in 
the combinational scattering (Raman) spectrum of an isolated molecule 
splits in the crystal into two vibrations. One of these belongs to the 
totally symmetrical type A,. The separation between split components 
of the vibrations must be small, since the corresponding vibrations do not 
have dipole moments. 

Nedungadi’s** experiments confirm the results obtained here. 

b. Crystals of the Calcite Type. Let us now consider crystals of the 
calcite type, which possess a rhombohedral lattice with a spatial sym- 
metry group D§, and contain two molecules per unit cell. Crystals of 
many carbonates of divalent metals (CaCO3, MgCOs, etc.) possess this 
type of lattice. The infrared spectra of such crystals have been 
repeatedly investigated.‘7—* 

Tt has been established that, in the wavelength region 3 to 15 y, infrared 
spectra of such crystals are identical. This was explained by the fact 
that in this wavelength region, the spectrum is determined by the vibra- 
tions of separate atoms within the ionic molecule CO; ~-. 

The binding of particles within the ion CO;~ is significantly firmer 
than the binding between ions in the lattice. Therefore, we may con- 
sider this ion an independent molecule. Furthermore, we may apply the 
methods developed above for molecular crystals in studying the change of 
intramolecular vibrations of this ion when it enters into an ionic crystal. 
The ion CO;y~ represents a plane equilateral triangle, in the center of 
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which lies a C atom (Figure 4); consequently, it may be related to the 
symmetry group D3. The molecule possesses the following symmetry 
elements: 1 threefold rotation axis C3, 1 improper rotation axis S;, 3 two- 
fold axes C, (1 Cy), 3 reflection planes o, (crossing C3), and 1 reflection 
plane o, (L C3). Table 5 lists the characters of the irreducible represen- 
tations for this group and shows the 
y characters for the transformations of 

coordinates 2, y, 2. 
Altogether, six internal vibrations, 
of which two pairs are doubly degen- 
x erate, are possible in the carbonate ion. 
The last three columns of Table 5 
indicate the N numbers of the internal 
vibrations for each irreducible repre- 
sentation, as well as the selection rules 
in the infrared spectrum and in the 
combinational scattering (Raman) 
Fie. 4. Coordinate axes for the car- spectrum. The vibration Ai, which 
bonate ion molecule. is inactive in the infrared spectrum, 
has, according to calculations (Ref. 50, 
p. 178), a frequency of 1069 cm! (9.42 2). The vibration Aj’, which 
is active in the infrared spectrum and which has an electric moment 
perpendicular to the plane of the molecule, has a frequency of 879 
em-! (11.4 «). Two degenerate vibrations EH’, with an electric moment 
lying in the plane of the molecule, have frequencies of 1415 cm™ 

(7.06 ») and 680 cm! (14.7 uw), respectively. 


TaBLE 5. CHaracrer TaBLE ror THE Ds, Point Grovr, INcLupINGc NUMBER 
or Norma VisrationaL Mopes and TuHerr SELectTion Russ 
FOR THE CarBONATE ION 


oT 
Combina- 
tional 
Dsa E on | 2CZ | 283 | 8C2 | 80. | N Infrared scattering 
(Raman) 

Ay 1 1 1 1 1 1 | 1 | Forbidden) Permitted 

A, 1 1 1 1) -1)| -1 0 | Forbidden | Permitted 

Ay’ 1 -1 1) -1 i1|{ -1 0 | Forbidden | Permitted 

z, Ay’ 1 | -1 1} ~1{ -1 1 | 1 { Permitted | Forbidden 
z+ty, EB’ 2 2} -1] -1 0 0} 2 | Permitted | Forbidden 
&” 2 —2)} -1 1 0 0} O } Permitted | Forbidden 


The ions of CO;~— are located in the lattice in such a way that their 
planes are parallel to each other (Figure 5) and perpendicular to the main 
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axis, which coincides with the diagonal of a rhombohedron. Figure 5 
shows only the ions of CO; ~—; the carbon atoms are indicated by white 


circles, the oxygen atoms by black ones. 


The calcite crystal belongs to the space group D§,, whose characters of 


irreducible representation are listed in Table 6. 


TasBLe 6. CHARACTER TABLE ror THE D§, Spack Group 


Ds, E 2s7 | 2C; 3C. | 30, 
Aig 1 1 1 1 1 
Aww 1 =4 1 =I =] 
Ave 1 1 1 a =] 
z, Aw 1 = 1 =I = 1 
E, 2 a 24 0 
x, y, Eu 2 1 =), 2 0 0 


Of the symmetry operations for the crystal, the following operations 
apply to the radicals. In operation S%, operation S; is applied to the 
radical and the radicals change places and turn by 60° about the z axis. 


In operation C's, operation Cz is applied to 
the radicals. In operation 7, operation S; is 
applied to the radical and the radicals change 
places and turn by 60°. In operation C2, 
operation C, is applied to the radicals. 
In operation o,, operations o, are applied to 
the radicals and the radicals change places. 

If we use these rules governing the appli- 
cation of crystal symmetry operations, as 
well as the tables of characters of the irredu- 
cible representations of the molecule (Table 5) 
and of the crystal (Table 6), we shall find 
the irreducible representations governing the 
transformation of the wave functions [(3.32)] 
for the excited states of the crystal. The 
results obtained are brought together in 
Table 7. 

Since the intramolecular vibrations of the 
ion CO;~- belong only to the irreducible 
representations A}, A;’, and E’, we shall now 
be interested only in the splitting of vibra- 
tions corresponding to these irreducible 
representations (first three lines of Table 7). 


t y 


x 


Fie. 5. Disposition of the car- 
bonate ions in the calcite 
crystal. 


According to Table 7, vibration A{, forbidden in the ion, remains for- 
bidden in the infrared spectrum as well. Furthermore, the permitted 
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vibration A5’ with vibrational direction along the z axis (perpendicular 
to the plane of CO;~—) is also permitted, with like polarization, in the 
infrared spectrum of the crystal. Finally, the degenerate vibrations LE’, 
permitted in the free ion (with a vibrational direction in the plane of 
the radical), remain degenerate in the crystal as well and possess the 
same polarization. 


TaBLe 7. CoRRELATION TABLE FOR SYMMETRY SPECIES OF CARBONATE JON 
VIBRATIONS IN THE MougcuLE (Dj) AND THE CrrsTaL (D$,) 


: Wave function| Irreducible | Vibration in Symmetry of vibration » 
Radical i : combinational scattering 
of crystal representation infrared 
(Raman) 
Al fv Aig Forbidden Totally symmetrical 
par Au Forbidden Forbidden 
Al per" Aw | 2 Forbidden 
° ef?” Axg Forbidden Forbidden 
FI er’ E, Lz Forbidden 
oF’ E, Forbidden Asymmetrical 
a ese Ary Forbidden Forbidden 
: a Aw || 2 Forbidden 
Al ea” Ai Forbidden Forbidden 
1 efi” Aig Forbidden Totally symmetrical 
RB oF” E, Forbidden Asymmetrical 
ee” EB. iz Forbidden 


Thus, despite the splitting of vibrational terms into two components, 
the character of the infrared spectrum does not change qualitatively in 
the crystal; new absorption bands do not appear and their polarization 
does not change. 

The same deductions may also be made for the combinational scatter- 
ing (Raman) spectrum: the spectral character of the crystal and of the 
isolated radical (calculated) must coincide. 

This result is caused by a particular orientation of the molecules in the 
crystal. If the molecular ions in the crystal had another mutual orien- 
tation, it would be necessary to expect qualitative changes in the spec- 
trum. This takes place when the same ion CO;~- enters into another 
type of lattice (aragonite). This case will be analyzed in Section 140. 


14. Crystals Containing Four Molecules per Unit Cell 


Normalization and orthogonalization [(3.29)] conditions for the wave 
functions for excited states, as well as equation (3.8), lead to the system 
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Taste 8. Corrricients Bk ror THE Case OF Four Mo.ecuLes per Unit CELL 


» | Be | Be | Be | Be 


1 1 1 1 1 
2 1 —l 1 —1 
3 1 -1 —1 1 
4 1 1 —1 -l 


of values for the B# coefficients given in Table 8. The wave functions 
for the excited states of a crystal will be: 


{= (4N)-A(xf + xf, + xf, + xf, 
f= (N)-4*xd — xh, + x4, — xf, 
Bi = (4N)-A(xh — xh, — xf, + xh) 
Of = (AN)-¥(xt, + xd, — x4, — xh.) 


| 


(3.36) 


Among crystals containing four molecules per unit cell, we first of all 
consider the crystals of the space symmetry group C3,. The molecules in 
the unit cell of crystals of this type may be arranged in two possible 
ways.5! 

1. If the molecules do not have a center of symmetry, all four mole- 
cules can change places during crystal symmetry operations. Crystals 
of this type comprise phenanthrene, 6-naphthol, and others. 

2. If the molecules do have a center of symmetry, they may occupy two 
pairs of independent centers of symmetry in the crystal, and the crystal 
symmetry operations may transpose among themselves only those mole- 
cules which belong to the same pair. Crystals of this type comprise 
stilbene, azobenzene, and others. 

We shall consider examples of crystals with molecules not possessing 
a center of symmetry and of crystals formed by molecules having a 
center of symmetry. 

a. Crystals of the Phenanthrene Type. The molecule of phenanthrene 
Ci.Hio (Figure 6) belongs to the point symmetry group Cy. It possesses 
the following symmetry elements: a twofold rotation axis C4, a reflection 
plane o* coinciding with the molecular plane, and a symmetry plane o? 
perpendicular to the z axis. The characters of the irreducible represen- 
tations for this group, together with the transformation properties of the 
components of the radius vector (x, y, z), are shown in Table 9. This 
table also gives the number of intramolecular vibrations, N, belonging to 
each irreducible representation. 

Operation E coincides in the crystal and in the molecule. Operation 
C. for the erystal corresponds to a transposition of molecule 1 and mole- 
cule 3 and of molecule 2 and molecule 4 and to a subsequent application 
of operation C¥ to the molecule. Operation 7 of the crystal corresponds 
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TaBLe 9. CHARACTER TaBLE FOR THE Co, Point Group, INcLUDING NUMBER 
oF Norma VisraTIONAL Mopzgs and THEIR SELECTION RULES 
FOR PHENANTHRENE 


Combinational 
Coy E Cy oF o? N Infrared scattering 
(Raman) 
y Ay 1 1 1 1 23 =| Permitted Permitted 
Az 1 1 -1 —1 Forbidden Permitted 
z B, 1 -1 1 —1 22 | Permitted Permitted 
B, 1 -1 -1 1 10 | Permitted Permitted 


to the transposition of molecules 1 and 4 and of 2 and 3 and to the subse- 
quent application of operation o* to the molecule. Operation o of the 
crystal corresponds to a transposition of molecules 1 and 2 and of 3 and 4, 
as well as to a subsequent application of the o7 operation of the molecule. 


Fic. 6. Coordinate axes for the phenanthrene molecule. 


If we use tables of characters of the irreducible representations for the 
symmetry groups of the molecule (C2,) (Table 9) and of the crystal (C3,) 
(Table 1), we can establish for the crystal symmetry group the irreduci- 
ble representations, according to which the wave functions for the intra- 
molecular vibrations transform. 

Table 10 shows the irreducible representations for wave functions of 
vibrational states of the phenanthrene crystal, the direction of the electric 
moment of the vibrations if the corresponding vibration is permitted in 
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the infrared spectrum, the symmetry of vibrations permitted in the 
combinational scattering (Raman) spectrum, and those irreducible repre- 
sentations for the molecular symmetry group to which the given excited 
states of the crystal belong. 


TaBLE 10. ConRELATION TABLE FOR SymMMETRY SPECIES OF PHENANTHRENE 


VIBRATIONS IN THE Mo.LEcuLE (C2,) AND THE CrystaL (C%,) 
Wave function} Irreducible Combinational 
Molecule ‘ Infrared scattering 
of crystal representation 
(Raman) 
et) Ay Forbidden | Totally symmetrical 
ri ett Ay | 5 Forbidden 
ca B, Forbidden | Asymmetrical 
eft Bu Lob Forbidden 
ed? Au |b | Forbidden 
i od? A, Forbidden | Totally symmetrical 
; pat By 1b | Forbidden 
Ca B, Forbidden | Asymmetrical 
@P B, tb | Forbidden 
B eP B, Forbidden | Asymmetrical 
: a5 Au || b Forbidden 
en A, Forbidden | Totally symmetrical 
a B, Forbidden | Asymmetrical 
B pe? Bu tb Forbidden 
: oF? Ay Forbidden | Totally symmetrical 
Cg Au || b Forbidden 


The wave functions for the excited states of the + electrons in the 
phenanthrene molecule belong to the irreducible representations A, and 
B,, Transitions from the normal level to both types of excited levels 
are permitted (there are no forbidden transitions) in the molecule. In 
the crystal, each of the molecular terms is split into four levels, and 
transitions under the action of light will be permitted in only two of the 
four split bands and will be induced by a light wave with an electric- 
vector direction along the b axis and perpendicular to the b axis. 

b. Crystals of the Aragonite Type. The crystals of this type possess a 
rhombic lattice with a V}®° space group (four molecules per unit cell). 
This structure is typical of many crystalline salts of inorganic acids. 
The infrared spectra of these crystals have been investigated repeatedly 
(unfortunately at a temperature not lower than 0°C). Table 11 gives 
the characters of irreducible representations for the V}® space group.: 
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TaBLE 11. CHaracrer TasLe ror THE Vi* Space Group 


Vis E CE cy C3 t Cd ov a* 
Aa 1 1 1 1 1 1 1 1 
Bi 1 ee ee 1 ge |||) eg 1 
ae 1 {|| ete he ee 1 a th ake vib» See 
Boy 1 —l 1 -l1 1 -l 1 -1 
Aw 1 1 1 1 -—1 —-1 —1 -1 
z,Bu| 1 =e ee ||) cet l pin a 
Ang ls 4 | (a ne ee 1 1 
y,Bu| 1 at dbs Czeatee Ill weed el Seed 1 


Figure 7 shows schematically the arrangement of molecules in the 
aragonite crystal.°2 When applying the symmetry operations of the 
crystal, the following operations apply to the ions of COz;-. In oper- 
ation C3, operation C2 is applied to the ion, and then the transposition of 
ion 1 and ion 4 and of ion 2 and ion 3 takes place. In operation C4, 


Fie. 7. Disposition of the metal ions and carbonate ions in the aragonite erystal. 


operation C, is applied to the carbonate ions, and transposition of ion 1 
and ion 3 and of ion 2 and ion 4 takes place. In operation Cz, operation 
C, is applied to the carbonate ions, and transposition of ion 1 and ion 2 
and of ion 3 andion4 takes place. In operation, the C3 crystal operation 
is applied, with subsequent application of operations o, to the carbonate 
ions. In operations o7 and o%, transposition of ion 1 and ion 4 and of 
ion 2 and ion 3 takes place, followed by the transposition of ion 1 and 
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ion 3 and of ion 2 and ion 4, respectively, with the subsequent application 
of operation o, to the ions. 

By using these rules for applying the operations of a crystal, as well as 
the table of characters of the irreducible representations of the molecule 
(D3,) (Table 5) and of the crystal (V}*) (Table 11), we shall find the 
irreducible representations according to which the wave functions [(3.36)] 
for excited states of the crystal transform. The results obtained are 
shown in Table 12. 


TsaBLE 12. CorRELATION TaBLE FoR SYMMETRY SreciES or CARBONATE ION 
ViBRATIONS IN THE MOLECULE (D3,) AND THE CrystTaL (V}5) 


Molecule Wave function Irreducible iepeaceg 
of crystal representation 
ef" Au Forbidden 
, eer Ao x 
A, ear Boy y 
oA Big Forbidden 
par’ Buy Zz 
A eae” Buy Forbidden 
7 gee” Ary Forbidden 
af” Au Forbidden 
ee” Ag + Bu zy 
E! er’ Aig + Big Forbidden 
ee’ Ar + Bi, Forbidden 
ee’ Ao + Bu ZY 


The intramolecular vibrations of the ion CO;~-, which belong to the 
symmetry group Aj’, are split in the crystal into four quasi-continuous 
bands of vibrations. From Table 12, it is apparent that in this connec- 
tion there must appear in the infrared spectrum only one extreme fre- 
quency corresponding to the first band (with a wave function 6#*") and 
that the vibrations will take place parallel to the z axis of the crystal. 
Each pair of degenerate molecular vibrations of the ion CO;~-, which 
belong to the irreducible representation £’, splits in the crystal (by virtue 
of the resonance interaction of vibrations) into four quasi-continuous 
bands of vibrations. Only the extreme vibrations from two of these 
bands will appear in the infrared spectrum. 

In view of the fact that the degeneracy of these vibrations is removed 
while still in the crystal, the extreme vibrations themselves must be 
doubly degenerate. 

Thus, according to the above, four absorption bands (possibly in the 
form of two narrow doublets), which correspond to a degenerate vibra- 
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tion in the molecular ion COy~-, must be observed in the absorption 
spectrum. Two of the bands (one for each component of the doublets) 
must be observed with the electric vector oriented along the x axis of 
the crystal and another two with the electric vector oriented along the 
y axis. The vibrations occurring in the absorption spectrum when the 
electric vector is directed along the 2 axis of the crystal must be single. 
As previously indicated, the investigations of infrared spectra of car- 
bonates were carried out at temperatures not lower than 0°C. In this 
case, the absorption bands are wide and the structural details are lost. 
Table 13 shows the results of the measurements for certain carbonates 
(Ref. 47, p. 276) possessing the structure of the aragonite crystal. In 
accord with theory, Table 13 shows that, for vibrations parallel to the 
2 axis of the crystal, there are no split components for all these crystals. 
For the vibrations lying in the plane perpendicular to the z axis, instead 
of the four absorption bands expected according to theory we observe in 
aragonite three absorption bands in the 1500-cm™! region and two in the 
700-cm-! region. In cerussite there are two bands in the 1400-cm7 
region and in witherite only one absorption band. However, even in 
this last case, the twofold splitting of the bands is clearly observed in 
the majority of overtones; therefore, Schaeffer and Matossi‘? conclude 
that ‘there cannot be any doubt that the frequencies really split in two.” 


Taste 13. Osservep INFRARED FREQUENCIES IN NATURAL CARBONATES 


Frequencies, cm~! (in parentheses, \ in yz) 
Crystal 
Vibrations 1 z Vibrations || z 
1550, 1490 1500 866 
2 (6.46), (6.70) (6.65) (11.55) 
Aragonite 707 712 
(14.17) (14.06) 
Garussite 1420 1370 833 
(7.04) (7.28) 
A : 1460 1460 861 
athens (6.85) (6.85) (11.61) 


eS 


A hypothesis on the deformation of the carbonate ions in the lattice, 
which in our opinion is very artificial and has little basis, was proposed 
to explain this splitting in two of the frequencies. It was assumed that 
the atoms of oxygen form a nonequilateral, isosceles triangle.*® Then a 
certain orientation (corresponding to the symmetry of the crystal) for 
these deformed radicals, relative to the axes of the crystal, was suggested. 
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According to the investigation conducted here, this splitting in two at 
low temperatures must result in the appearance of four absorption bands 
and is explained by the removal of the degeneracy in the field of the 
crystal and by the resonance interaction of the vibrations. 

Furthermore, according to theory, in the 10,000- to 11,000-cm—? region 
of the spectrum of the crystal, one must expect an absorption band 
corresponding to a vibration in the free radical that is inactive in the 
infrared spectrum. 

c. Benzene Crystal. Crystalline benzene, with four molecules per unit 
cell, belongs to the V}5 symmetry group. Figure 8 shows the projection 
of the arrangement of molecules into 
a plane perpendicular to the 0 axis. a 
The plane of the molecular ring is 
approximately perpendicular to the 
b axis and forms an angle of 40° with 
the c axis. The crosshatched mole- 
cules lie in a plane situated at dis- 
tances of 6/2 from the plane in which 
the noncrosshatched molecules lie. 

The wave function for the normal 
state of the crystal, as well as the 
wave function for the normal state of 
the molecule, belong to the totally 
symmetrical representation Ay,. 

Knowing the symmetry of wave 
functions for the excited states of the Fie. 8. Disposition of the projection of 

the four benzene molecules onto the 
molecule, we can, by means of group plane perpendicular to the 6 axis of the 
theory, determine the symmetry crystal. 
properties of the corresponding wave 
functions of the crystal. Consequently, we can determine the selection 
rules and the polarizations of electronic transitions in the crystal. 

The characters of the irreducible representation for the space symme- 
try group V}® of the crystal coincide with the characters of group Ds, 
(Table 2). The characters of the irreducible representations for the point 
symmetry group Dg, of the benzene molecule appear in Table 14. 

Let us take into consideration the fact that identity operation H and 
inversion ¢ coincide in the molecule and in the crystal. To the symme- 
try operations of the crystal C$, C$, o*, and o° there correspond trans- 
positions of molecules 1 and 2 and of 3 and 4; the symmetry operations 
for the molecules C2, C4’, iC2, and 7Cy’, respectively, are subsequently 
applied to the molecule. The application to the molecules of the symme- 
try operation Ci corresponds to the symmetry operation for the crystal 
C3. Now, by using the table of characters of the irreducible represen- 
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TaBLe 14. CuaracTer TABLE FOR THE De, POINT GRouUP 


Dir E| Cy | 2C3 | 2C, | 3C, | 3C;’ | cE | iC | 21Ca | 27, | 3¢C) | 3205 
Ay 1 1 1 1 1 1 1 1 1 1 1 1 
Ai 1 1 1 l 1 a} —-1} —1]} —-1} -1] -1] -1 
An 1 1 1 1) -1} -1 1 1 1 1) -1} -1 
Zz, An | 1 1 1 1} —-1]} —1 | -1] —1 ] -1] -1 1 1 
Bi, 1] =1 1 | -1 1{ -1 1] -1 1| -1 1{-1 
Bu 1|-1 L| -1 1/ —-1] -1 1] -1 1] -1 1 
Ba, 1{| -1 tL} -1] -1 1 1] 1 1} -1] -1 1 
Bou 1] -1 1] -1 ] -1 1| -1 1] -1 1 1] -1 
Ex, 2 2) -1] -1 0 0 2 2{ -1] -1 0 0 
Ex 2 2} -1] -1 0 0] -2| -2 1 1 0 0 
Ey, 2} -2] -1 1 0 0 2) —-2} -1 1 0 0 
(x, y), Hu | 2 | -2] —1 1 0 0} -2 2 1) -1 0 0 


tations for the symmetry group of the benzene molecule, we can deter- 
mine the irreducible representations for the corresponding wave functions 
of the crystal. The results, which pertain to the three possible types of 
wave functions for the excited states of a molecule, are shown in Table 15. 


TasLe 15. ConreELATION TaBLE FoR SYMMETRY SPECIES FoR ELECTRONIC 
States or BENZENE IN THE Mouecure (Dex) AND THE Crystau (V}5) 


depeeecible Irreducible 
representation | Wave function cancoeantaciene Selection rules 
of excited state of crystal P 16 and polarization 
of molecule of group Vi 
3 and @F™ Bou Vibration of electric vector paral- 
Bu lel to b axis 
oe? and BF Aw Forbidden 
@3™ and HF Bu Vibration of electric vector paral- 
Bou lel to ¢ axis 
}3™ and 6B Agu Same, parallel to a axis 
E ef and a Bay + Bru Same, parallel to b and ¢ axes 
ae BF” and &F Ae + Ai Same, parallel to a axis 


Thus, from Table 15, we find that: 

1. The forbidden molecular transition to the level Bi, in the crystal 
becomes permitted and must occur when the electric vector of the light 
wave is directed along the b axis of the crystal. 

2. To the forbidden molecular transition to the level Bo, in the crystal 
there corresponds a doublet; the electric vector of one of the components 
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of the latter is directed along the ¢ axis, and that of the other along the 
b axis, of the crystal. 
3. The permitted doubly degenerate molecular term FE}, in the crystal 
is split and must appear as @ triplet with components a, b, and c. 
According to studies by A. F. Prikhotko, V. L. Broude, and VY. 8. 
Medvedev®* on the absorption spectrum of the benzene monocrystal in 
polarized light, there is a doublet of 37,811 and 37,786 cm! with sharply 
polarized components. This doublet is subsequently repeated in the 
spectrum with a period of 920 cm-}, corresponding to the frequency of 
the totally symmetrical vibrations in the benzene molecule. A compari- 
son of these data with theory enables us to draw the following conclusions: 
1. Monocrystalline laminae of benzene, as used in the experiments, 
possess a surface which coincides with the plane ac of the crystal lattice. 
2. To the electronic spectrum in the near ultraviolet there corresponds 
the molecular term Boy. 


CHAPTER IV 


GENERAL EQUATIONS DETERMINING THE ENERGY 
STATES OF MOLECULAR CRYSTALS 


15. Introduction. Statement of the Problem 


In molecular crystals the molecules preserve their individuality. It is 
therefore convenient to study the energy states of such crystals by start- 
ing from the energy states of isolated molecules. 

When investigating the optical behavior of a crystal in the visible and 
ultraviolet regions of the spectrum, the excited state of a molecule will 
imply both electronic and electronic-vibrational excitation. This excited 
state of the molecule will be described by the wave function 


o, (4.1) 


where the index f indicates the number of the excited state. 

When investigating the optical behavior of a crystal in the infrared 
region of the spectrum, the index f will define the state of intramolecular 
vibration in the crystal. The wave function [(4.1)] must depend on the 
coordinates of the optical electrons and of the nuclei of the molecules, 
and is a solution of the wave equation 


Hol = Erg’, (4.2) 


where H is the operator of the internal energy of the molecule. 

We shall assume that the molecular problem [(4.2)} is resolved, i.e., 
that the eigenfunctions yf and the eigenvalues for the energy E’ are 
known. We shall then investigate the energy states of the entire crystal 
when it is made up of identical molecules. 

At present there are no sufficiently satisfactory methods for solving the 
molecular problem [(4.2)]. All existing methods for calculating the elec- 
tronic states of complicated molecules (valence-bond method, * molecular- 
orbital method, 5*.5¢.46 antisymmetrized-molecular-orbital method,5758 and 
others) have little foundation and yield very crude results. M. V. 
Volkenstein, M. A. Elyashevich, and B. I. Stepanov®? and L. 8. Mayants® 
achieved substantial results in the calculation of molecular vibrational 
states. The methods of investigating molecular vibrations applied by 
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these workers are extremely simple and convenient. These authors suc- 
ceeded in calculating the vibrational spectra of numerous molecules; how- 
ever, it has not yet been possible to extend the same methods to aromatic 
molecules. 

The difficulties that arise when calculating the intramolecular states 
will not be discussed here. We shall study the changes of energy states 
for a large number of molecules disposed at the points of a crystal lattice, 
in comparison with the energy states of infinitely separated (free) mole- 
cules. The energy states of free molecules are here assumed to be known 
(albeit from experimental data). 

When studying light absorption, it is necessary to consider (1) which 
electromagnetic frequencies are absorbed by the crystal and how they 
depend on the direction of the propagation and of the polarization of 
impinging waves and (2) what happens subsequently to the absorbed 
energy—does it again appear in the form of electromagnetic radiation 
(luminescence), or is it absorbed by the crystal and transformed into heat, 
i.e., distributed over all degrees of freedom of the crystal? 

The processes for the transformation of the electronic excitation energy 
of molecular crystals into heat may take place in two ways: (1) intra- 
molecularly and (2) intermolecularly. The intramolecular deactivation 
of the molecules in a crystal, i.e., the transformation of electronic exci- 
tation energy into vibrational energy of a molecule (if such transformation 
takes place), will be assumed known. We shall be interested only in the 
intermolecular deactivation of molecules, i.e., in the transformation of 
excitation energy of molecules into lattice vibration energy. 


16. Quasi-stationary Excited States of a Crystal 


Let us consider a molecular crystal, the unit cell of which is determined 
by the three noncoplanar vectors a1, a2, and a;. Let us assume that in a 
unit cell there are o molecules; the positions of their centers of gravity 
relative to the unit cell corner nearest the coordinate origin are deter- 
mined by the vectors r. (a = 1,2, ...,0¢). The index a characterizes 
the number of the molecule in the cell and the spatial orientation relative 
to the vectors ai, a2, and az. 

The position of a desired molecule in the crystal relative to the coordi- 
nate origin will then be determined by the radius vector 


3 
teot= ) am; + Ya, (4.3) 


where n = m1, M2, N3; the nm: are whole numbers. 
Let the crystal have the form of an oblique-angled parallelepiped with 
axes Nyai, Noa, and N3a3, so that the total number of unit cells in the 
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crystal is equal to N = N,N.N; (here and henceforth the number N will 
be assumed to be very large), and the over-all number of molecules is 
equal to oN. 

We shall consider all molecules of the crystal to be identical and each 
one to contain p atomic nuclei, as well as S optical electrons determining 
the optical behavior of the molecule. 

The state of the crystal will be characterized by the wave function 


Wr, B, t), (4.4) 


which depends on the internal r coordinates of the molecules (coordinates 
of all optical electrons and relative coordinates of the atomic nuclei in 
the molecules), on the coordinates and spatial orientations R of all mole- 
cules in the crystal, and on the time #. This function must satisfy the 
time-dependent equation of quantum mechanics: 


ovr, BR, t) 


th 7 = (Te + Hy)V(r, BR, t) (4.5) 
where Tx is the kinetic energy operator for the rotational and displace- 


ment motion of all molecules. 


ene Be (4.6) 


t=1 


and Ho = Hr, 0/(dr), R] is the energy operator for the crystal of fixed 
molecules. This operator may be written in the following form: 


0 I 
Hy= > Ana (- 2) + 3 > Vaama(r, 2), (4.6a) 


na,mp 


where H,,.[r, 0/(dr)] is the energy operator for the internal state of the 
molecule a, which is located in the nth unit cell of the crystal, and Vaama 
is the interaction operator of the two molecules. 

In solving equation (4.5) we shall proceed from the zeroth approxi- 
mation, which was analyzed in Chapter IIT and in which the molecules 
were assumed to be rigidly fixed (Mf — o) at the lattice points. This 
permits us to indicate the limits of applicability for the zeroth approxi- 
mation and to take into consideration in the first approximation the 
motion (vibrations and rotations) of the molecules in the lattice. 

Thus we shall assume that the solution of the equation for fixed mole- 
cules in the crystal is known: 


HSi(r, R) = EM R)ei(r, R), (4.7) 


where @f(r, R) is the wave function for the electrons of the crystal when 
the molecules are fixed. This wave function forms such a crystal state 
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when »; molecules are in the f;th excited state, when v2 molecules are in 

the f.th excited state, etc., and the other molecules are in a normal state. 
In particular, for crystal states which do not contain excited molecules 

(vy = 0, f = 0), as we have seen in Chapter III, the wave functions are 


&§ = 3(r, R) (4.8) 
and the energy of the crystal 
E(R) = f ®§*H @$ dr = 5) ES, + V(R) (4.82) 


is determined by formulas (3.2) and (3.4). 

The wave function for the state of a crystal containing one molecule 
in the fth excited state depends on the number u of the zone of the 
excited states of the crystal, as well as on the value of the wave number 
of excitons in the zone [(3.22)]. 


bf (k, r, R) = (¢N)-* >. Beyt eck, (4.9) 


and the energy is 
Et(k, R) = | f*(k, 7, R)A,&f(k, 7, RB) dr 


= » ES, + AE, + V(R) + Df a(R) + » zi > MS a mal) ee-™, 


(4.9a) 


If there were (simultaneously) two excited molecules in the crystal, the 
wave function of this state would have the form 


Off (k, k’, 7, R). 


Bp! 


We shall seek the solution of (4.5) in the following form: 
V(r, Rk, t) = 8) >, bo, (t) Ab, exp (—tw ost) 


+ J) diel, )Ag(e, R)B{(h, r, R) exp (—twf¥t) +--+. 4.10) 
xf, ke 

In (4.10) we did not extract the terms corresponding to the excited states 
of a crystal with two, three, etc., excited molecules. These terms play a 

small role if the energy lifting the crystal to an excited state is small. 
The squares of the moduli of the coefficients bj,(t) will characterize the 
probability that at the moment of time ¢ there are no excited molecules in 
the crystal. The state of the crystal is determined by the wave function 


Aos(B)®5(r, R) 
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and the total energy of the crystal is equal to 
Voz = hw oz. 


The function Ao,(#) characterizes the vibrational state of the molecules 
in the lattice and satisfies the equation 


[Tr + EQ(R)]Ao(R) = UorAoe(R). (4.11) 


The squares of the moduli of the coefficients bf*(k, ¢) determine the 
probability that at the moment of time ¢ there is in the crystal an exciton 
(corresponding to the excitation of one molecule) in the uth zone with a 
wave number k, and the lattice vibrations are characterized by the wave 
function Af#(k, R), which satisfies the equation 


[Te + El(k, R))AM(k, R) = UH APC, R), (4,12) 
where UL = haf. 


If we substitute (4.10) in (4.5), multiply by Ao. ?* exp (dwo.t) and by 
ALRwel*(k’, r, R) exp (tuft), and then integrate over r and R, taking 
into consideration (4.11) and (4.12), we obtain the system of equations 


jp, Goel) _ 


a b£( Kt) [O2' [4] exp [—7(wfk# — woe’)t], (4.18) 


rfkp 


where 


[0x’|(a] = i Ave( R) vii Ate(k, R) dR 


+ » i Ava(R) Pfie(s ) Sa dAte(t, R) ap (4.14) 


vit(R) = / @5*(r, R)T 2bL(k, r, R) dr, (4.142) 
ri) = FF i 5O*(r, Ry SPE TB ap es (4.146) 


and 
Pie eb ete, u = Doo(t) [Ox |{2*']* exp [—i(wor — wk’) t] 


Be > bieck, ie" |i] exp [~t(ott — of F')4], (4.15) 


ofl 
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where 
erkisl = f ageqeyter ante aR 


f Ki 
+ » / Ate ce yrgye tm) SABO ae, (4.16) 


yikes = | BL (k', r, R)TRbL(k, 7, R) dr, (4.164) 
Tek ska(g) = -; i(k’, r, R) , R) SSF dr (4.160) 


From equations (4.13) and (4.15), it follows that if, at the moment of 
time t = 0, the system was in one of the excited states A{*(k)@1(k, 7, R) 
with the energy U{*, then it can spontaneously go over into the state 
Aoz'3, corresponding to the absence of excited molecules in the crystal, 
with a probability calculated per second of 


2 
(Ox'[Plie) = FF [[On' ell? (Ue — Ue), (4.17) 
or into the state A{#’(k’)®f(k’) with a probability (per second) o 
(Eee LP ie) = 5 lege gee (Uae — UE), (4.18) 


Naturally, the probabilities of the converse transitions will be determined 
by the same formulas. Thus, in the crystal, molecular excitation energy 
may be spontaneously transformed into lattice vibration energy without 
a change in the total energy of the system (U{* = U{##') or into exci- 
tation energy with other values of » and of the wave number k of the 
exciton, also with the conservation of the total energy of the system 


(Ute = see), 


The results obtained above are correct only when the probabilities of 
the transitions [(4.17) and (4.18)] are sufficiently small. Only under 
these conditions is it possible to solve equation (4.5) by proceeding from 
the zeroth approximation for fixed molecules. 

If, during lattice vibrations (and upon excitation of the molecules), 
the translational symmetry of the crystal is not disturbed, the matrix 
element [(4.16)] and, consequently, also the probability [(4.18)] will be 
equal to zero when k # k’. When, however, k ~ k’, the transitions from 
one zone pu of excited states to anothee zone p’, or font one excited molecu- 
lar state f to another state f’, will be but little probable if, during the 
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vibrations of molecules, the mean displacement, which is 


AR = Rit +7) — RO 
during the time 
_ h 
a ne 
satisfies the condition 


a oF" aR| < |B — BLY, (4.19) 


Similarly, a condition for the small probability of (4.17) is the inequality 
a okt 


AR | « |Ei* — Be), (4.20) 


where AR = Rit +7) — Rit), r= a 

This condition [(4.20)] is especially important and significant, since it 
determines the smallness of the probability for the transition of molecu- 
lar excitation energy to lattice vibrational energy (and for converse 
transitions). 

The probability for transitions of molecular excitation energy to lattice 
vibrational energy (to heat) may also be evaluated in the following 
manner. 

In molecular crystals the quantum of lattice vibrational energy ho, is 
considerably smaller than the molecular excitation energy hwy. If, now, 
the excitation energy of the molecule includes an electronic excitation, 


= ~ 10-3 — 10-4; 


wom 


if there is only an excitation of intramolecular vibrations, 


Or ~ 107! — 1072. 
rg 


When JT = 0, the excitation of the molecule at the expense of lattice 
vibrational energy is generally impossible. At higher temperatures, such 
transitions are also only slightly probable, since they assume a simultane- 
ous transfer to one molecule of the energy of thousands of phonons (or of 
100 to 10 phonons in the case of a solely vibrational excitation of the 
molecule). However, if the temperature is small, then the displacements 
of molecules from the equilibrium positions are also small. 

The matrix elements [(4.14a) and (4.14b)] entering into (4.14) depend 
on the change of wave functions for the stationary molecular states, 
when the coordinates of the molecules change. These expressions will be 
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slowly changing functions of the displacements of molecules from equi- 
librium positions. 

The matrix elements given by (4.14) depend on the matrix elements 
given by (4.14a) and (4.140) and on the wave functions A(R). If, upon 
excitation of the molecules, the molecular equilibrium positions in the 
crystal and the potential energy as a function of molecular displacement 
from the equilibrium position do not change, the wave functions A(R) 
will be orthogonal for the various vibrational states of the crystal. By 
resolving (4.14a) and (4.140) into a power series of displacements, and by 
substituting the series obtained in the integrals which determine (4.14), 
we shall see that, for the molecule to be able to absorb the energy, for 
instance, of 103 phonons, it is necessary either to take the thousandth 
term of the decomposition by powers of displacements or to consider the 
thousandth order of disturbance if we limit ourselves merely to the terms 
proportional to the displacements. 

If the displacements are small relative to the distances between mole- 
cules in the equilibrium position, these probabilities will be small.f In 
order for large displacements to take place, a fluctuating accumulation 
of lattice vibrational energy in a small region of the crystal is necessary. 
At low temperatures, such fluctuations are extremely improbable. 

Thus, at low temperatures, the transformation of lattice vibrational 
energy to molecular excitation energy is practically eliminated. In other 
words, the molecular motion in a lattice that does not contain excited 
molecules may be considered adiabatic; i.e., one can set expression (4.17) 
equal to zero and consider equation (4.11) as the equation of the station- 
ary states of the system. 

The situation changes essentially for transitions of the reverse type— 
transformations of molecular excitation energy to lattice vibrational 
energy. The state preceding the transition is thermodynamically one of 
nonequilibrium. We are here dealing with a large energy fluctuation in 
the region of the excited molecule. 

However, this alone is insufficient to transfer molecular excitation 
energy to lattice vibrational energy. Because of resonance interaction, 
the excitation energy will migrate from one molecule to another. There- 
fore, the displacement of molecules from equilibrium positions (corre- 
sponding to the nonexcited crystal) to new equilibrium positions may 


+L. D. Landau*! took the general case of a system consisting of a quantum part 
and of a part whose state may be quasi-classically described as being in a potential 
field determined by the energy of the quantum part of the system. Landau showed 
that in this general case the integral, like the integral given by (4.14), is exponentially 
small if, when the coordinates for the system’s quasi-classical part change, one does 
not obtain intersection points for the polydimensional potential curves, which deter- 
mine the motion for this part of the system in the corresponding states. 
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attain large magnitudes only if the excitation is capable of being delayed 
on one molecule for a time sufficient for the molecular displacement (this 
time is of the order of a half period of molecular vibrations). 

In this last case, even at low temperatures, the displacements may 
attain significant magnitudes. 

As a result, the probability of multiquantum transitions increases, even 
when we limit ourselves to the harmonic lattice vibrational model, which, 
for large molecular displacements, is incorrect. 

An exact quantitative evaluation of the probability for such transitions 
is extremely difficult, since one cannot use perturbation theory and since 
one must allow for the anharmonicity of lattice vibrations. 

Thus the conditions for the applicability of the adiabatic approxi- 
mation to the problem of energy states in molecular crystals differ essen- 
tially from the conditions for the applicability of this approximation to 
other problems (with a small number of degrees of freedom). 

The major work of Born and Oppenheimer®® gives the quantum- 
mechanical basis for the applicability of electronic potential-energy 
curves to determining nuclear movement in molecules. In this work, 
it was shown that the adiabatic approximation is accurate when the 
square root for the ratio of the mass of the electron to the average mass 
of the atoms making up the molecule is small. Approximately the same 
criterion is indicated by Pauli (Ref. 63, p. 143): ‘The nuclear vibrational 
frequencies are small in comparison to electronic frequencies.” 

In the case of molecular crystals (and certain other systems with a 
large number of degrees of freedom), these conditions governing the 
accuracy of the adiabatic approximation—i.e., the possibility of deter- 
mining (1) intramolecular (electronic excitation together with intra- 
molecular vibrations) states for fixed molecular positions and (2) the 
motion of molecules in a potential field, which constitutes intramolecular 
energy as depending on the positions of molecules in the lattice—are only 
necessary but by no means sufficient. 

In molecular crystals, the ratio of lattice vibrational frequencies to the 
frequencies corresponding to the change in intramolecular state is very 
small. In spite of this, the adiabatic approximation may be utilized only 
when the temperatureft of the crystal is low enough so that the average 
thermal energy is less than the intramolecular excitation energy. 

When analyzing the states of a crystal with excited molecules, even at 
a temperature of absolute zero, the adiabatic approximation cannot be 
applied. 


Tt When the temperature rises, there is an exponential increase in the probability of 
obtaining displacements, corresponding to the intersection points of polydimensional 
potential curves, determining the motion of molecules in a crystal for two correspond- 
ing states. 
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Such excited crystal states are not rigorously stationary ones, since a 
spontaneous redistribution of the crystal’s energy takes place; the energy 
of the intramolecular excitation either fully or partially transforms into 
lattice vibrational energy. 

Only in the case in which the time for the transfer of excitation to the 
neighboring molecule is small by comparison with the time for molecular 
displacement to new equilibrium positions will the connection between 
the intramolecular excitation and the lattice vibrations be small, and the 
probabilities given by expressions (4.17) and (4.18) will be near zero. The 
corresponding excited states will be almost stationary (quasi-stationary). 

However, it must be kept in mind that this stationary condition is 
disturbed above all for energy transfers, which are linked with a change 
in the wave number of the exciton (within the zone of excited states). 

This nonstationary condition leads to two manifestations: (1) the 
broadening of levels, which is connected with the shortening of their 
lifetime, and (2) the change of velocity for the transfer of excitation 
through the erystal. If the velocity of the movement of excitation 
decreases, there is a significant increase in the probability for a non- 
radiative transfer of excitation energy to lattice vibrational energy; 1.e., 
the nonstationary condition increases in proportion to the transformation 
of excitation energy to heat. 


17. Free and “Localized” Excitons 


As was indicated in the preceding section, two limiting cases must be 
distinguished when analyzing the energy states of molecular crystals. 

1. The excitation transfers from one molecule to another so quickly 
that, under a change in the forces of interaction between neighboring 
molecules upon excitation of one molecule in a crystal, displacement of 
molecules to new equilibrium positions does not have time to occur. 

2. The excitation transfers from one molecule to another so slowly 
that the molecules do have time to occupy new equilibrium positions. 
A local deformation, which travels through the crystal together with the 
excitation, arises in the crystal. 

In case 1, at low temperatures, the “adiabatic” approximation may be 
used with certain reservations; in case 2, this, generally speaking, cannot 
be done. 

The excitation transfer time 7:, as we shall see later, is determined by 
the matrix element Mf, na {(3.10)], by the band of excited states, by the 
wave vector of excitation waves, and by the geometry of the unit cell. 
Thus it depends on the character of the molecular excitation, on the dis- 
tance between molecules, and on their mutual orientations. 

The displacement time rz of molecules from the old equilibrium posi- 
tions depends on the change in the forces of interaction of a molecule 
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with neighboring ones upon excitation of the molecule, this magnitude 
being characterized by a matrix element D£, [(8.13)]; on the mass of the 
molecule; and on the latter’s moment of inertia. For large molecules, 
the time rz may be less than 7;, even if Df, ~ M1, ns. 

In molecular crystals, apparently, the fulfillment of the inequality 


Ta > Te (4.21) 
corresponding to case 1 above, as well as the fulfillment of the inequality 
Ta < Tt (4.22) 


corresponding to case 2, are both possible. 

If the inequality given by (4.21) is fulfilled, the excitation will be 
transferred to the neighboring molecule of the crystal with local lattice 
deformation. This case will correspond to a nonequilibrium state of the 
lattice. The motion of excitons which is not accompanied by a local 
lattice deformation will be called free-exciton motion. 

If the inequality given by (4.22) is fulfilled, the lattice passes by means 
of a local deformation into a new equilibrium state. The local defor- 
mation will accompany the movement of the exciton. This crystal state 
will be called a “localized”? exciton. One must, of course, keep in mind 
that this exciton is not localized in the full sense of the word. However, 
the excitation waves corresponding to the ‘localized’ exciton possess a 
velocity which is small compared with the velocity of free excitons. 

In molecular crystals and under certain conditions, the changeover of 
free excitons to “localized” ones (and vice versa) is possible. 

In the following sections we shall consider the case in which free 
excitons appear in a crystal (when T ~ 0). Here, as we have seen, the 
corresponding states may be approximately considered quasi-stationary. 


18. Equations Determining the Motion of Molecules in a Crystal 


If, in expression (4.11), we substitute the value #9, where 


ES = [ ®$*H.8§ dr = > ES, + V(R), 
we obtain an equation determining the wave functions for molecular 


motions in a crystal when all molecules are not excited. 


[Tr + V(R)]Aoz = 82 Ace (4.23) 
where &2 = Uo — yy EN (4.23a) 
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is the energy of rotational and translational molecular vibrations in a 
crystal which does not possess any excited molecules. 


V(R) = 6 YY [ f le8a(L)[2V namale®(II)|? dr 
na,mBp 


— J f PFD e2a ID Vaamoe25(I1)e29(1) dr]; (4.24) 


here ¢ne(I) and gmg(II) are the wave functions of the molecules na and 
mB for a certain disposition of the electrons in the molecules, ¢,.(II) and 
gmp(1) are the same functions upon permutation between the molecules 
na and m6 of any pair of electrons, and the sum over u extends over all 
possible permutations of electron pairs between molecules na and m8. 

Equation (4.23) represents an equation of translational and rotational 
molecular vibrations in the crystal lattice, with potential energy V(R). 
The coordinates R (displacements and rotations of molecules) are meas- 
ured from the equilibrium values corresponding to the minimum potential 
energy V(R) in a nonexcited crystal. Decomposing the potential energy 
V(R) into powers of displacements from the equilibrium positions and 
limiting ourselves to the terms which are quadratic relative to these devi- 
ations, we can, by introducing normal coordinates, go over to a system of 
independent equations. Each one of the latter will represent an equa- 
tion of a harmonic oscillator. Consequently, the solution of equation 
(4.23), corresponding to the energy 


= » hag;(N oj + 44), (4.25) 
q 


may be written in the form 


Ave = [[ ¥%., (4.26) 
a 


where Wy, is the wave function of the harmonic oscillator. The vibra- 
tional state x of a crystal in this case will be determined by the assign- 
ment of numbers N,; characterizing the number of phonons which possess 
a wave number g and a polarization j(j = 1,2, ..., 6c). In Chapter 
III, we determined the wave function [(3.7)] for the excited state of a 
crystal in which one molecule is in the fth excited state and all others 
are in a normal state for a certain disposition (characterized by the 
coordinates &) of the molecules in the crystal: 


B= (GN)-* Y dnaxhas (4.27) 
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and the energy for the corresponding state of the crystal is 


by Anat» ma( Fe) 
E!(R) = ye 2 + AB, + V(R) + Di(R) + ™ 


AmB 
(4.28) 
After substituting (4.28) in equation (4.12), we find the equation which 


determines the motion of molecules in a crystal possessing one molecule 
excited to the fth state: 


ou Mo. ma( BR) Ona 
Ir, + V(R) + Df ais + aia PO head Af = 8fAf (4.29) 
where 2) Ee, — AES, (4.30) 


is the energy of molecular vibration in a crystal possessing one excited 
molecule, AE4, = EL, — EX. is the energy for the excitation of one mole- 
cule, Df,(R) has the same meaning as (3.13), and the matrix element 
Mims coincides with (3.10). 

The matrix elements Dig(R) and M1, , in equation (4.29) of the 
excited crystal play the role of supplementary energy to the potential 
energy V(R). Asa result of the appearance of these terms, the equa- 
tions for molecular motion in an excited crystal differ from the equations 
for motion in a normal crystal. The frequencies of normal vibrations 
change, and the equilibrium positions are displaced in relation to the fre- 
quencies and to the equilibrium positions in a nonexcited crystal. 


CHAPTER V 


FREE EXCITONS IN A MOLECULAR LATTICE 


19. Zeroth Approximation 


Let us assume that the inequality (4.21) for a given excited state of a 
crystal is fulfilled and that the crystal is at absolute zero. 

The motion of molecules in the excited state of a crystal will be deter- 
mined by equation (4.29). Ina nonexcited crystal, the molecules execute 
zero-point vibrations near the equilibrium positions (R = 0), determined 
from the minimum of the potential energy V(R). 

The matrix elements D/, and M‘, ..s, which enter into equation (4.29) 
and are functions of the molecular positions in the crystal, may be decom- 
posed into power series of the displacements R relative to the equilibrium 
positions of the molecules (in a nonexcited crystal): 


Dig(R) = Dig(0) + Dig(1) + - °°, 
Min ma( FR) vd ME ame(0) + ME, ma(1) + er 


aDf 
t = — za! mB 
Di,(1) ) Re (32a) 


toy! 


pi aM? (5.2) 
if er za’ na,mp : 
Moa mp(l) > Rx ( OR’ , 


ana’ 


and D5,(0) and M/, ng(0) are the values of the matrix elements Df, and 
MM‘, mp for the equilibrium configurations corresponding to the nonexcited 
erystal; x runs through six values corresponding to the six degrees of 
freedom of the molecule. 

At low temperatures, which are of principal interest here, the displace- 
ments of molecules from equilibrium positions are small in comparison 
with the distances between molecules. Actually, at a temperature of 
absolute zero, the average quadratic deviation from the equilibrium posi- 
tion will be equal to 


(5.1) 


where 


Lh 
2M iw: 


for the 7th vibrations; w; is the circular frequency of the corresponding 
vibrations. Thus, for the magnitude of the molecular displacement 
71 


2 
to = 
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from the equilibrium positions in a naphthalene crystal (M = 128, 
a ~ 6- 10#"), we obtain 
Zo ~ 6-107 cm. 
Therefore, in the zeroth approximation, only the first terms may be 
left in (5.1); thus equation (4.29) assumes the form 


DY Moma) On 


ln, + V(R) + Df,(0) + caer ee coma Aff = 8fAl, (5.3) 


Then, in (5.3), the expression 
eT Mi oma (0) Gna 


AmB 


will have the dimensions of energy and will not depend on R. If we 
equate this expression to e, we obtain the equation for determining this 
energy: 


~ame + JY Mia mp(O0)@ne = 0. (5.4) 


Equation (5.3) may now be written 
[Te + V(R)JAL = (8% — DL,(0) — «) AP. (5.5) 


Comparing this equation with equation (4.23) for the translational- 
rotational vibrations of a crystal, we are satisfied that 


Af = Aj,, (5.6) 

& = & — Df,(0) —«. (5.7) 

Thus in the zeroth approximation, as was to be expected, the excitation 
of the crystal is not accompanied by a change in the vibrational state of 


the crystal. 
At the same time, the excitation energy of the crystal, 


AE! = UL — Una, 


may be calculated with the aid of expressions (4.23a), (4.30), and (5.7), 
and is equal to 


I 


AE! = AEL, + Df,(0) + &, (5.8) 


where AES, is the excitation energy of a separate molecule; ¢/ is deter- 
mined by a system of homogeneous equations (5.4) coinciding with the 
system of equations (3.9), which were studied in Chapter III. Conse- 
quently, in the zeroth approximation the excitation energy of the crystal 
already differs from the excitation energy of one molecule by two terms 
Df ,(0) and e{(k). 
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Thus the crystal excitation energy obtained here in the zeroth approxi- 
mation coincides with the excitation energy of the crystal in the case of 
rigidly fixed molecules in the lattice points. Consequently, the results 
of Chapter III are correct if, upon excitation, free excitons appear [the 
inequality given by (4.21) is fulfilled] and if we can limit ourselves to the 
zeroth terms in the decompositions [(5.1)]. In Chapter III, we analyzed, 
with the aid of group theory, the properties of the solutions obtained. 
Here we shall again analyze some of the previous results, without, how- 
ever, resorting to group theory. The results we shall obtain here are in 
some respects clearer, and may in principle yield not only qualitative, 
but also quantitative, data on the spectrum. 

Thus the excitation energy of a crystal, 


AE{(k) = ABZ, + Dig(0) + i(k), (5.9) 


will depend on the wave number k of the excitation waves (excitons). 
Furthermore, to one excitation of a molecule in a crystal, there will 
correspond o quasi-continuous bands of excited states, which, generally 
speaking, may overlap. To each one of these bands of excited states in 
a crystal there will correspond its own wave function, which depends on 
the wave number k. 


1 
ViN(h, 1, B) = ALB = Tay ARR) Y) eBExL, — (610) 


where B# are the solutions of the system of equations (3.15) for the 
corresponding value «¢,(k). 

The normal state of a crystal, from which we measure off the energy 
terms, is described by the wave function 


Wi(r, R) = Ape, (5.11) 
where ®° is determined by (3.2). 


20. Energy Terms of Several Molecular Crystals 


Let us consider the solution of equations (3.15) and (3.17) for the 
particular case of crystals of the monoclinic-prismatic system in the 
spatial symmetry group C3, with two molecules in the unit cell (¢ = 2). 
In crystals of this type, the molecules are at the centers of symmetry 
of the crystal; therefore, the matrix L.s(k) will be real and symmetrical. 

Equation (3.17) will assume the form 


Luk) — « Lye(k) | _ 


Lox (k) Bgcy =e" 
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Since Ly,(k) = Loo(k) and Ly. = Dei, we obtain two solutions: 


e(k) = Ly(k) + Li(k), 


2 
ex(k) = Ly(k) — Li(k). (5.12) 
Substituting the values Leg from (3.16), we obtain 
¢,(k) = » Maroc + » Mn2,01e™. (5.18) 
n=O n 


If » = 1, the plus sign must be taken, and if u = 2, the minus sign. 
Substituting (5.12) in equations (3.15), we obtain 

By 1 ifp=1 
Be | 1 if p = 2. 


Therefore, to the two energy zones of the excited states [(5.13)], there 
will correspond two systems of wave functions: 


Wi (Ie, r, R) = (2N)-4AS,(R) D (xhs  xhodei. (5.14) 


The matrix elements Mnroms for the transfer of excitation from one 
molecule to another diminish with the distance not more slowly than 
Rz3me- However, when calculating (5.13), one cannot limit consider- 
ation to interaction with the nearest molecules, which surround the 
excited molecules, since, with the increase of the distance, the number 
of molecules interacting with a given molecule also increases (as R?). 
At the same time, we must keep in mind that the matrix elements essen- 
tially depend on the mutual orientation of induced molecular dipole 
moments. Therefore, interaction with a more remote molecule will occa- 
sionally be greater than with a close one. The signs of the matrix ele- 
ments are also determined by the orientation of the molecules. 

For a rough estimate of the magnitude of (5.13), one may use the 
approximate expression 


ex(k) = 2{M, cos ka + M, cos kb + 2Ma 
+ [M, 4 2(Mae + Ma)] coske + -- -}. (5.15) 


Here M, is the matrix element for the transfer of excitation to the nearest 
(first) molecule situated on the a axis of the crystal (see Figure 9). Fig- 
ure 9 shows the disposition of the molecular centers of gravity in a mono- 
clinic crystal with two molecules in a unit cell. The molecules having 
an orientation identical with that of the first molecule are designated by 
shaded circles. The molecules (second) having an orientation which is 
different from the orientation of the first are designated by white circles. 

The matrix element M, in (5.15) determines the transfer of excitation 
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to the nearest (first) molecule, located on the b axis. The matrix ele- 
ment /, is for the transfer of excitation to the nearest (first) molecule, 
situated on the ¢ axis. The matrix 
element M., is for the transfer of 
excitation to the nearest second mole- 
cule located in the center of the unit 
cell face ab. The matrix elements 
Ma, and M‘,, are for the transfer of 
excitation to the nearest second mole- 
cules, located at the centers of the 
faces ab, which do not coincide with 
the plane, where the primary mole- 
cule is. 

If in the unit cell there are more 
than two molecules, a direct solution 
of the system of o equations (3.15) 
and (3.17) isdifficult toachieve. For 
crystals in which the molecules 
occupy the centers of symmetry, the 
unknowns B# will be real and equal 
aM any an absolute magnitude. a Fie. 9. Disposition of the molecular cen- 
is then easy to write directly ¢ in- ters in a monoclinic crystal with two 
dependent and orthogonal solutions molecules per unit cell. 
for this system of equations. Thus, 
for crystals of this type containing four molecules per unit cell, the 
values Bt may be written in tabular form: 


J 1 1 1 
—1 —1 

rn ee a ea (5.16) 
i Moe ee 


Here the number of a row corresponds to the index » and the number of 
a column to the index a. 
Using the values of B# obtained, and utilizing (8.15), we can write the 
solutions for e,(k): 
+ + + 
e(k) = > Marne + | _ Maron is Mn3,01 4 Mnaor |e, (5.17) 
nO n 
fe = = 
where uw = 1, 2, 3, 4. 
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To each energy zone y, there will correspond its own wave function: 
+ + + 
1 a + = 
WL 9, B) = ee ABR) » ha xhe T xb 7 oxhe | et. (6.18) 
: + 


1+ 


The deduced system of terms of a molecular crystal does not determine 
the absorption spectrum. We still must determine the selection rules, 
i.e., the rules indicating those energy-zone sublevels to which the tran- 
sitions under the action of a light wave are possible. We shall do this 
in the following section. 


21. Transitions under the Action of Light 


The electric field intensity (in a crystal) of a plane electromagnetic 
wave falling on the crystal is 


E = Epei(er-e), (5.19) 


Here Ep is the amplitude, w is the vibration frequency multiplied by 27, 
and Q is the wave vector. Under the action of this wave, the crystal 
will pass into an excited state. 

To calculate the probability of the transition under the action of the 
electric field [(5.19)], we must add the operator for the perturbation to 
the time-dependent equation of quantum mechanics; this yields 


weit = — & Ee (or—wt) grad. (5.20) 


Here ¢ is the electronic charge, u is the electronic mass, and h is Planck’s 
constant divided by 27. 

Then the probability (Pf) for the transition in one second from a 
normal state of the crystal Y° to an excited one W“ will be proportional 
to the square of the matrix element for the operator w from (5.20): 


2 
Ph = F |(uflwl0)|? 8(AE{ ~ he), (5.21) 
where (uf|wlO0) = f¥Le*wwW? dr dR. (5.22) 


In the matrix element [(5.20)], the integration is carried out over all 
internal coordinates of the molecules dr and the coordinates determining 
the spatial distribution of molecules dR. 

Substituting in (5.22) the value w from (5.20), and taking into con- 
sideration that 


iE wie grad W° dr = — © : WHA dr, 
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we obtain 


(uf|wl0) = eBof Wie*e!(r) Wo dr dR. (5.23) 


We shall be concerned with light absorption in the infrared, visible, and 
ultraviolet parts of the spectrum, where the wavelength is significantly 
greater than the dimensions of the unit cell of the crystal. Therefore, 
we shall disregard the phase difference between different molecules in 
one unit cell. In other words, we shall substitute for e'@ at the point 
where the electron is, the corresponding value e*@ at a certain point of 
the elemental cell of number 7, in which one of the molecules becomes 
excited, 

Let us now substitute into (5.23) the values for the wave functions 
and W/ from (5.10) and (5.11). If we consider also that in the zeroth 
approximation there is no coupling between the molecular vibrations in 
the lattice and their optical excitation (consequently, the integration over 
the external molecular coordinates may be carried out independently), 
we obtain 

(uf|wl0) = (6N)-* ¥' e'(0-! Ey Y Be f xttrb? dr. (5.24) 
n a 


lf we take into account the phase difference of the light wave within the 
limits of the unit crystal cell, the matrix element [(5.24)] becomes some- 
what complicated: 


(uf|w|0) = (oN)-¥ Yee ) Buel f xftrb°dr. (5.25) 


In our case Qra «1; therefore, we shall utilize the simpler expression 
[(5.24)]. 

In the integral entering into (5.24), the integration is carried out over 
all internal coordinates of the molecules. By virtue of the translational 
symmetry of the crystal, which in our case (free excitons) is not dis- 
turbed upon light absorption, these integrals will not depend on n. 
Furthermore, substituting the values x, and S°, we obtain 


Sxftr6° dr = oltre’, dr = [oftres dr. 
Now, (5.24) may be written 
(uf|w|0) = (0N)—%eEo (Be f oltre? dr) Y ee», (5.26) 


The summation over 7 always yields zero if the following conditions are 
not fulfilled: 
0, (a) 


22'S uaa Goiee. 1 eae 


The az? are the vectors of the reciprocal lattice. 
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The selection rules [(5.27)] express the law for the impulse conservation 
of the light wave and of the excitation waves produced, as well as the 
possibility for transmitting this impulse to the lattice as a whole (case b). 

For infrared, visible, and ultraviolet light, we can neglect the wave 
vector for the light wave (Q ~ 0); then the selection rules [(5.27)] will 
assume the form 


0, 
k = | oe (5.28) 


If conditions (5.27) are fulfilled, then 


oc 


(uflw|0) = Je eEy » Be i gitree dr. (5.29) 


a=l 


Expression (5.29) determines the matrix element for the transition to an 
excited uf state of the crystal and shows that in the first approximation 
such a transition takes place only when the matrix element for the dipole 
operator of a separate molecule is not equal to zero, i.e., when the corre- 
sponding transition is permitted in the molecule. If the dipole transition 
in the molecule is forbidden, then in the dipole approximation the tran- 
sitions to all yu excited-state zones are also forbidden. In this case, it 
becomes necessary to use (5.25) and to take into account the change in 
light wave phases within the limits of the unit cell (we formerly dis- 
regarded this change) in order to obtain the probabilities for the tran- 
sitions of a higher multipole type. These probabilities will be less than 
for dipole transitions by an order of magnitude 


me dimension of the molecule \”# 
— light wavelength 


(5.30) 


The order of multipole character is / for the corresponding transition 
in a molecule. 

We must, however, keep in mind that the distances between separate 
uw zones are also proportional to (5.30). These distances are determined 
by the matrix elements for the resonance interaction of translationally 
nonequivalent molecules in a crystal. In other words, for forbidden 
molecular transitions, the crystalline structure is less significant. 

If the dipole transition in a molecule is permitted, the transition in a 
crystal, generally speaking, will also be permitted. However, depending 
on the mutual orientation of electric moments for intramolecular tran- 
sitions in molecules in a unit cell, and depending on the polarization and 
direction of light propagation relative to the crystallographic axes of the 
crystal, transitions will take place to various zones obtained by the 
splitting of molecular transitions in the crystal. 
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Thus, upon light absorption (zeroth approximation of the theory), an 
excitation does not occur to any arbitrary level of the quasi-continuous 
excitation zones of energy: 


Bt : 
ef (k) = > By e Mi, nt). (5.31) 
1 n 


On the contrary, transitions occur only to those levels that correspond to 
the value of the wave number k, satisfying the conditions given by (5.27) 
or the approximate conditions given by (5.28). Simultaneously, we note 
that the energy values turn out to be identical, not only when the con- 
ditions for case a in (5.27) are fulfilled, but also when conditions for case b 
are fulfilled. 

The results obtained here pertain to a case in which, at the moment of 
light absorption, free excitons arise in a crystal. For this to occur, it is 
necessary for the group velocities (proportional to k, for small k) of the 
excitation waves to be sufficiently large. 

Upon excitation of excitons by light (k = Q), this condition may be 
violated for certain values of the coefficients Bt. The corresponding 
bands of excited states in a crystal will not be realized. 


22. Transitions under the Action of Light in Crystals Containing Two 
Molecules per Unit Cell} 


Among the molecular crystals of this group are the naphthalene and 
anthracene crystals (the structures of which have been thoroughly 
studied**-*5), as well as a series of other crystals. It is the custom in 
this case to designate the vectors (a1, a2, and as) of a unit cell of a crystal 
by the letters a, b, c, respectively. The vector determining the position 
of two molecules in a unit cell of a crystal is given by 


ri = 0, 


_a lb (5.32) 
ee ae 


If we use the functions given by (5.14) for the excited states of the crystal, 
we obtain two matrix elements for the transitions to an excited state: 


ciflwjo) = [Rete ( f etteyt dr +f oftrot dr), 
W (5.33) 
(2f|w|0) = Va eEy (J of* re? dr — i of* res ar). 


+ The contents of Sections 22 to 24 have been published in part previously. 
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Each of them will have its own corresponding excitation energy: 


AE{ = AE fs + Dhg(0) + YY (Mir or + Mhr ore’, 


(5.34) 
AEi 


AEfs + Dig(0) + YY (Mésor — Mioode'. 


Here Mato, is the matrix element for the transfer of excitation from the 
first molecule of the zeroth unit cell of the crystal to all first molecules 
located in other unit cells; M{,2,0, is the matrix element for the transfer 
of excitation from this same molecule to all second molecules. 

If both molecules in one unit cell have an identical orientation, only 
one of the matrix elements [(5.33)], corresponding to the plus sign, will be 
different from zero for the corresponding direction of Eo. Therefore, in 
such a case, to each transition in an isolated molecule there will corre- 
spond one transition in a crystal; this latter transition will merely be 
displaced relative to the molecular transition by the magnitude 


h Aw = Dig(0) + >! (Maio + Mie ee. 


If, however, both molecules possess different orientations, both matrix 
elements [(5.33)] may be different from zero. In this case, to each tran- 
sition in a molecule there will correspond in a crystal two absorption 
bands determined by (5.384). The positions of these bands will, in 
general, depend on the direction of light propagation in the crystal 
(direction of vector Q). The emergence of one or another of these 
bands, or simultaneously of both, will depend on the polarization of the 
light wave vector Eo. 

Crystals of the type studied are usually produced in the form of fine, 
thin plates whose planes coincide with the cleavage plane of the crystal. 
This plane contains the monoclinic axis b of the crystal and the a axis 
perpendicular to it. When this experiment is conducted, the light ray is 
usually directed perpendicular to this plane of the crystal. 

Let us assume that the light falls perpendicularly onto the face ab of 
the crystal and that the projection onto this face of the electric vector 
for the transition moment of a molecule forms the angles a and —a 
(a < 45°) with the b axis for the two molecules. This case is realized, 
for example, for the electronic transition in a naphthalene and an anthra- 
cene molecule when the vector for the transition moment is parallel to 
the intermediate molecular axis. Then, with the aid of (5.33), we see 
that the frequency corresponding to the energy AE{ given by (5.34) will 
be absorbed if the electric vector for the light wave is parallel to the 
b axis of the crystal. The frequency corresponding to AE of (5.34) will 
be absorbed also if the light wave vector is parallel to the a axis of the 
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crystal. The intensity ratio of both of these absorption bands is given by 


ay = cot? a. (5.35) 


The distance between these absorption bands, measured in reciprocal 
centimeters, is 


aX 
Av=z > Mo, ef. (5.36) 


Here c is the velocity of light. 

When the projection onto the face ab of the transition moment vector 
in the molecule forms angles y and —y (y < 45°) with the a axis for two 
molecules, respectively, the frequency AE{/h will be absorbed if the elec- 
tric vector of the light wave is parallel to the a axis. The frequency 
AE{/h will be absorbed also if the electric vector of the light wave is 
parallel to the b axis of the crystal. The intensity ratio of both bands is 


I(1) 


T(2) = cot? Y- 


This case applies to an electronic transition in the naphthalene and 
anthracene molecule when the vector for the transition moment is 
parallel to the long axis of the molecule. 

Knowing the dependence of the crystal excitation energy [(5.13)] on 
the wave number k for excitation waves produced, we can determine the 
group velocity for the propagation of these waves in the crystal: 


U,(k) = = grads e,(k). (6.37) 


Of course, for zones possessing different values of » and belonging to one 
excited state of the molecule, the group velocities will be nonidentical. 

If the light falls perpendicular to the face ab of the crystal, then for 
excitons appearing under the action of a light wave and by virtue of the 
selection rules [(5.27)], 


ka = kb = 0, ke = Qc sin 8, 


where @ is the angle between vectors a and c of the crystal unit cell. 
The group velocity of free exciton motion will be (Qe sin 8 « 1) 


U(Q) ~ — 2@er sin PIM. fu SM ga AON ga ick 295, 
(5.38) 
U(Q) = — 20e' sin’ B iat, = OM gre Mina bea, 
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These two expressions correspond to the two values for the excitation 


energy [(5.34)]. 
The ratio of these velocities is 


Ui, ay M, + 2(M ae + M ax’) + # SS a 
U» M, = 2(M ave + Mw’) + 6S 


Depending on the sign and on the relative magnitude of the matrix ele- 
ments entering into (5.39), this ratio may assume different values and 
may be at times substantially different from unity. 

Excitons moving with a lower velocity will interact with the lattice 
more strongly; such states will be less stationary. In certain crystals 
the velocity of travel for excitons, corresponding to one of the excitation 
energies [(5.34)], may be so small that the inequality given by (4.22) 
will be realized; this leads to the formation of “localized’’ excitons. 

We see from (5.38) that the group velocity essentially depends on the 
dimensions and form of the unit cell and, in particular, on the mono- 
clinic angle 6. In lattices with angle @ close to 90°, other conditions 
being equal, the group velocities are extremely high. Therefore, the 
light absorption bands in such lattices must be narrow. The narrow- 
ness of the absorption bands for the phenanthrene crystal!® is possibly 
explained by the fact that, in this crystal, angle 6 is close to 90°. 

If the light falls perpendicular to the face ac of the crystal, and if the 
projection onto this face of the electric moment vector for the molecular 
transition forms an angle 6 for both molecules with the direction of the 
erystal’s a axis, only the frequency corresponding to the excitation energy 
AE\ from (5.34) will be absorbed. This case must apply to the naphtha- 
lene crystal. It would be very desirable to conduct a special experiment 
to investigate the absorption spectrum when the light impinges on a 
monocrystal in the foregoing manner. 

If the molecular transition projections onto the face ac of the first and 
second molecule in the unit cell form various angles with the a axis, we 
can again observe two absorption frequencies corresponding to one tran- 
sition in the molecule. 


(5.39) 


23. Transitions under the Action of Light in Crystals Containing Four 
Molecules per Unit Cell 


Among the crystals of this type are phenanthrene, azobenzene, and 
others. In crystals of the monoclinic system, four molecules usually 
occupy sites determined by the radius vectors 


_ an _a,b _¢ _a_b c 
rea(n=On=F+ on = gu $48 +5) 


If we use the coefficients given by (5.16), we obtain, with the aid of 
(5.29), the following matrix elements for electronic transitions under the 


FREE EXCITONS IN A MOLECULAR LATTICE 83 


action of light: 


N = 
(uf\w|0) = ae cE J oltre? dr - ghtre3 dr 


tf Te: a5 fon dr }- (5.40) 


For each of these transitions, the corresponding energy for the exci- 
tation of the crystal is 
AEI => AES, + Df £,(0) 
+ + + 
+3 | Moi Mose 7 se 5 
n 
+ — i 


ML, o, |e. (5.41) 


Depending on the orientation of molecules in a unit cell and on the direc- 
tion and polarization of light, one will observe one or another of the 
four absorption bands corresponding to the excitation energies [(5.41)). 
Knowing the arrangement of molecules in a unit cell as well as the direc- 
tion of the electric moments for the transitions in these molecules, we can 
obtain from (5.40) qualitative and quantitative data on the number of 
split components, on the magnitude of the splitting, and on the displace- 
ment of the absorption bands in a crystal for various polarizations and 
directions of light propagation in the crystal. 


24. Molecular Theory of Pleochroism 


In white light, certain crystals of both organic and inorganic com- 
pounds appear colored, with the color depending on the direction of 
propagation and polarization of the impinging light waves. This phe- 
nomenon is called pleochroism or polychromism. To a larger or lesser 
degree, it characterizes all double-refracting crystals; however, for many 
crystals, the phenomenon is observed in the spectral regions not per- 
ceived by the eye. 

In phenomenological theory, pleochroism is described by introducing 
the tensor for dielectric constants and the tensor for electric conductivity 
for each light frequency. As far as we know, a molecular (nonphenome- 
nological) theory of pleochroism does not exist at the present time. 

The results of the theory for light absorption by molecular crystals 
set forth in this chapter, and in Chapter III as well, may be considered 
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as a first attempt toward forming such a theory for molecular crystals. 
As we know, precisely for molecular crystals, the phenomenon of strong 
pleochroism is little understood, even when the molecules of the crystal 
possess high anisotropy so that the light absorption of a given frequency 
may take place only if the projection of the light wave electric vector 

onto a given molecular axis differs 


from zero. Actually, in molecular 
oe crystals containing more than one 
molecule per unit cell, the molecules 
are found to have a diverse orientation. 


! It is therefore often impossible to find 
| such directions for the light wave elec- 
tric vector as would produce a zero 
value projection onto the correspond- 
ing directions for all molecules in a 
| unit cell. This may be easily seen, 
| | for instance, in Figure 10, which sche- 


matically depicts a typical disposition 

of molecular projections onto one of 

= the crystallographic planes. If a 
molecule absorbs light of frequency » 


. ait Ree a a when the light wave electric vector is 
ae pa ee ere cevrtillle. parallel to the long (or short) molec- 
graphic plane. ular axis, the intensity for the ab- 

sorption of light of frequency », 
although changing with polarization of the impinging light, apparently 
will not be equal to zero. 

In his monograph, “Structure of Molecules’ (Ref. 67, p. 161), M. V. 
Volkenstein states: ‘‘The question always arises as to whether the differ- 
ence observed in absorption spectra along and across the crystal axis 
is a characteristic of the crystal as a whole, occurring upon interaction of 
molecules in a lattice, rather than a characteristic of separate molecules.” 
In this work, we have taken into account the resonance interaction 
between molecules, which occurs at the moment of light absorption. 
This in turn provides an answer to the above question and leads to 
qualitative results; with the aid of these results, we can, in a uniquely 
defined way, solve the problem concerning the number of split compo- 
nents in a crystal for each molecular transition, as well as their polari- 
zations, if the number and disposition of molecules in the crystal unit cell 
is known. In addition, if the wave functions for the normal and excited 
states of the molecules are known, we can also determine the magnitude 
of the splitting, i.e., quantitatively calculate the pleochroism, of such a 
crystal. 


CHAPTER VI 


EXCITATION OF FREE EXCITONS SIMULTANEOUSLY 
WITH THE EXCITATION OF LATTICE VIBRATIONS 


25. Solution of Equations Determining the Vibrations of Molecules 
in the Lattice in the First Approximation 


We previously disregarded the connection between excitation waves 
(excitons) and lattice vibrations. Thus, in the decompositions [(5.1)] of 
matrix elements D and M into powers of displacements from the equi- 
librium positions, we utilized only the first terms and found the solutions 
of equation (4.29) for free excitons in the zeroth approximation. 

Let us now consider the terms of first order [(5.2)] relative to displace- 
ments from the equilibrium positions. Then, taking (5.4) and (5.10) 
into account, equation (4.29) may be written in the following form: 


[Tr + V(R) — &{ — Dfg(0) — fk) Al = —Ki,01)AZ (6.1) 
In this, 


Kil) = Dia(1) + > Miamp(1) i aaa (6.2) 


The uth root of equation (3.17) is e,(k). 

The ratios of unknown coefficients B4/B% are determined by equa- 
tions (3.15) for the corresponding values «,(k). 

If we take into account (5.2), as well as the fact that the xth (x runs 
over six values corresponding to the six degrees of freedom of the mole- 
cule) displacement of molecule @’ in the n’th cell may be expressed 
through coordinates a,; (see Section 4) as 


za’ = (gN)—4 > e? (agje' + aj,’ (6.3) 
qi 
wherej = 1,2, ... , 60 and gq, > 0, we obtain 
Ki,(1) = —  Wi,(qin’) (agje" + ake’). (6.4) 
gin’ 
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In the above, 


oop F ; fe) 
Wi,(ain’) = (oN) > oar | ade KL | (6.5) 
Ki,(0) = Dig(0) + > MM a.ma(0) ae (6.6) 
Equation (6.1) may be written in the following form: 
[H + Ki,(1)]Af = 6.Af, (6.7) 
where 

H = Tr + V(R) (6.7a) 
and & = &£ — Df,(0) — «(k). (6.76) 

When Ki(1) = 0, equation (6.7) becomes 
HAP = 8c Al. (6.8) 
Here & = & — Df,(0) — e(k). (6.8a) 


Equation (6.8) was studied in Chapter V [cf. (5.5)]. Its solution deter- 
mines the translational-rotational vibrations of molecules in a nonexcited 
crystal, since by virtue of (5.6) we have Af° = Aj.. The vibrational 
state of the molecules in a crystal may be described by indicating the 
number Ny; of phonons for each value of the wave number q’ and of the 
Bolerieation 7. Then 


Af? = AS, = [| W3,,.(aey”). (6.9) 
q" 


The lattice vibration energy is 


82 = ) Eve (6.10) 
[Nai 
Here Ein,; = hwg( No + ). (6.104) 


The operator K/,(1) in Saionen (6.7) will be considered as a small pertur- 
bation. Then, knowing the solution [(6.9)] of the simplified equation 
(6.8), we can obtain the approximate solution, which allows for the 
influence of K{,(1). 
The energy of the system in the first approximation does not change: 
& = 8&2, (6.11) 
since the average values (i.e., the diagonal matrix elements), relative to 
the wave functions for the simplified equation (6.8), of the perturbation 
operator K{(1) are equal to zero. Besides, by virtue of (6.70) and (6.8a), 
we have 
& = 8, 
where &% is determined by (5.3). 
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The change in energy could occur [for operator AJ,(1)] only in the 
second approximation, or in the first if we consider the higher-order terms 
in the displacements from the equilibrium positions in the decompositions 
{(5.1)]. Inasmuch as the magnitude of displacements from equilibrium 
positions depends on the temperature, in this case, owing to the energy 
change of the system, a displacement of the absorption band, which 
depended on the temperature, would occur. In this section we shall con- 
cern ourselves with low temperatures (T ~ 0). We can then limit our- 
selves to considering only A/,(1) and the first-order perturbation theory. 
In this approximation, as we indicated above, a displacement of the 
absorption band does not take place; however, the wave functions change. 

In the first approximation, the wave functions will be: 


I] Yar) KD I] ¥e1, daev 
An = I] Wrrei( Gai) + » | ah J (I ) i (I a's ) 
a 


Nai gi > (Exe. ona ER ari) 


(6.12) 


In (6.12) the summation » must be taken as follows: 
Nj 


i oe 


Ni N'q15, N'asig Nagi 


Each Nj, runs over whole numbers 0,1, ... , ~. 


J 


Using (6.4), we shall introduce the following designations: 
| (I Wires) KEL) (1 Vie dary) = MG(+) + AG (-). (6.18) 


In the above, 


M,(—) = — 2, Wi. gin'yew [ (Il Ves) Oy (I Vine, dda’), 
AL.(+) = - Ym ty (qjn’ ein’ / (i Vy a) Qgj (II Wrens Ags ). 


Because of the orthogonality and normalization of wave functions Wy,, in 
the integrals, only the integral over the normal coordinates a,; remains 
from the products of these functions. Thus we obtain 


AL (+) = — Yo Whi(gin')e [ VE ,dai Wires Ue, (6.14) 


qin’ 


Me(—) = — Y Whe(gin'e er [Wide irg Ug (6.15) 
gin’ 
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The integral A/,(+) differs from zero only if 


NG = N,— 1; 
then 
; hNg 
f oe S (ain'\eian’ . {it ai, 
ALC +) », Wi.(qjn Je a Jae (6.16) 
qin’ 
If Ni, = Ni, +1, then 
sey age WAN gre LY 
aX, f / iqn eee | SL 
» Whalajnt ort PCa (6.17) 
qin! 


If, while taking into account (6.16) and (6.17), we substitute (6.13) in 
equation (6.12), and if we consider that 


Eo. ~Ee,.) =| ~hew NG = Natl, 
ae Navi Mens) | +hug;, if Nis = Na im 1, 


7" 


we obtain the wave function of equation (6.7) in the first approximation: 


Wi.(gin’) 
Av=]] vs,, Mika) ( gia /N,; [| ve 
1 I] Nai a ee V 2hM gj, (« VN ai ai I] Nqi-\ 
@ qin 


— ena’ Vis + has bt Yan) (6.18) 


If we acknowledge that the distribution of phonons ee to possible 
lattice vibrations qj is, in the case of thermal equilibrium, determined by 
the Bose-Einstein statistics, then 


(Vado = [exp (Fi?) — 1] (6.18) 


For absolute zero (Nqj)avg = 0, and the wave function [(6.18)] adopts a 
simpler form: 


At =|] ve, — ¥ Wen) iw TT ve 6.20 
T= ll sa V 2hM gj; I] ss Se 


a gin! 


26. Absorption of Light at Low Temperatures. Distribution of 
Intensity within an Absorption Band 


The wave function for an excited state of a crystal (or T = 0) may be 
obtained by multiplying (6.20) by the wave function 4: 


Af on Be 
z etn —* yf. 6.21 
a > By X* (6.21) 


vie 
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Here A‘ is determined by (6.20); x4, is determined by (3.5); and B# are 
the solutions for the system of equations (3.15) for the corresponding 
value of u, the latter being the numbers of the zone of excited states. 

The wave function for the normal state of a crystal (for T = 0) will 
have the form [®° is determined by (3.2)]: 


Ve(r, R) = &° [| ¥,,. (6.22) 

q 
The probability of the transition under the action of a light wave 
{(5.19)] from the normal state of the crystal [(6.22)] to the excited state 


[(6.21)] will be proportional to the square of the matrix element for the 
perturbation operator of [(5.20)]. Let us designate this matrix element as 


(uf|wl0)! = [WH*(k, r, Ryw¥?(r, R) dr dR. (6.23) 


Substituting w from (5.20) as well as the wave functions given by (6.21) 
and (6.22), we obtain 


(uf|w|0)? = (uf|wl0)° 
+ > vet Bf Oster VLA) Y eer, 6.24 


Yi.(g/) = =| l] Vii ves I] vs dR, (6.25) 


(ufjwl0)? = Ee i mf Bocods a (6.26) 


Here 


The matrix element given by ees exactly coincides with the previously 
studied matrix element [(5.26)] for the transition without excitation of 
the lattice vibrations (without production of phonons). 

The second matrix element in (6.24) indicates the possibility of molecu- 
lar excitation with simultaneous lattice vibrational excitation (production 
of one phonon). This matrix element differs from zero only if the follow- 
ing conditions are fulfilled: 


0 
SEU ties yg are on 


For the spectral region which we are considering, the wave vector of the 
light wave (Q = 0) may be disregarded. We then obtain the conditions 
for the possibility of excitation of an exciton simultaneously with phonon 
excitation in the following form: 


Me | 2 quote. (6.28) 
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If conditions (6.27) are realized, the matrix element responsible for the 
exciton excitation simultaneously with the phonon excitation will be 
equal to 


cuflwioys, = [Mea () 3 | etteretar)- vaca. 6.29) 


The excitation energy for the crystal will here be equal to 
AET(qs) = AETe + Dig(0) + e{(—a) + hors, (6.30) 


where AES, is the excitation energy for one molecule and D‘,(0) is the 
difference in interaction energy between the excited and normal molecule 
on the one hand and all other molecules on the other hand. This energy 
is calculated for an equilibrium disposition of molecules in a nonexcited 
crystal; hw,; is the energy of the excited phonon: 


e(—q) = Y S Mia mp(O)e*. (6.31) 
a=1 n 


For crystals of the monoclinic system, each containing two molecules per 
unit cell, expression (6.31) can be written in the following form, if we 
limit ourselves to considering the resonance interaction with only the 
nearest molecules surrounding the excited molecule: 


ef(—q) = 2{M.cosqa + M,cosqb + 2Ma 
+ [M. + 2(Mae + Mas')] cosqe + +--+}. (6.82) 


The values of all magnitudes entering into (6.32) are determined in (5.13). 
From formula (6.30), it can be seen directly that two inequalities may 
be realized: 
AE‘ (qj) > AEL(0) (a), (6.33) 
AEL(qj) < MELO) (0). 
In these, 


AES(O) = AEL, + D£,(0) + ef(0). 


The above expression gives the excitation energy for a crystal without 
lattice vibrational excitation. 

Which of the inequalities given by (6.33) is realized, a or b, depends 
on the zone number uy, on the sign of the matrix elements entering into 
(6.31), on their relative magnitude and direction, and on the magnitude 
of the wave vector q for the phonon. Thus 


AE{ (qj) — AB{(O)  hevgs (6.34) 


The result obtained is nontrivial. The frequency difference for two 
absorption lines (excitation with lattice vibration and without lattice 
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vibration) is not equal to the vibrational frequency. Moreover, the fre- 
quency of the light absorbed with simultaneous excitation of a phonon 
may in certain cases be less than the frequency of light absorption with- 
out phonon excitation. In other words, even at absolute zero, the fre- 
quencies of the light absorbed upon simultaneous excitation of an intra- 
molecular state and of lattice vibrations may be situated on the long-wave 
side of the frequency that corresponds only to an intramolecular exci- 
tation of the crystal. 

Furthermore, in contrast to the case in which excitation of the crystal 
occurs without excitation of phonons, and owing to the impulse conser- 
vation law [(5.27)] that the transition take place only at a determined 
point of each energy zone, with a simultaneous excitation of excitons and 
phonons, the impulse conservation law leads only to the following require- 
ment: The sum of impulses for the exciton hk and for the phonon hq must 
be equal to the light-quantum impulse hQ. Since q possesses all possible 
values (compatible with a cyclic condition), there exists a series of possi- 
bilities for light-quantum absorption, which are caused by a diverse distri- 
bution of the light-quantum impulse between the impulse of the exciton 
and that of the phonon. Owing to this, there appear, not absorption 
lines, but absorption bands, even when thermal motion is absent. 

These absorption bands must be continuous. However, since the 
matrix elements Y%,(qi, 7) entering into the transition matrix element 
[(6.29)] depend essentially on the wave number q, and since the density 
of states for the various values of q are not identical, the continuous 
background will have a complicated distribution of intensity. This effect 
is apparently the main factor determining the width of molecular mono- 
crystal absorption bands at a temperature near absolute zero. 

To determine the form of the absorption band (the distribution of 
absorption intensity inside the band), it is necessary to calculate the 
probability of a transition from the initial state (with energy E°) to a 
state with a given energy F corresponding to light absorption of fre- 
quency w; the latter is determined from the condition 


E = E° + ho. (6.35) 
The probability of a transition in the unit of time P to a final state 
with given energy will be determined by the transition matrix element 


[(6.29)], upon fulfillment of the selection rules [(6.27)] with the help of the 
well-known formula (Ref. 68, p. 269): 


-= y i [(uf|w|0)3,%0(E) do. (6.36) 


The differential d@2 = sin 3 dd dg is the element of the solid angle charac- 
terizing the direction of the phonon’s wave vector; the integration is 
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carried out over all possible directions q; the summation is over various 
polarizations of the phonon: 


p(B) = hag? $4. (6.37) 
(a9) 


The derivative dq/dw and q are calculated as functions of w, %, and gfrom 
the equation 
ho = ABLy + Dhg(0) + e{(—a) + hea (6.38) 


where ¢{(—q) is determined by (6.31). 

With the aid of (6.36), we can calculate the probability for the absorp- 
tion per volume unit (absorption coefficient) of one light quantum per 
second. 

The number of quanta of light of frequency w impinging in one second 
on 1 cm? of crystal surface, perpendicular to the direction of light propa- 
gation, is equal to E?c/4mrhw, where c is the velocity of light. Therefore, 
the absorption coefficient upon excitation of a crystal to the pth zone of 
the fth excited state will be equal to 


Tosa) = Hage fost |0), PoC) a 


Here V is the volume of the crystal absorbing the light. Substituting 
the value for the transition matrix element from (6.29) and (6.37), 
we obtain 


Lapa) = SON ahe (, ye i fe poeta) f tvicanlat $8 an 
(6.39) 


Here N is the number of molecules in a molecular crystal of volume V, 
o is the number of molecules per unit cell of the crystal, w is the frequency 
of the light absorbed, e is the electronic charge, c is the velocity of light, 
Ep is the electric field amplitude of the light wave, Y{k(q7) is determined 
by (6.25), and q is the wave number of the phonon. 

The wave number q for each polarization 7 of the phonon is determined 
from the frequency of absorbed light with the aid of the relation given by 
(6.38). 

Thus, to calculate the absorption coefficient [(6.39)], it is necessary to 
know explicitly the dependence of frequency w,; on the wave number q 
and on polarization 7 of the phonon. One must know also the depend- 
ence of the wave vector q on the frequency of absorbed light w for each 
direction determined by the solid angle d2. The dependence of the first 
type we; = w,;(q7) is determined by the lattice dynamics. Unfortunately, 
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for the majority of crystals, this dependence is unknown. Calculations 
have been made only for a small number of very simple lattices. How- 
ever, the dependence of q on w and dQ is determined by (6.38), and in 
principle it may be calculated if the dispositions of molecules in a crystal, 
as well as the wave functions for the molecular states (for free molecules), 
are known. 

Although it is at present practically impossible to use formula (6.39) 
to determine the distribution of intensity in the band, it is possible to 
make a few general qualitative observations. 

The statistical weight of states (for small q’s) is proportional to q’. 
Furthermore, the matrix elements Y/, depend on q in an obvious and a 
nonobvious way. At the same time, the matrix element increases with 


Fie. 11. Schematic representation of absorption of light in the region of an absorption 
band of a crystal. 


the increase of q, since for small q the phase difference for vibrations of 
neighboring molecules is small. Therefore, the derivatives dK/,/dR2*, 
entering into Y/, through W/,(qjn) [cf. (6.5)] and allowing for the changes 
of matrix elements Df, and M‘/,., upon displacement of molecules from 
equilibrium positions, will also be small. Thus it is to be expected that 
the absorption coefficients will attain generally different maximum values 
for different directions and polarizations of the absorbed light. 

Figure 11 represents schematically the possible dependence of the 
absorption coefficient on the frequency of impinging light. This figure 
enables us to draw the following conclusions, which are apparently con- 
firmed by the experiments of A. F. Prikhotko.1432_ In the latter experi- 
ment the absorption of light by naphthalene monocrystals of thickness 
varying from 1.5 to 0.05 mm was investigated, at a temperature of 20.4°K. 
Thick monocrystals absorb in a broad spectral range since even a small 
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absorption coefficient (calculated per unit volume and corresponding to 
minima of the absorption) is sufficient to totally absorb light. The 
dashed lines a, b, and ¢ in Figure 11 indicate the absorption coefficients 
sufficient for total light absorption for the corresponding thickness of a 
crystal. Line @ corresponds to the thin crystal, line ¢ to the thickest one. 
As the thickness of a crystal decreases, a total absorption will be observed 
only for frequencies corresponding to the intersection of straight line b, 
and subsequently of a, with the curve of absorption. Thus, for very thin 
laminae, the absorption spectrum will represent a series of sharply defined 
isolated absorption bands. 

Figure 12 shows (from A. F. Prikhotko) the portion of an absorption 
spectrum for naphthalene monocrystal at 20.4°K for two crystalline 
thicknesses (the direction of the electric vector for the impinging light is 
parallel to the crystal b axis). 


thick 


31620 31955 32255| 32673 32948 cm"! 
32411 


Fie. 12. Absorption spectra of monocrystals of naphthalene at 20.4°K. (After 
Prikhotko.13.32) 


The distances of the subsidiary absorption bands from the basic absorp- 
tion band depend on e{(—q) and will, generally speaking, be different for 
the various zones uw and for the various molecular excitations f. In other 
words, each basic spectral absorption band, corresponding to a given 
electronic-vibrational excitation of the molecules of the crystal, will be 
surrounded by its own subsidiary absorption bands. The latter are 
caused by the interaction of excitons with lattice vibrations. Con- 
versely, in the fluorescence spectra of crystals at low temperature, one 
can directly observe lattice vibration frequencies, which will appear in 
the form of satellites near the basic lines. 


27. Change in the Structure of an Absorption Spectrum with 
Temperature Rise. Influence of Translational Symmetry 
Violations in a Crystal 


In the preceding section we studied the structure of a band of light 
absorption by a molecular crystal at low temperatures (7 ~ 0). Let us 
now consider how the appearance of the spectrum will vary when the 
temperature rises. 

In this case we must first of all consider the fact that the average value 
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of the numbers of phonons for each type of lattice vibration wave will not 
be equal to zero and will increase with the temperature: 


(Noses = |exp () ~ ie 


This phenomenon leads to an increase in intensity for the absorption 
bands caused by the production of phonons and was investigated in the 
preceding section. The probability of the transition [(6.36)] found for 
the case T = 0 must be multiplied by the factor (Nq;)are + 1 to obtain 
the probability of transition for T #0. Moreover, we must keep in 
mind that the magnitudes entering into the matrix element (uf|w|0)/, also 
depend in a nonobvious manner on the temperature, since the Y/,(qj) 
depend on molecular displacements and since the magnitude of the latter 
varies with temperature rise. 

Furthermore, for T ~ 0, the matrix elements for the transitions dur- 
ing which the phonons are absorbed, together with the matrix element 
responsible for the production of phonons, will be different from zero. 
The probabilities for such transitions depend on the factor (Naar, and 
on the matrix elements, whose magnitudes also increase with temper- 
ature rise. 

Thus the increase in temperature leads to a general increase in the 
intensity of absorption bands, which are caused by molecular excitation 
in a crystal with simultaneous excitation of lattice vibrations. It further 
leads to the appearance of supplementary absorption bands, which are 
generated by the excitation of a crystal already having excited lattice 
vibrational states. The intensity of these supplementary bands will 
increase with temperature rise. 

Upon further temperature rise, the molecular displacements become so 
great that one must, in the decompositions [(5.1)], take into consider- 
ation the terms proportional to the quadratic and higher powers in 
molecular displacements from equilibrium positions. This leads to the 
production and absorption simultaneously of two or more phonons. At 
the same time, the displacement of basic absorption bands may also be 
observed, for the reason mentioned in Section 25. However, when the 
processes of the second and higher orders become sufficiently probable, 
the decomposition into powers of displacements itself loses meaning; it 
therefore will be impossible to apply perturbation theory. 

All results obtained above pertain to an ideal periodic lattice. How- 
ever, the ideal periodicity of a lattice is violated for a number of reasons: 
presence of admixtures, presence of isotopes, and finally, presence of 
thermal density fluctuations. Furthermore, it is also necessary to con- 
sider the finiteness of the crystal dimensions, i.e., the influence of the 
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surface.t If monocrystals possess sufficiently large dimensions com- 
pared with the lattice constant, the influence of the surface will not be 
significant. However, the violation of periodicity may lead to significant 
changes in the absorption spectrum. 

Conditions (5.27) and (6.27) were a result of the crystal translational 
symmetry. These conditions expressed the fact that, in a crystal under 
the influence of a light wave, only excitons with a wave vector equal to 
the wave vector of the light wave are excited if the lattice vibrations do 
not participate. If the lattice vibrations do participate, those excitons 
are generated that correspond to the differences in the wave vector of 
the light wave and in the wave vector of the lattice vibration. If the 
periodicity of the lattice is violated, a change in conditions (5.27) and 
(6.27) occurs. The excitation of excitons with other wave-number values 
as well becomes possible. We must expect that, as a result, wide absorp- 
tion bands appear in a quasi-stationary model. If a crystal does not 
contain various isotopes or admixtures, the periodicity will even then be 
violated by fluctuations. This leads us to believe that the violation of 
periodicity in a crystal with temperature rise is still another reason (in 
addition to those mentioned in the preceding section) for the widening of 
absorption bands. In particular, this effect of absorption-band widening 
must play an important role in the absorption of light by liquids. 

Here we shall limit ourselves merely to these qualitative considerations. 
A quantitative accounting for the influence of factors on absorption spec- 
tra causing the violation of crystal translational symmetry must be the 
subject of a special study. Let us merely observe that a series of very 
interesting works by I. M. Lifschitz7® has already been devoted to the 
problem of the optical behavior of nonideal crystal lattices in the infrared 
region. 

If the molecules possess degenerate energy states, and if the degener- 
acy is not removed by the field of the crystal lattice (for T = 0), then for 
certain molecular symmetry classes a splitting of degenerate states (upon 
temperature rise) may occur. This results from the change of symmetry 
of a crystal unit cell caused by certain types of lattice vibrations. Appar- 
ently the effect of the splitting of absorption bands will be significant 
only for comparatively high temperatures in a crystal. 

+ The effect of the influence of a molecular crystal surface on the excitation of intra- 


molecular vibrations was investigated in the work of this writer and of V. M. Agrano- 
vich of Kiev State University.* 


CHAPTER VII 


“LOCALIZED” EXCITONS IN A MOLECULAR LATTICE 


28. Excitation Energy of a Crystal in the Case of “‘Localized”’ 
Exciton Formation 


As shown in Chapter IV, upon excitation of a molecule in a crystal 
there arises a local violation of the equilibrium molecular arrangement, 
since the excited and the normal molecules act with different forces upon 
their neighbors. If the excitation remains localized on one molecule, 
the change in interaction forces between the excited molecule and its 
neighbors will result in translational and rotational vibrations of mole- 
cules relative to a new equilibrium position. If, however, the excitation 
transfers from one molecule to another with a period shorter than the 
time necessary for a noticeable displacement of molecules to new equi- 
librium positions, then the local deformation of the crystal will not have 
time to take place. 

If the excitation transfer time is considerably greater than the time of 
molecular displacement, i.e., the inequality ra < 7 is fulfilled, then the 
movement of excitons will be accompanied by a local lattice deformation. 
As a result, the velocity of exciton motion decreases still more, and its 
effective mass increases. The excitons whose motion is accompanied by 
the motion of a local crystal deformation we have called “localized”’ 
excitons. In this chapter, we shall study the excitation energy of a 
crystal upon formation of “localized”’ excitons. When calculating the 
excitation energy, one must keep in mind that, upon fulfillment of the 
inequality given by (4.22), the adiabatic approximation is fulfilled only 
for a nonexcited state (T ~ 0). Thus the normal state of a crystal will 
be described by the function 


V(r, R) = Aod?. (7.1) 


Here Ao satisfies equation (4.23) and is given by (4.26), and ©} coincides 
with (4.8). 

For excited states of a crystal, it is impossible to discard the operator 
[(4.14)] when studying the solutions of equation (4.5). Generally a rigor- 
ous solution of this equation is hardly possible for this case. However, 
approximate solutions of this equation can be obtained. 

97 
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Owing to nonfulfillment of the adiabatic conditions, it is not possible to 
seek solutions for equation (4.5) in the form of stationary states for which 
one or several molecules in the crystal are excited, the remaining mole- 
cules not excited, and all molecules in some state of motion. Such a state 
will not be stationary in the sense that, in the system, transitions will 
take place (with conservation of the total energy) from excited molecules 
to lattice vibrations (for T ~ 0) until thermodynamic equilibrium is 
established over the degrees of freedom for the system. 

In classical language, when the inequality given by (4.22) is satisfied 
for the molecular excitation, the system must be considered as one with 
great anharmonicity in the region of the excited molecule. Because of 
this anharmonicity, an energy transfer takes place from some normal 
vibrations to others. We may also formally describe this process of 
energy transfer in the language of a harmonic model, considering that 
the vibrations, corresponding to the molecular excitation energy, occur 
with positive damping and that the lattice vibrations occur with nega- 
tive damping. Considerable damping of the vibrations of the oscillator, 
corresponding to the excited molecule, will result in a wide absorption 
band. 

This classical analogy may facilitate finding an approximate solution 
of equation (4.5). When investigating the problem of light absorption, 
we shall initially look for quasi-stationary states (stationary during a 
short time) of equation (4.5); i.e., we shall disregard the terms causing 
transformation of excitation energy to lattice vibrations. Then we shall 
consider the influence of these terms as the reason for the widening of 
excited levels, owing to the short lifetime of such states, and as a mecha- 
nism for transforming absorbed energy into heat. 

Inasmuch as a state with a given total energy U of a system may be 
produced by numerous means according to the distribution of energy 
among the various degrees of freedom, we shall seek the solution of equa- 
tion (4.5), corresponding to energy U, in the form 


Wir, R,t) = [ a(t) AG(R) #5 (r, Rk) + o b(nat) A{(na)xh, + °° | eae 
(7.2) 


where the first term represents the states of a system without excited 
molecules (the excess of energy over zero-point energy for the system is 
distributed only among lattice vibrations), and the second sum corre- 
sponds to the case in which there is one excited molecule in a system 
and the remaining ones are in a normal state. Subsequent terms are of 
little importance if energy Uis small. 4, is the wave function for inter- 
nal molecular states of a crystal, for which one molecule is excited on 
site na; A{(na) determines the vibrational state of molecules in a lattice 
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with local deformation in the region of the lattice point na. By virtue 
of the equivalence of molecules in a crystal, this state is degenerate. 
Therefore, we shall seek the solution in the form of a sum of solutions 
corresponding to the localization of excitation in various molecules of the 
crystal. Coefficients a and b(n), the squares of whose moduli are pro- 
portional to the probabilities of the corresponding states for the system, 
will, generally speaking, be slowly changing functions of time. 

Let us now substitute (7.2) in (4.5), multiply by x4%, and integrate 
over all internal coordinates of molecules dr. By virtue of x%, being 
orthogonal to $$ and to all x4,, if m8 = na, we shall obtain 


b(ms, t)[8! — VI(R) — Trl A{(ms, t) = > b(na, t) A{(na) M4 


naxtmp 


+ » b(na, t)rha ma + ar? — Af(ms) 


no ,mBp 


obmnB, 2) (7.3) 


where =U — » oa — ALL». (7.4) 
ma 
The energy of molecular vibrations in a lattice is &’. The term 1’(R) is 
given by 
Vi(R) = V(R) + Di(R). (7.5) 
Here V(R), Dig(R), and M‘, .,g have the same values as in (4.24), (3.13), 
and (3.10), respectively.t Also, 


1 aA{(n axhg 
Tha,mp = Al(na) | xat0d. dr — MD ar M; oatee) ®) xe 5 ee dr, (7.7) 


° ° * ° 1 dAG 4 OP° 
T AG / xe! RPS dr — h? >, Mae, xa aR, dr. (7.8) 


We shall consider as an initial state, after light absorption by a crystal, 
that for which |a|/? = 0 and |b(ma)|? = 1. Furthermore, we shall con- 
sider the solutions of (7.3) only for time intervals small in comparison 
with the lifetime of this initial state. Then a = 0, [db(mf, t)]/at = 0, 


t Here we assume that functions ¢° and y/, may also be written in the form of the 
product of corresponding wave functions for separate molecules [see (3.2) and (3.5)]. 
However, the qualitative conclusions obtained by us do not depend on this assumption. 

If we eliminate this assumption, then instead of (7.5), we must use 


V1(R) =i4 > i eV nana’ Xap dr (7.6a) 


nana’ 


and Ma mB = | Vg VineumBXna dr. (7.6b) 
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and one can disregard the last two terms in (7.3). This consideration 
does not permit us to take into account the width of the absorption line, 
because the latter is precisely dependent on the probability decrease of a 
state for which 

[b(ne)|? = 1 and jal? = 0. 


We shall subsequently modify the solution of our problem to be able to 
consider also a time which is comparable with the time of transition to a 
new state (i.e., with a magnitude the inverse of the transition proba- 
bility). This allows us to obtain a formula for the width of the absorp- 
tion line. 

Thus, after limiting consideration to only small time intervals and 
introducing the definitions 


Min = Mhamg + Gants (7.90) 
f 
V"(R) = VCR) + qaenh, (7.9b) 


we may reduce equation (7.3) to the form 


b(m)(e — Vir(R) — Tr}Af(mB) = J. d(na)Af(na)bamg (7.9) 
nase 
Here A{(m@) determines the state of lattice vibrations prevailing when 
one molecule is excited at lattice point m6 and the remaining ones are in 
a normal state. 

Let us suppose that the localization of the excitation takes place in 
molecule m§; then all b(na) will be small at the first moment, if na ~ m8. 
If, besides, the matrix elements M‘, ,, are sufficiently small, then, dis- 
carding in the zeroth approximation all their products by b(m«), we obtain 
an equation determining the vibrational state of a crystal when the exci- 
tation is localized in molecule m@: 


[Tr + V"(R))At(mB) = 8, Ai(ms). (7.10) 


In this, 8&2 is the lattice vibrational energy for localization of the exci- 
tation in lattice point m8. 

Equation (7.10) is analogous to the equation for translational and 
rotational molecular vibrations in a nonexcited crystal [cf. (4.20)], but 
with a different potential energy: 

1 


Tne mp . 
Al(ms) (7.11) 


Vir(R) = V(R) + DER) + 


The minimum for this energy will correspond to the new disposition of 
molecules in the vicinity of the excited molecule, in comparison with the 
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case of a lattice in which the excitation is not localized or where it is 
absent altogether. 

We shall henceforth designate the deviation of molecules from these 
new equilibrium positions by means of R’. It is obvious, then, that 
Rf = 0 may be determined from the variational condition 


\8V/7(R)] aro = 0. 


Let us decompose the potential energy [(7.11)] into a power series of dis- 
placements relative to the new equilibrium positions R’ = 0 and limit 
ourselves to quadratic terms of the decomposition. We can thereby 
obtain eigenfrequencies in a crystal in the harmonic approximation with, 
generally speaking, new normal frequencies in comparison with the nor- 
mal frequencies of the nonexcited crystal. At the same time, the wave 
functions for molecular vibrations in a lattice, upon localization of exci- 
tation in the region of molecule m8, may be written in the form 


Ad(mg) = II Vi,,,(m8). (7.12) 


qi 


Here W¥,,(m8) represents the wave functions for harmonic oscillators, 
corresponding to normal lattice vibrations relative to new equilibrium 
positions, which arise upon molecular excitation in the mth lattice point. 

The transition from a system of connected homogeneous equations 
(7.9) to a system of independent equations (7.10) determining the 
degenerate (cN-multiple) system of functions A{(m@), which belong to 
one energy &°, is justified if the following inequality is fulfilled: 


Mime R’) 
lec <1, (7.18) 


As a rule, the inequality given by (7.13) is fulfilled in (chemical) solu- 
tions where, owing to large distances between absorbing molecules sepa- 
rated by a solid or liquid solvent, the excitation energy transfer from 
one absorbing molecule to another is rendered difficult. However, this 
inequality may also be fulfilled in molecular crystals for certain excited 
molecular states. Generally, however, it is not fulfilled, and it is inad- 
missible to discard the right-hand portions of equations (7.9). 

Owing to the presence of the right-hand portion in equations (7.9), the 
oN-multiple degeneracy of energy &/° is removed. 

The system of equations (7.9) is considerably simplified if we take into 
account the translational symmetry of a crystal and select for the coeffi- 
cients b(ma) in (7.3) [determining (in a quasi-stationary state) the wave 
function for the excited state of the crystal 


W(r, R) = ¥) b(na) Af(na)xh (7.14) 
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where only one molecule is excited in the fth state] the following 
expression: 
b(na) = Be, (7.15) 


Here the wave vector k has the same discrete values as in (3.14). 
Substituting (7.15) in (7.9), we obtain a system of o equations deter- 
mining the energy for system &/ as well as the unknown coefficients B,: 


Cc 


Bphe— ) Balfe(k) = 0, (7.16) 
a=l 
where 
Ae = &f — 8°, (7.17) 
Lf,(k) =D" [ A{(OB) MI, og Af(na)et™. (7.18) 


The matrix Leg(k) is Hermitian: 
Lae(k) = Lj. (k), 


and for crystals in which each molecule is the center of symmetry, this 
matrix is real and symmetrical relative to the indices a and 8. 

The setting equal to zero of the determinant, formed from coefiicients 
with unknowns, is the condition for the solvability of the system of 
homogeneous equations [(7.16)]._ This condition provides an equation of 
degree o for determining Ae (as functions of the wave number k). 


[Ldg(k) — Ae daa| = 0. (7.19) 
The roots of this equation determine o functions 
Aci = Aci (k), w=1,2,...,¢. (7.20) 


Using (7.17) we obtain 
&f = 6f° + Aef (kk). (7.21) 


Taking into consideration (4.23a¢) and (7.4), we obtain the excitation 
energy of the crystal in the following form: 


AFL (I) ing = AELy + 67° — 6° + Acf(k). (7.22) 


Thus, when “localized” excitons are formed, the crystal excitation 
energy will be determined by the excitation energy of one molecule AEZ,, 
by the change in lattice energy (&/° — &°) upon formation of localized 
deformation in a crystal in the vicinity of an excited molecule, and by the 
energy Ae{(k) connected with a movement of excitation together with the 
local crystal deformation and depending on the index yu, which runs 
through o« values, as well as on the wave number k of the excitation 
waves. Owing to the presence of the last term, the excited state of the 
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crystal corresponding to one molecular excited state will be characterized 
by » energy zones, just as in the case of free excitons. 
The separation between zones and their width will be determined by 
the matrix elements 
JAG(08) AL, op Af(na) dR. (7.23) 


For example, in the case of crystals containing two molecules per unit 
cell, 
Aé,2 = Ey, (0k) = Dy(k), (7.24) 
where 


Lui(k) = 5 et f A{(01) D159 Af(nl) d (7.25) 


Li(k) = [ Af{(1) Mf, A{(02) dR 
+ Yi" eft [ Af(O1) M5 9, Af(m2) dR. (7.26) 


Consequently, the separation between zones (for k = 0) is equal to 


21110) = 2) f A{(01) 1145 9 Af(n2) dR. (7.27) 


Thus, in the case of “localized” excitons, both the zones and the dis- 
tances between them will be considerably smaller than in the case of 
free excitons, where the corresponding magnitudes are determined only 
by the matrices MZ, .s. If we consider only the first terms of the decom- 
position [(5.1)] for matrices MZ, ng in powers of displacements, then 


SA{(08) ME, op Ai(na) dR ~ M2 o8(0)Sna.08, (7.28) 
where Snap = fA{(08) Af{(na) dR. (7.29) 


The integral Sroog represents an overlap integral for wave functions 
corresponding to crystal vibrational states for which there is a local defor- 
mation in the region of the Ofth lattice point and, correspondingly, in 
the region of the nath lattice point. 

In the zeroth approximation (when the mass of molecules in a crystal 
is infinitely large), functions A{(08) and A{(na) may be represented in 
the form of 6 functions: 


Aj(08) = I] 5[Roa(t) — Regs (2)], 
(7.30) 

Af(na) = I LF [Rre(t) can RA), 
where Rog(t) — Rfg(z) is the displacement of the ith normal coordinate 
from among the 6¢N normal coordinates, determining the position of 
molecules in a crystal away from the aquilibrun value Rfg(z), corre- 
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sponding to a local crystal deformation in the O@th lattice point; 
Rra(t) — Rf(z) is the displacement of the 7th normal coordinate away 
from the equilibrium value Rf8(z), corresponding to a local crystal defor- 
mation in the rath lattice point. 

If we utilize (7.30), then 


di if Py = mg 


Spyme = SPy,mep = | 0 if Py < ms. (7.31) 


If we take into account the finiteness of the molecular mass for a crystal, 
we obtain, in lieu of 6 functions, the harmonic oscillator functions: 


AS°(mg) = Nez | exp (- aa) | Hid), (7.32) 


where x2(2) = h/(2M.w;) is the average quadratic deflection for zero-point 
vibrations with frequency w;, & = [Rna(t) — R%G(2)]/[2x0(2)], and H,(£) is 
the Chebyshev-Hermite polynomial of the nth degree. 

In the zeroth approximation, by virtue of (7.28) and (7.31), we obtain 


SAL(08) Mt, op Al (na) dR = 0. (7.33) 


Consequently, the width of zones and the separations between them are 
equal to zero in this approximation.t In reality, of course, this is too 
coarse an approximation; however, it is to be supposed that both the 
splitting of zones and their spread have small value because of the small 
magnitude of matrix elements [(7.23)]. 

The physical meaning of the smallness of the integrals [(7.23)] can be 
described approximately as follows. For an excitation to pass from one 
molecule to another in the case of a “‘localized’’ exciton, it is necessary 
that the wave functions for vibrational states A(08) and A(na), corre- 
sponding to the localization of excitation in two molecules, overlap each 
other, i.e., that Sogne ~ 0. If Sosne = 0, the excitation, together with 
the local deformation, cannot pass and remains frozen” in one molecule. 

To each of the » bands of the excited states [(7.22)] there will corre- 
spond its own wave function 


W4,.(, 1, RB) = wa » eit » Af (na) BxL,, (7.34) 


where coefficients B# are the solutions of homogeneous equation (7.16) 
for the uth root Ae of equation (7.19). The wave function [(7.34)] is 
normalized in such a manner that [Bz]? = 1. 

{In this approximation the group velocity for excitation waves, corresponding to 


localized excitons, is also equal to zero, but the effective mass of excitons is equal to 
infinity. 
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The functions found by us describe the behavior of a crystal only dur- 
ing a very short time, since the corresponding states are not stationary 
states. 

In order to calculate the width of the absorption band which arises as 
a result of the corresponding excited-state levels being nonstationary, 
we can use a method applied by Weisskopf and Wigner’! for calculating 
the natural width of a spectral line. 

Inasmuch as we shall be concerned with the change of crystal energy 
levels in time, we must consider the solution of the time-dependent 
equation (4.5). 

According to the preceding investigation for time ft ~ 0, we may con- 
clude that the velocity of localized exciton travel is very small; therefore, 
phenomena connected with the propagation of excitation over the crystal 
(splitting and polarization) do not play any noticeable role. Conse- 
quently, in order to study the problem of transformation of molecular 
excitation energy into heat energy (lattice vibrational energy), it suffices 
to consider, as a possible solution of equation (4.5), the wave function 


U(r, Rt) = a(t) A5(R)%(r, Re + B(na)xhq(r, R) AX(R, nae, 
(7.35) 


where A{(nqa) satisfies equation (7.10) and characterizes lattice vibrations 
for localization of intramolecular excitation on molecule na; hw» is the 
total energy of this crystal state, which we shall here consider an initial 
one; A§(R) satisfies equation (4, 23) and determines lattice vibrations 
when there is not a single excited molecule in a crystal; and hw is the 
total energy for the corresponding state of the crystal. 

Let us substitute (7.35) in (4.5), then multiply the expression obtained 
by ®}*A6(R) and by x{X Af(na); let us now integrate over variables r and 
R. We then obtain a system of equations: 


th a) = b(t, na) MPe ort 
t (7.36) 
ay db(t, os ne) a(t) MP eo—eot 
where 
Mf? = » / Acta) a a dR +f ov’ Al (na) dR, 
(7.37) 
off) ° OxX,a 
rv?(4) = en 


h? Qf 
pee” / a > un ye ar, (7.374) 
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Let us suppose that for ¢ = 0, b(0, ne) = 1 and a(0) = 0. Then we 
can determine the change of a(t) by taking into account for b(¢, na) the 
value which satisfies the initial conditions 


b(t, na) = eH", (7.38) 


where 7 is the average lifetime of the initial state. 
Substituting (7.38) in the first equation of (7.36), we obtain 


MSE (elieo-e)-VANIt — 1) 


ho wo — w +1/(27) (7.39) 


a(t) = 


The value found for a(t) is substituted in the second equation of (7.36); 
then 


2 J — elie—eo)+l/ (re 2 
= th |AA7?|? a | 


woo + i/) M/*|? 8(w0 — a). (7-40) 


In this manner, the transition is accomplished without energy change 
(w = wo). There exist a very large number of states, determining lattice 
vibrations, with the very same energy hw; these states are distinguished 
by a different number of phonons and by the direction of the propaga- 
tion g and polarization j of these phonons. Furthermore, we are not 
concerned with the transition to one of these states, but to any of them. 
In view of all this, the total value of 7 may be obtained by integrating 
(7.40) over all these states in the range of values for Aw. If we designate 
the density of phonon states with wave vector q and with polarization j 
in the frequency range dw through p,;(w), then 


1 2 
22" fim pu(oo) do. (7.41) 
i 


Thus, 1/r is equal to the total probability for the transformation per 
unit time of intramolecular vibration energy into heat. If we designate 
this probability by Pr, we obtain 

Pr= (7.42) 


The width of level, which is determined by the finiteness of the lifetime 
of the state, will be equal, in energy units, to 


AE ~" = APr, (7.43) 


where Pr must be taken from (7.42) and (7.41). 
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29. Comparison of Crystal Excitation Energy upon Formation of a 
‘“Localized’’ and Free Exciton 


Let us now compare the crystal excitation energy upon formation of 
a “localized” exciton 


AEL(K)ice = AES, + &/° — & + Act. (k) ~ AES, + &° — 8° (7.22) 
with the excitation energy upon formation of a free exciton 
AES(k) = ABEL, + Dig(0) + ef(k). (5.15) 


The fundamental distinction between excitation energy values amounts 
to the dependence of excitation energy 
on the wave number of the exciton. 
By virtue of (7.28), the following in- 
equality is always fulfilled: 


O =~ [Aci ie(k)| « lef(k)|. (7.44) 


Furthermore, D/,(0) in (5.15) deter- 
mines the energy difference for the 
interaction between either an excited 
or a normal molecule and its neighbors 
[usually DZ,(0) <0]. This energy is 
calculated relative to equilibrium dis- 
positions of molecules in a nonexcited 
crystal. Now 6&/° — &°determinesthe Fic. 13. Dependence of excitation en- 
same diference for anew equilibrium GY, of et on eee ores 
configuration, plus the lattice vibration —ocalized exciton (I). 

energy change. In view of all this, 

and since the transition to a new equilibrium state cannot be accom- 
panied by an energy increase, 


DS Sher 8°, (7.45) 


Excitation energy 


k 


Let us now determine the energy difference between a free and a 
“localized” exciton: 


AES (k) — AE{(E)ic = Dig(0) — (8° — 8°) + e(k). (7.46) 


The difference [(7.46)] may be either positive or negative, depending 
on the type of crystal, on the excited molecular states f, on the number of 
u zones for the excited states, and on the wave vectors of free excitons. 

Case A. The excitation energy of the free exciton is larger than the 
excitation energy of the “‘localized”’ exciton. This case is schematically 
represented in Figure 13. Section I of the curve pertains to the case of 
‘localized’ excitons; section II of the curve refers to the case of free 


108 DAVYDOV: THEORY OF MOLECULAR EXCITONS 


excitons. Owing to (7.45), in order to realize case A, it is sufficient that 
the following inequality be fulfilled: 


ef(k) > 0. (7.47) 


Thus, if (7.47) is fulfilled, the excitation of free excitons to the pth zone 
of excited states with wave number k will require more energy than the 
excitation of ‘ocalized”’ excitons (small k). 

Case B. The excitation energy of the free exciton is smaller than the 
excitation energy of the “localized” exciton. In order to realize this case, 
it is necessary that 


Di gO) — (8° = 8°) < — i(k), (7.48) 
ie., that d(k) <0 


and that ¢/(k) be sufficiently large in absolute magnitude. Figure 14 

shows schematically the relative disposition of excitation levels for both 

types of excitons (I, “localized” excitons; II, free excitons) for this case. 

We must note that, in all works 

known to us?}7273 in which exciton 

I concepts were utilized, the following 

\ was always assumed: Free exciton levels 

are disposed higher than the level of a 

\ “localized” exciton. In this manner, 

case B, as indicated here, was not an- 

Oe alyzed. Yet, as we shall see later, 

this case may be realized compara- 

tively often and leads to interesting 
phenomena. 


Excitation energy 
- 
= 


k * 30. Excitation of ‘Localized’ Excitons 


Fie. 14. Dependence of excitation by Light. Structure of the 
energy of excitons on wave number. Absorption Band 


Case B: Localized exciton (I) great ; par 
than free exciton (II). ee In order to determine the probabilities 


for the transition from a normal state 
W° [(7.1)] to an excited state W{4,, [(7.34)] under the influence of a 
light wave [(5.19)], one must consider the matrix element 


(uf |w|O)roo = CEof Wfe* eit (r) h° dr dR. (7.49) 
Substituting (7.1) and (7.34) in (7.49), we obtain 


(uf|t|0) roe = AT ) ei(Q—-k)n 


x [> Bs | [] wine) - T] s. aR | eft (eee ar] (7.50) 
a=] a a 
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For the excitation of “localized” excitons, the correct selection rules are 
again 
oskee|® (7.51) 
+2ra}, 1 = 1, 2,3. 


Upon their fulfillment, the matrix element [(7.50)] is equal to 


(ufhol0n. = af >) Be | I] (filma yh, ae | olt(eed dr 


(7.52) 


Since wave functions V4,,(na) and WR, belong to diverse potential 
fields with diverse equilibrium positions, the integral 


} [] eu (va)¥hu] dR, (7.58) 


entering as a factor into ae, will be different from zero if N4; = Nq- 
In other words, together with the excitation of molecules, several phonons 
will be simultaneously excited (or absorbed at high temperatures). It 
must be emphasized that this result is obtained in the case in which 
“localized’”’ excitons are already formed in the zeroth approximation, 
whereas in the case of free excitons, the possibility of multiquantum 
phonon absorptions and emissions arose only in the highest approxi- 
mations of the theory. 

Thus, upon “localized” exciton excitation, wide absorption bands are 
to be expected owing to the combination of molecular excitation with 
lattice vibrations and owing to the fact that excited states are non- 
stationary. Because of the closeness of excitation energy zones, corre- 
sponding to various values of u, and because of the above-mentioned 
broadening of absorption bands, separate bands will not be permitted 
but will occur in the form of one wide absorption band. The distribu- 
tion of intensity in the band may depend on the direction of light propa- 
gation and polarization (there will appear bands corresponding to various 
values of u), as well as on the mutual disposition of polydimensional 
potential-energy surfaces as functions of the disposition of molecules in 
a lattice. 

Since, upon local deformation of a crystal, a change in equilibrium 
position of molecules takes place, a strong probability arises for con- 
siderable changes in vibrational quantum numbers. At the same time, 
noticeable intensity fluctuations must be expected in the band, since the 
value of the integral [(7.53)] changes with the change of AN,;, not mono- 
tonically, but in a rather irregular manner. This is the result of the 
chance disposition of wave-function nodes and antinodes corresponding 
to the initial and final vibrational states of molecules in a crystal. 


CHAPTER VIII 


TRANSFORMATION OF ABSORBED ENERGY INTO 
HEAT AND RADIATION 


31. Introduction 


If a crystal has passed into an excited state, and if there is no change in 
the composition of the solid as a result of a photochemical reaction, then 
such a crystal state is thermodynamically not in equilibrium. The tran- 
sition of a crystal to a thermodynamically equilibrated state may take 
place in two ways: 

1. The molecular excitation energy is transformed into heat; i.e., a 
distribution of excitation energy takes place among degrees of freedom 
which are connected with the molecular vibrations in the lattice. 

2. The crystal excitation energy is re-radiated in the form of electro- 
magnetic waves; luminescence of the crystal takes place. Sometimes 
energy of the same wavelength is radiated (resonance radiation); how- 
ever, more often, energy of lower frequencies is radiated, and a portion is 
transformed into heat. 

In this work, we shall concern ourselves only with pure substances 
which do not contain admixtures of foreign materials. Luminescence of 
pure substances has been comparatively poorly studied. The opinion 
already has been expressed” that, in a pure state, only those substances 
can luminesce in whose unit cell there is a noncompensated electron layer 
screened from surrounding neighbors. To such substances belong the 
rare-earth-element salts, uranyl salts, and others. However, a series of 
works has shown that many crystals of inorganic and organic compounds, 
although not containing foreign admixtures, produce luminescence. 

8. I. Vavilov showed” that luminescence of solid nitrogen is a typical 
example of luminescence of a pure substance in a solid crystalline state. 
In this case, the diatomic nitrogen molecules, forming a molecular-type 
crystal in the solid state, constitute the emitter. 

Especially convincing in this respect are the works of I. V. Obreimov, 
A. F. Prikhotko, and their collaborators.518 In their work, they showed 
that there is a bright luminescence even in completely pure monocrystals 
of naphthalene, anthracene, phenanthrene, and other organic compounds. 

The problem of transformation of excitation energy into heat in a solid 
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initially excited as a result of photon absorption was first investigated 
theoretically by J. I. Frenkel!%-?! for atomic lattices of inert gases. By 
means of qualitative reasoning, Frenkel pointed out the possibility of 
transitions in a solid which were not accompanied by radiation and which 
consisted in the emission or absorption of any number of phonons, or in 
particular, of a number such as would correspond to the transformation 
of the total excitation energy into heat. According to Frenkel, the 
possibility of such transitions arises from the following: The correspond- 
ing wave functions for vibrational states of a normal and an excited 
crystal are unlike owing to the displacement of equilibrium distances 
between atoms in a crystal at the moment of light absorption. This is 
especially pronounced for localized (‘‘adherent” as J. I. Frenkel calls 
them) excitons whose emergence is accompanied by local lattice 
deformation.t J. I. Frenkel considered that this transformation of 
excitation energy into heat always takes place when light is absorbed 
by a solid. According to Peierls,?? in certain cases a crystal must emit 
the absorbed energy in the form of ‘‘modified resonance”’ radiation and 
not transform it into heat, i.e., behave like a “diffuser.” Frenkel 
criticized Peierls’ work and showed (Ref. 21, p. 6538) that in all proba- 
bility there are no ‘‘diffusers’” within the meaning of Peierls’ theory. 
Besides, Frenkel proceeded from the assumption that, at the moment of 
light absorption, free and ‘‘localized” excitons are formed. The former 
pass into a state of “adherent” (“localized”) excitons. ‘Localized’ 
excitons subsequently transform their energy into heat. The trans- 
formation process of a free exciton into a “localized” one is, in Frenkel’s 
opinion, a necessary intermediate stage of the transformation of excita- 
tion energy into heat. 

Apparently, for some molecular crystals the conclusions of Frenkel 
are essentially correct. However, one cannot unconditionally apply 
them to all molecular crystals. 

In this chapter, we shall study the problem of the transformation of 
absorbed energy into heat for the case of molecular crystals consisting of 
polyatomic molecules. In such crystals, the process of transformation of 
molecular excitation energy into heat may take place in two different 
ways: 


1. Conversion of excitation energy (including electronic excitation) 
to lattice vibration energy with the aid of interaction forces between 
molecules 


t The probability of excitation energy (of a “localized” exciton) transforming into 
heat was calculated incorrectly by Frenkel. When calculating this probability, 
Frenkel did not consider that the wave function of the electronic states of atoms 
changes when the interatomic distances change. In this work a more rigorous expres- 
sion [(7.4)] was found for the probability of this conversion (Section 28). 
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2. Intramolecular conversion of electronic excitation energy into 
vibrational energy of the atoms which make up the molecules, and a 
further conversion of vibrational energy of atoms in a molecule to 
vibrational energy of molecules in a lattice 


The first portion of method 2, i.e., the transformation of electronic 
excitation energy into atomic vibrational energy, was called intramolecular 
deactivation by A. N. Terenin (Ref. 76, p. 74). A. N. Terenin maintains 
that “the transformation of all electronic excitation energy into potential 
energy of atomic nuclear motion is characteristic for polyatomic mole- 
cules.” Here we shall not study the reason for this phenomenon. We 
shall merely analyze the influence of the interaction between molecules 
on the transfer of the vibrational energy of atoms in a molecule to lattice 
vibrations. We shall also study the direct transition of molecular elec- 
tronic excitation energy to lattice vibrational energy. 

Experiments have demonstrated that, at low temperatures, the vibra- 
tional portion of energy for molecular electronic vibrational excitation is 
extremely rapidly (within a period of less than 10-8 sec) transferred to 
lattice vibrations. Among the most convincing works of this kind are 
those of I. V. Obreimov, A. F. Prikhotko, and their collaborators! 17.18 
on the study of naphthalene, anthracene, and phenanthrene absorption 
spectra and luminescence, as well as on the comparison of these spectra. 
In these works, they showed that at low temperatures the absorption 
spectrum corresponds to the electronic transition plus the vibrational 
frequencies of the excited molecule; the luminescence spectrum cor- 
responds to transitions from the lowest vibrational state of an excited 
molecule. Thus the vibrational energy of nuclei can pass over to the 
lattice vibrations in a period of time less than the lifetime of a molecule 
in an excited state. These experiments indicate beyond any doubt that 
the vibrations of nuclei, in any case those of the totally symmetrical 
type, i.e., those which combine with electronic transitions, quickly lose 
their energy to lattice vibrations. 

The total intramolecular excitation energy may be represented as the 
sum of vibrational and electronic excitation energy; i.e., one may assume 
the absence of interaction between both types of molecular excited states. 
This assumption is approximately justified in molecules, since intra- 
molecular vibration frequencies are considerably smaller than fre- 
quencies corresponding to the energy differences of two electronic states. t 

In such an approximation, one may consider the problem of the trans- 
mission of intramolecular excitation energy to lattice vibrations sepa- 

} To study intramolecular deactivation processes it is, of course, impossible to 
limit oneself to this first approximation; it is necessary to take into account the possi- 


bility of energy transition between both types of excited states in a molecule, i.e., pass 
over to the next approximation. 
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rately for the electronic excitation of molecules and for the excitation 
of intramolecular vibrations in a molecule. 

Therefore, we shall henceforth assume that the excited state of a 
molecule is the state which corresponds either to a vibrational or to a 
purely electronic excitation. 

When we study the question of transfer of intramolecular excitation 
energy to lattice vibrations, we must distinguish two cases: (1) Upon 
excitation of a crystal, free excitons are formed. (2) Upon excitation of 
a crystal, ‘localized” excitons are formed. 


32. Behavior of Free Excitons after Their Formation 


As was noted above, according to J. I. Frenkel, free excitons must 
after their formation turn into “localized” excitons. In F. Seitz’s 
work,” these views found the following graphic interpretation. In 
optical transitions excitons appear with wave vectors equal to the wave 
vector of light. Their group velocity is generally small but may be 
sufficient for a free exciton to form. However, the group velocity of an 
exciton must decrease when the wave number decreases. Therefore, if 
at the moment of formation the exciton’s speed is sufficiently great, then 
within a certain time, as a result of giving its kinetic energy to the lattice 
vibrations, the exciton loses its velocity and becomes “‘localized.”’ 

In these observations regarding the necessary transformation of a free 
exciton into a “localized” one, the above-mentioned authors seemingly 
overlook the following circumstance: The exciton energy loss is not always 
linked with a decrease in the wave number and speed of the exciton but, 
on the contrary, may sometimes lead to an increase of the speed of the 
exciton (see Figure 14). 

A second important circumstance which was not taken into account 
consists in the following: The energy of the free exciton may be less than 
the energy of the “‘localized” exciton (see Section 29). In this case, the 
free exciton will not be able to pass spontaneously into the state of a 
“localized”’ one, unless it obtains the necessary energy from heat fluctua- 
tion. At low temperatures, this last process is not very probable. 

We must, of course, also keep in mind that a free exciton state is not 
rigorously stationary, since there is a probability (albeit very small) for 
the direct transformation of exciton energy into lattice vibrational energy. 
However, this process plays practically no role at all for free excitons, 
since usually the free exciton gives up its energy in the form of electro- 
magnetic radiation during a considerably shorter time. Thus the 
luminescence of states which correspond to free excitons may occur only 
if the excitation energy of the free exciton is less than that of the “‘local- 
ized” exciton, and if, when wave number k increases, the energy of the 
free exciton decreases (the exciton has a negative effective mass). Fur- 
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thermore, the probability for transformation to a “localized” exciton 
state (or exciton state of other value u) is less than the probability for 
photon radiation. 

When an exciton with negative “effective” mass forms, the transfer of 
energy to lattice vibrations will be accompanied by an increase in the 
wave number of the exciton. As a result of this, the transformation of 
excitation energy into radiation will be rendered difficult, since the wave 
number of the exciton will no longer correspond to the wave number of 
the light wave. The lifetime (with respect to radiation) of the exciton 
may increase. 

Let us now find the probability for the transformation of a crystal 
from an excited state having a free exciton to the same excited state with 
a “localized” exciton. The wave function for the excited crystal with 
one free exciton is given by (5.10): 


wi = Ae Aste) », eit o Bixle). (8.1) 


asl 


The wave function for the excited state of a crystal with one localized 
exciton is given by (7.34): 


UH, = a » ems » BeAt(nad te (8.2) 


n 


If both these functions represented stationary states, ie., if they were 
exact solutions of equation 


(Tre + Ho) = UY, (8.3) 
then there would be no spontaneous transitions between them. In 
reality, neither solution is stationary; therefore, there can take place 


between them spontaneous transitions whose probability in 1 sec will be 
determined by 


2Qr 
Pr= = 


It follows from (8.4) that the transition probability differs from zero if 
the law of conservation of energy for the system as a whole is satisfied: 


Vio = Utree- (8.5) 


2 
| UO(Pe + Ho) Wty de dR 8(Uwe — Unes). (8-4) 


If the condition given by (8.5) is fulfilled, the probability for transforma- 
tion in 1 sec from one excited state corresponding to excitation with a 
localized exciton into a state corresponding to excitation with a free 
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exciton (or the reverse), will be equal to 


9. 
P, = = a2 (8.6) 
where Q = [wee'* (Te + Ho) Vf dr dR. (8.7) 


The probability P; depends on the degree of overlap of wave functions .A° 
and A{(m) for vibrational lattice states without local crystal deforma- 
tion and with local deformation in the region of one of the molecules of 
the crystal. Furthermore, P; depends on the derivatives of wave func- 
tions with respect to molecular displacements from their equilibrium 
positions. 

Besides transformation of free exciton excitation energy into lattice 
vibrational energy or into “‘localized’’ exciton energy (and lattice vibra- 
tional energy), transitions of free excitons from one excited-state zone 
into another are also possible. The probability of this type of transition 
was determined in Chapter IV [expression (4.18)]. 


83. Behavior of ‘Localized’ Excitons after Their Formation 


As we noted previously, excitons with local lattice deformation may 
form directly at the moment of light absorption or may result from free 
excitons if the free exciton energy is equal to, or greater than, the “‘local- 
ized” exciton energy. The probability of this transition is determined by 
expression (8.6). 

The “localized” excitons formed do not correspond to stationary crystal 
states. Because of the violation of adiabatic conditions, there will take 
place spontaneous transitions into states, in which all of the excitation 
energy will be distributed over the lattice vibrations. The probability 
of this process was determined in Section 28. 


h ox OX,e dA! (na) 
a (- Be fare Ske dr) Aloe) oR oe 
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We must keep in mind that the processes of energy transfer from the 
excited molecule to lattice vibrations are determined not only by a 
change in wave functions for vibrational states of a crystal, by which 
multiquantum transitions become possible, but also by the fact that 
corresponding states of a crystal are not quasi-stationary [by the non- 
zero value of the expressions in parentheses within the integrals of (8.8)]. 
If the adiabatic conditions were fulfilled (if there were stationary states), 
the above-mentioned change in wave functions for molecular vibrations 
in a lattice would result only in the absorption of light quanta with more 


Paj(w) d2. (8.8) 
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energy, and the surplus of energy over that necessary to excite molecules 
would be transferred to lattice vibrations. The process of total trans- 
formation of light into heat must necessarily pass through nonstationary 
states of a crystal. 

Since the lattice vibrational frequency differs less from intramolecular 
vibrational frequencies than from electronic excitation frequencies, we 
must assume that the probability of intramolecular vibrational energy 
transformation into heat is greater than the corresponding probability 
for electronic excitation energy. This apparently explains the rapid 
transition (mentioned at the end of Section 31) of the vibrational portion 
of electronic-vibrational molecular excitation energy to lattice vibra- 
tional energy. 

Thus, upon formation of “localized” excitons, as a result of such states 
being nonstationary, an excitation energy transformation to heat, i.e., 
the disappearance of the “localized” exciton, will take place within a 
certain time, if, of course, during the localized exciton’s lifetime (with 
respect to the transformation into head), 


1 lie oe (8.9) 


the exciton is not subjected to other transformations. Such transforma~- 
tions consist of: 


1. Disappearance of the ‘“‘localized” exciton due to radiation 
(luminescence) 

2. Transformation of a “localized” exciton into a free one 

3. Disappearance of a “localized” exciton as a result of intramolecular 
deactivation 


Luminescence from “localized” states is possible if the lifetime of the 
“localized” state rr with respect to the transfer of its energy to lattice 
vibrations is greater than, or comparable to, the lifetime of the localized 
state with respect to luminescence 7,. The latter is determined by the 
relation 

fie (8.10) 
YBa] (ufl0]0) ioel® 
where (uf|wl0),.. must be taken from (7.45). 

It seems that, contrary to the generally accepted opinion, the trans- 
formation of a ‘‘localized” exciton into a free one is also possible in 
certain cases. At low temperatures, a necessary condition for this 
process is the requirement that the excitation energy of the “localized” 
exciton be greater than, or equal to, the excitation energy of the free 
exciton. In Section 29 we showed that this may take place (case B). 
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In this connection, if the lifetime of the “localized” exciton with respect 
to the transformation into a free exciton state (r, = 1/P,) does not exceed 
the lifetime with respect to the transformation into heat rr of the 
“localized” exciton energy and the lifetime with respect to luminescence 
tt, then such transitions will occur and may subsequently result in 
luminescence (see Section 34). 

In some cases, apparently, an essential role is also played by the intra- 
molecular deactivation of excitation energy. The probability of these 
processes may be determined in principle when studying luminescence 
spectra of greatly rarefied vapors. t 


34. Necessary Conditions for the Occurrence of Luminescence of 
Solid Bodies Not Containing Admixtures 


Let us analyze the conditions that could favor the occurrence of 
luminescence of a solid. As was noted above, both free and ‘“‘localized”’ 
excitons may arise as a result of light absorption in a crystal. The 
luminescence will arise from the radiative transitions between excited 
states and the ground state. However, these processes will compete 
with a series of other processes which lead to the disappearance of both 
types of excitons. 

Let us determine the ratio of the number of excitons emitting their 
energy in the form of radiation to the total number of excitons formed. 
This ratio is equal to 

1 Pit Pr+ Pot Pi 


where Py is the probability that in 1 sec the exciton emits its energy in 
the form of luminescence 
Pr is the probability that in 1 sec the exciton gives up its energy 
to (thermal) lattice vibrations 
Pp is the probability that in 1 sec there takes place an intra- 
molecular deactivation which results in the destruction of the 
exciton 
P, is the probability that in 1 sec a transition of excitons of a given 
type (free or “‘localized”) to another type takes place 
Among possible processes of exciton disappearance we are not con- 
sidering molecular ionization processes due to heat fluctuations, which 
transfer the electron over to a free zone. These processes are essential 
only at comparatively high temperatures and for sufficiently highly 
excited molecular states. 


(8.11) 


t In a molecular crystal, however, the role of intramolecular deactivation may be 
less than in an isolated molecule, since, as a result of the transmission of excitation 
along the crystal, the lifetime of the excited state for one isolated molecule will be 
shortened, although the lifetime for the crystal emission is greater than that for the 
free molecule. 
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For luminescence to be observable, it is necessary that 
Pr~ (Prt+ Po t+ Py. (8.12) 


In the case of localized exciton formation, usually Pr>> Px; thus 
such states, as a rule, do not luminesce. If in some cases such lumi- 
nescence takes place, it must decrease when the temperature rises, since 
Pr increases with temperature rise. Actually, Pr [(8.8)] depends on the 
magnitude of molecular displacements from equilibrium positions (see 
Section 18). As a result of heat fluctuations, fluctuation-type displace- 
ments will be superimposed upon the displacements brought about by 
the exciton localization. If we designate the fluctuation energy neces- 
sary to attain a sufficient magnitude of displacements (corresponding to 
the intersection of potential surfaces) by Fr, the probability of Pp will 
increase according to the exponential law 


Pp rw e-Erike, 


In the case of free exciton formation, as has been previously noted, Pr 
does not play an essential role. Therefore, if the intramolecular deacti- 
vation is small, and if the probability of P; is small, a noticeable lumi- 
nescence will be observed. 

If the excitation energy of free excitons is greater than that of ‘‘local- 
ized” ones, luminescence is hardly possible. If, however, the excitation 
energy of free excitons is less than that of “localized”’ excitons, then 
a luminescence with a large quantum yield will be observed at low 
temperatures. In this case, luminescence will be promoted by transitions 
from “localized” states to free exciton states. 

Luminescence must decrease even in this case when the temperature 
rises, Since the number of free excitons will decrease at the expense of 
thermal transfers to localized states, which will subsequently give off 
their energy to lattice vibrations. 

When an exciton moves through a crystalline lattice, a portion of 
exciton energy is spent on emitting phonons. At the same time, for 
excitons with negative effective mass, the exciton wave number increases 
(k—k’ > Q). The effect indicated may be evidenced experimentally 
in two relations: 

1. The lifetime of a crystal excited state may increase, since optical 
transitions are permitted only when the law of photon- and exciton- 
impulse conservation is observed. If k > Q, the transition to radiation 
will be possible only upon simultaneous emission of a photon and a 
phonon, whose composite quasi impulse is equal to the exciton quasi 
impulse, or as a result of the presence of inhomogeneities in the lattice 
(crystalline surface, fissures, isotopes, heat fluctuations, etc.). The 
violation of lattice periodicity makes optical transitions possible even 
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when k ~ Q; therefore, at temperatures different from zero, an increase 
in lifetime is scarcely possible. 

2. The position of the absorption band corresponding to the origin of a 
free exciton with wave number k = Q will not coincide exactly with the 
luminescence band corresponding to the disappearance of an exciton 
with wave number k # Q. 

Thus we arrive at the following conclusion: If a solid is composed of 
molecules which, in a high-f excited state, luminesce (small intra- 
molecular deactivation), one will observe the luminescence of such a solid 
at low temperatures, provided that the magnitude Pr is small for 
‘localized’ excitons, and provided that, free exciton excitation with 
an excitation energy smaller than the localized exciton excitation energy 
is possible in the solid. 

A weak dependence of wave functions, which determine the optical 
behavior of a molecule, upon the distances between the molecules of a 
erystal and the deformation of a molecule is conducive to small values 
of Pr. Such weak dependence occurs in a series of cases: 

1. Solids containing atoms of rare-earth elements, uranyl salts, and 
others. The optical behavior of such solids is determined by internal 
electron transitions within 4f (or 5f) levels in atoms. The bonding of 
these atoms with other atoms depends on external electrons. In other 
words, in this case optical electrons differ from electrons which participate 
in interaction forces between atoms. As aresult, when distances between 
atoms and ions change, only the surfaces of external electron distributions 
are deformed; electron states of (internal) electrons of the 4f shell are 
little affected. Conversely, a change in the state of internal electrons of 
the 4f shell tells but little concerning the change of interaction forces 
between atoms in a crystal lattice. In this manner, the stationary 
character of states increases, since the derivatives 06°/dR; and 0°6°/aR? 
will be small, and since, consequently, the probability Pr will also be 
small. Absorption spectra of such crystals must not be strongly dis- 
placed in relation to the spectra of corresponding gases, inasmuch as the 
displacement magnitude of the absorption band is determined by the 
change in interaction forces between molecules upon excitation of one 
of the molecules. Crystals of this type should display intense lumines- 
cence (especially at low temperatures). 

2. Solids formed from molecules of organic compounds which contain 
multiple (double, triple) bonds. To these solids belong molecular 
crystals of the aromatic series. In molecules of these crystals, the r 
electrons of double (triple) bonds are involved in the optical properties, 
while the o electrons play the essential role of the formation of intra- 
molecular bonds. In addition, molecules of such compounds usually 
contain a large number of hydrogen atoms (or of other groups of atoms), 
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which are situated on the periphery of the molecule and which are bonded 
only by o electrons; the latter do not take part in the optical excitation 
of the molecule. In a certain sense, hydrogen atoms (CH; groups and 
others) play the same role as the external electrons of the atoms in rare- 
earth elements. They accept all influences of the surrounding molecules 
and protect the w electrons from perturbations. Thus the dependence 
of molecular wave functions for 7 electrons on molecular displacements is 
diminished. This leads to a decrease in matrix elements [(8.8)] and, 
therefore, to an increase in the lifetime of the corresponding molecular 
states also. Consequently, the probability of luminescence increases. 
Actually, all these substances display a strong luminescence, especially 
at low temperatures. If, upon light absorption, excitons with negative 
effective mass form, the lifetime of excited crystalline states may exceed 
that of the corresponding excited states of an isolated molecule. This 
last effect may be observed only in large monocrystals which do not 
contain inhomogeneities and which are at low temperature. 

Investigations by the Soviet school of luminescence of 8. I. Vavilov 
have established the scientific basis for classifying luminescence accord- 
ing to various characteristics. §. I. Vavilov’> proposed the division of 
luminescence into three types—spontaneous, forced, and recombinational 
—which are distinguished by their mechanism. This division is not 
linked with the aggregation state of the substance. V. L. Levshin”’ 
proposed dividing the luminescence of solids into two types, differing 
completely in their external characteristics and in their radiation 
mechanism. One fluorescence was named by him crystalline fluorescence. 
This fluorescence was observed only in a crystalline state; it did not have 
any direct link with the absorption spectrum; it possessed a complicated 
damping character and was accompanied by photoconductivity. The 
second fluorescence was called by V. L. Levshin fluorescence of discrete 
centers. It was observed in all states of aggregation; it has a direct 
connection with the absorption spectrum, possesses an exponential 
damping character, and is not accompanied by photoconductivity. 
V. L. Levshin relates fluorescence in crystals of aromatic compounds to 
fluorescence of discrete centers. 

The above-mentioned division is extremely essential and necessary in 
scientifically establishing the investigation of luminescence in solids. 
However, it seems that, on the basis of the present work and of certain 
experimental data relative to luminescence spectra of molecular crystals 
at low temperatures (works of I. V. Obreimov, A. F. Prikhotko, and 
their school), the following conclusions may be made: luminescence of 
molecular crystals may essentially depend on the crystalline structure of 
the solid. The luminescence of such solids often takes place during 
transitions to a normal state from excited crystal states, which corre- 
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spond to the occurrence in a crystal of free excitons. Such excited states 
may arise only in crystalline solids, and they depend on the lattice 
structure. 

Therefore, luminescence of discrete centers should not always be under- 
stood to be a luminescence which takes place completely in localized 
fluorescence centers. Sometimes excitation energy may execute a con- 
siderable path in a crystal, passing from one ‘‘center” to another prior 
to being radiated. It is necessary in this connection that a direct trans- 
mission of excitation energy from one center to another occur. This 
process of energy transmission affects relatively little the character of 
the luminescence. Therefore, if we do not take into consideration these 
small changes, we can sometimes identify such a luminescence with the 
luminescence of the center, as found in a free state. 

For crystalline fluorescence, excitation energy is transmitted by the 
electron over the crystal, and the change from a crystalline state to other 
states of aggregation may destroy the capacity for luminescence. 


CHAPTER IX 


COMPARISON OF THEORY WITH EXPERIMENT 


In some of the preceding chapters we pointed out the connection 
between results obtained and experimental data. In this chapter we 
shall consider supplementary experimental data. 

It must be noted here that available experimental material is often 
only preliminary and incomplete. The reasons for this are, on the one 
hand, the great experimental difficulties that arise when carrying out 
experiments of this kind, and on the other hand, the absence of theory, 
which often causes an insufficiently rational foundation of the experiment. 
For instance, many measurements were made in nonpolarized light and 
at relatively high temperature. 

We hope this work will serve as a certain stimulus to setting up addi- 
tional experiments which, on their part, will lead to a fuller develop- 
ment of the theory. 


35. Combinational Scattering (Raman) Spectra 


Despite the fact that the discovery of combinational scattering 
(Raman scattering) in 1928 by L. I. Mandelstam and G. 8. Landsberg 
was made on crystalline quartz, combinational scattering (Raman 
scattering) of solids up to the present time has had considerably less 
study than the scattering of liquids and gases. This is explained first 
of all by the great experimental difficulties of investigating combinational 
scattering (Raman scattering) of monocrystals. One must have trans- 
parent monocrystals of sufficient dimensions. The investigations must 
be conducted in polarized light and at low temperatures. Secondly, it 
seems to be explained by the prevailing, long-standing opinion that 
combinational scattering (Raman scattering) of crystals is not different 
from scattering of liquids and gases. 

For many vibrational frequencies, the differences between combi- 
national scattering (Raman) spectra of gases, liquids, and solids are 
actually very small. However, as follows from the investigation con- 
ducted here, significant deviations from this rule may be observed. 

As was shown in the preceding chapters, upon transition from a gaseous 
state to a crystalline one, a displacement of molecular frequencies, a 
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splitting (if, in a unit cell, there is more than one molecule or if, in a 
molecule, there are degenerate states and the symmetry of the crystal is 
lower than that of the molecule), and a widening of lines as a result of 
interaction with lattice vibrations take place. 

If the corresponding molecular vibration in an isolated molecule is not 
accompanied by a change of the dipole moment of the molecule, all 
changes enumerated above are very small and, except for the widening 
of lines, may often not be observed under present experimental conditions. 

Thus intramolecular vibrational frequencies, observed in combinational 
scattering (Raman scattering), may be divided into two groups: (1) vibra- 
tions noticeably affected by crystalline structure (these vibrations must 
be active in the infrared spectrum) and (2) vibrations not noticeably 
affected by crystalline structure. Sometimes, for one and the same 
molecule, vibrations may belong to both above-mentioned groups. Sev- 
eral researchers have already noted this circumstance. Thus G. Herz- 
berg, in a monograph devoted to infrared and Raman spectra of poly- 
atomic molecules, writes (Ref. 50, p. 584): ‘Worthy of our attention is 
the fact that in certain cases various vibrations of one and the same 
molecule give non-identical relative frequency displacements.” 

Table 16 shows frequency displacements for certain substances as 
observed in the combinational scattering (Raman scattering) of gas, 
liquid, and solid, as well as frequencies observed in the infrared spectrum, 
if these frequencies coincide with combinational scattering (Raman 
scattering) frequencies. Although the molecule C.H, possesses a center 
of symmetry and should not have coinciding vibration frequencies in the 
infrared spectrum and the combinational scattering (Raman) spectrum, 
nevertheless, owing to the chance coincidence of certain vibrational 
frequencies in the combinational scattering (Raman) spectrum, larger 
frequency displacements occur upon transition from gas to crystal than 

TaBLe 16. DispLacementT or C. S. (RaMAN) FREQUENCIES IN CONDENSED 


Puaszts rok MoLecutark VIBRATIONS, CORRELATED WITH 
SIMULTANEOUS INFRARED ACTIVITY 


Combinational scattering (Raman) ' 
=| COMPRDONEING ic ocutice 
Molecule econ : infrared 
= Liquid Solid source 
Gas », em7} frequency 
Av = vg — vt| Av = vg — vgs 
HCl 2886 101 118 2886.7 [78], [79] 
C8, 665 —1.5 (80) 
CoHe 992 Oy, jiill aivtact, “decors 
CoHe 3069 oe) ae Ca 3073 [81] 
CcH. 3099 9.0 10.0 3099 
CeH6 3045 10.0 11.0 3045 
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for vibrations that are only present in the combinational scattering 
(Raman) spectrum. 
Other experimental data were mentioned in Section 10. 


36. Infrared Spectra 


Studies of changes in infrared absorption spectra, with changes in 
the state of aggregation for a substance, encounter numerous experi- 
mental difficulties. There are comparatively few works devoted to these 
problems. 

From all that was stated above, it is clear that, in order to discover the 
corresponding changes, one must conduct experiments at the lowest 
possible temperature, with thin absorbing layers, and certainly in 
polarized light. We do not know of any works where all these conditions 
were fulfilled. 

Of great interest are the investigations undertaken by A. F. Prikhotko, 
I. V. Rodnikova, and K. Shabaldas to clarify the changes introduced in an 
infrared spectrum upon transition of certain substances to a crystalline 
state. Unfortunately, since these authors conducted their investiga- 
tions in the photographic infrared region, they were only able to observe 
spectra corresponding to high combinational frequencies. Besides, the 
investigations were carried out in nonpolarized light. Although these 
works were interrupted by the war and were not completely finished, 
the results published by I. V. Rodnikova are still of interest. 

The spectra of ethylene were studied at a temperature near the melting 
point. The molecule of ethylene CoH, belongs to the point symmetry 
group Ds, (the characters of irreducible representations for this group 
were given in Table 2). This molecule possesses 12 normal vibrations 
which are distributed among the irreducible representations of the mole- 
ewe’s symmetry group in the following manner: 


3A 1, + 2A 29 + Bi, + Au + 2B + Ag + 2Boy. 


The molecule has a center of symmetry; therefore, the vibrations active 
in the infrared spectrum (2B, + Ao + 282.) should not appear in the 
combinational scattering (Raman) spectrum. Experimental data cor- 
roborate the absence of coincidences in both spectra of gases. °° 

As has already been indicated in the preceding section, when there is 
absence of coincidences of vibrational frequencies in the infrared spectrum 
and the combinational scattering (Raman) spectrum, one must expect 
small changes in the combinational scattering (Raman) spectrum when 
the state of aggregation changes. The combinational scattering (Raman) 
spectrum of the ethylene molecule well confirms this rule. All combi- 
national scattering (Raman) lines observed in a gas are also observed 
in a liquid with very small displacements (Ref. 50, p. 326). 
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Vibrations observed in the infrared spectrum behave quite differently. 
In her doctoral dissertation, A. F. Prikhotko, on the basis of the analysis 
of spectra for a crystal, a liquid, and a gas in the photographic infrared 
region, concludes: ‘“‘The absorption spectrum of the ethylene crystal is 
remarkable in that it has a large number of satellites and a series of new 
absorption lines, which are absent in a gas and which appear in a liquid 
and mainly in a crystal.” 

Figure 15 represents schematically (according to data furnished by 
I. V. Rodnikova’s work®) absorption spectra of a gas (without rota- 
tional structure) and of a crystal. We see in the absorption spectrum 
of the crystal a large number of new absorption bands in comparison 
with the gas spectrum. A portion of these absorption bands must be 


Gas, T=20°C \ 


Crystal, T=— 180°C 


95 100 105 110 115 120 140 
7, cm7}, hundreds 


Fia. 15. Schematic representation of the absorption spectrum of the ethylene molecule 
in the photographic infrared region. (After Rodnikova.) 


attributed to the combination with small lattice vibrational frequencies. 
However, the genesis of other absorption bands is apparently linked to 
the splitting of absorption bands in a crystal; this splitting was analyzed 
in this work and is caused by intermolecular interaction. 

Since ethylene molecules do not have degenerate vibrations, it is quite 
impossible to explain the observed splitting (as is frequently done) by 
the removal of degeneracy in the crystalline field. 

According to X-ray studies, * an ethylene crystal at —175° belongs to 
the rhombic system with two molecules per unit cell. Therefore, a 
splitting of absorption bands into two components must be expected. 
Figure 15 clearly shows the splitting (noted in the previously cited work®) 
of the absorption band at 14,100 cm~ for a gas into two components in a 
crystal of 13,972 and 14,014 em (Av = 42 em~). These two com- 
ponents should have a different polarization. 

The methane spectrum in the photographic infrared region was investi- 
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gated at 20.4°K in the laboratory for crystals of UFTI [Ukrainian Photo- 
Technical Institute]. A. F. Prikhotko notes in her dissertation that the 
spectrum of the crystal is richer in lines than the spectrum of the gas. 
A. F. Prikhotko points out the removal of degeneracy in the crystalline 
field of degenerate molecular vibrations as a possible reason for this. 
In view of the incompleteness of this work (there are no measurements 
taken in polarized light), it is at this time still difficult to say whether 
one can explain all new absorption bands merely by the removal of 
vibrational degeneracy in the isolated molecule. Apparently this is not 
so, if only because, owing to the neutrality and absence of a static dipole 
moment in a methane molecule, it is not possible to expect noticeable 
absorption-band splittings as a result of the removal of degeneracy. 

Besides the splitting by a crystalline field of degenerate vibrations in a 
free molecule, other forms of absorption-band splitting in crystals were 
not known. Therefore, an attempt was made to explain every splitting 
by means of this first phenomenon. And, if it was not possible to do 
this, an attempt was made to explain the splitting by a different kind 
of polymerization and association of molecules. Without denying the 
possibility of molecular polymerization and association in some cases, 
we believe that splitting phenomena may frequently be explained by the 
theory developed in the present work. 

The difference in position and in relative intensity of vibrational bands 
in a liquid, a solid, and a vapor state is noted by other authors as well. 
Thus, in a recently published work of Richards and Thompson,* the 
infrared spectra of more than 20 substances (in a solid and liquid state) 
were investigated. These substances consisted of nonpolar molecules 
(in order to exclude effects connected with the association of molecules). 
The authors came to the conclusion that ‘‘the bands of a solid are sharper; 
many wide bands, existing in a liquid, split into two or more components 
upon transition to a solid state; in many cases, albeit not always, there 
are more bands in a solid state than in a liquid.” Unfortunately, the 
authors conducted their measurements at relatively high temperatures 
and in nonpolarized light. 

It would be very desirable to conduct a series of experimental investi- 
gations of crystalline absorption spectra at low temperatures and in 
polarized light, especially for crystals consisting of molecules that do not 
have degenerate vibrations in a free state. 


37, Absorption and Luminescence Spectra in the Visible and 
Ultraviolet Regions 


a. Naphthalene Crystal. Henri and de Lazlo® and de Lazlo®* investi- 
gated the spectrum of light absorption Gin the ultraviolet region) by 
naphthalene vapor. A fuller investigation, together with an analysis of 
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the vibrational structure of the spectrum, was made by A. F. Prikhotko.** 
When analyzing the absorption spectrum in a layer of gas 100 cm long, 
Prikhotko established three electronic transitions and the vibrational 
frequencies which combine with them. Table 17 lists the results 
obtained. 
TABLE 17. OBSERVED ELECTRONIC TRANSITION FREQUENCIES IN 


NAPHTHALENE, WITH STATE SYMMETRIES, POLARIZATIONS, 
AND COMBINING VIBRATIONAL FREQUENCIES 


Electr ae Trreducible Direction of Vibrational 
transition, < vibration of s - 
= representation : frequencies, cm~! 
em7! electric vector 
32,455 Bu, y 473, 703, 1002, 1181, 1415, 1449, and 
1603 
35,905 Bu 2 493 and 1397 
About 44,000 Bu Zz 


The irreducible representations corresponding to these electronic 
transitions, as well as the direction of electric-vector vibrations in the 
molecule, are shown as the basis of a comparison with this author’s 
theoretical calculation.44 The same results were obtained in a later 
theoretical work by Coulson®’ for the first two electronic transitions, 
which are the only ones that interest us further. 

The absorption spectra for solutions of naphthalene in hexane, studied 
by de Lazlo,®* furnish the same picture as the gas spectrum. The spec- 
trum of the solution is merely displaced in relation to the spectrum of the 
gas. The displacement of the first transition is unlike that of the second. 
For each group of lines in a gas, de Lazlo was able to show in the solution 
the corresponding absorption maximum. Figure 16 juxtaposes the 
spectrum of naphthalene gas and that of a solution in hexane (according 
to de Lazlo). Thus the interaction of the solvent molecules with the 
naphthalene molecules does not result in the emergence of new absorption 
lines, but merely displaces the spectrum. This was to be expected, 
since, aS a result of the absence of resonance interaction, the excitation 
energy of a naphthalene molecule in solution AH, will differ from the 
excitation energy of free molecule A#,, only by the energy difference for 
the interaction of an excited and a normal naphthalene molecule with 
the surrounding molecules of the solvent: 


AE, = AE, + D. 


A considerably more complicated picture is observed in a crystal. 
The absorption spectrum of naphthalene monocrystals at the temperature 
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of liquid hydrogen is now relatively well known, thanks to the works of 
A. F. Prikhotko.™*3* Here we shall merely indicate the most funda- 
mental features of the spectrum which are necessary for a comparison 
with theory. 


34 ~~ 6 7 
aie we ON 
B EF, 
A G) 123 ca: 
Vapor, pt. 1 


Solution in hexane 


31 33 35 37 39 41 43 45 
v, cm7}, thousands 


Fig. 16. Schematic representation of the absorption spectrum of naphthalene in the 
near ultraviolet. (After de Lazlo.) 


Fig. 17. Disposition of the two molecules of naphthalene in the unit cell of the mono- 
clinic crystal. 


Naphthalene crystallizes** in a monoclinic-prismatic system C%,, with 
a = 8.29 A, b = 5.97 A, c = 8.68 A, and 8 = 122.7°. In the unit cell 
there are two molecules whose disposition is schematically shown in 
Figure 17. The long molecular axis z forms, with axes a and 6 of the unit 
cell and with the perpendicular to planes ab, angles of 115.3°, 102.6°, and 
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28.7°, respectively. The intermediate axis of the molecule forms, with 
the same axes of the cell, angles of 71.2°, 28.8°, and 69.1°. 

Crystals obtained from vapors represent fine leaves whose plane 
coincides with the plane ad of the crystal. When studying the absorption 
of light, the latter is usually directed perpendicularly to the plane of the 
leaf. The vibrations of the electric vector for the light wave are in this 
case either parallel to axis a of the crystal (@ component) or parallel to 
axis 6 of the crystal (6 component). 

The absorption spectrum of very fine (fraction of a micron) naphthalene 
monocrystals at a temperature of 20.4°K represents the sum total of a 
large number of absorption bands; some of these bands (especially in the 
acomponent) have a very small width. Between bands of strong absorp- 
tion, a background of continuous absorption is usually observed. This 
background becomes stronger in proportion to the increase in the thick- 
ness of the crystal, especially in the short-wave portion, and, for thick- 
nesses beginning with 1 mm, turns into a continuous absorption beginning 
with a frequency of 31,476 cm—. New absorption bands simultaneously 
appear in the long-wave portion of the spectrum. 

It is to be noted that the study of the naphthalene absorption spectrum 
was usually limited to the study of absorption bands in the long-wave 
portion, and that of continuous absorption in the short-wave portion, of 
the spectrum. It was only work with very fine layers of monocrystals 
which enabled A. F. Prikhotko to obtain absorption bands in the short- 
wave portion of the spectrum. 

The first electronic transition in the molecule (32,455 cm—), together 
with its vibrational structure in the crystal, yields corresponding series 
in both spectral components (a and 6); in relation to the absorption in a 
gas, these series are displaced by 500 cm~ toward long waves. These 
absorption bands must be determined by the excitation in the crystal 
of “localized’”’ excitons. According to (7.22), the energy of the crystal 
is in this case equal to 


AEL, = ABY, + 8° — 8°, 


since Ae, (k) ~~ 0 for “localized” excitons. In this manner, the change in 
lattice energy upon excitation of a “localized”’ exciton will be determined 
from the condition 

gle — g° 


ps -1 
Oahe 500 em7!. 


However, to.the first electronic transition in the molecule in a crystal, 
there correspond weak and extremely long wave absorption bands: 
29,944 cm~ in a b component and 29,931 cm! in an acomponent. The 
splitting (Av = 18 cm7?) and the sharp polarization of these split com- 
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ponents indicate that both of them appear as a result of the formation of 
free excitons in a crystal. In this case, by virtue of (5.34), the crystal 
excitation energy will be equal to 


AE} 


ABEL + D+ is (ADL oy + M3501), b component, 


AE} = AEL + Di+ 5 (Aft o — Mio, a component. 


If we equate the frequency difference to the experimental magnitude 


AE} — AE} = 2) M1, 9, = 13 - Qehe, 
we obtain 
by Mio o, = 8 X 10-* ev. (9.1) 


The displacement of the absorption band in relation to the gas is equal 
to 2518 em7; consequently, 


D+" Mh o = —0.31 ev. (9.2) 


The ratio of the intensities of both components of the splitting, according 
to (5.35), must. be 
I(1) 


73) = cot? a = 7.2. 


aie 


In A. F. Prikhotko’s experiments, the intensities were not measured; 
however, it is noted that the line in the a component (29,931 cm-) is 
very faint; in the b component (29,944 cm), it is sharp and more 
intense. 

Utilizing Robertson’s data concerning the dimensions of the naphtha- 
lene unit cell and the disposition of molecules therein, one can roughly 
estimate the sum of matrix elements [(9.1)] in theory, if one limits oneself 
to taking into account the interaction of the nearest molecules and if, for 
the matrix elements, one uses the approximate expressions (3.20). 


YY! Maaor = 3.6 X 104| fo! (te? dr]. (9.3) 


Inasmuch as (9.3) is positive, it indicates that the disposition of split 
components is qualitatively in agreement with the experiment. 

Equating the calculated value of (9.3) with the experimental one of 
(9.1), we obtain for the square of the matrix element of the transition in 
a molecule the value 


[fol (r)¢° del? = 2.2 x 10-8 em*. (9.4) 
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Using (9.4), one can estimate (by order of magnitude) the oscillator 
strength for the transition in a molecule 


f~001. 


To each molecular excitation consisting of the first electronic transi- 
tion, there must, in a crystal, correspond a split absorption band; one 
component of the latter will be in the b spectrum, the other, considerably 
less intense, will be in the a component. Insofar as the matrix elements 
of the electronic-vibrational transition and D? differ from the correspond- 
ing magnitudes of a purely electronic transition, the magnitude of the 
splitting and the displacement of corresponding absorption bands will 
also be different. Therefore, the distances of the vibronic bands from 
the corresponding purely electronic absorption bands will not be rigorously 
equal to the vibrational frequencies observed in a free molecule. This 
actually occurs for absorption bands observed in a b component. 


Vibrational frequencies Vibrational frequencies 
in @ gas, cm} in a crystal, em~* 

473 400 

703 724 

1002 1051 
1181 1282 
1415 1383 
1603 1463 


In a crystal, all these bands are of very low intensity in the 6 com- 
ponent; in the a component, there are very weak absorption traces, the 
position of which is difficult to measure. 

To the second transition in the molecule (35,905 cm~), a large number 
of wide absorption bands correspond in a crystal; among them, an 
absorption band of 33,736 cm—! stands out by its extremely high intensity; 
it is comparatively easy to establish from the latter band combinational 
repetitions with frequencies of 518 and 1398 ecm-'. On this basis, 
Prikhotko identifies precisely this absorption band, with the frequency 
35,905 cm—}, observed in naphthalene vapor. A characteristic peculiarity 
of all these absorption bands is their large displacement in comparison 
with corresponding bands of the gas spectrum: 


62° — &° = —2169 em—!: 2rhe. 


A crystal absorption-band displacement large in comparison with gas 
absorption bands indicates a significant change in interaction forces upon 
excitation of molecules in a crystal. Therefore, absorption is always 
linked to lattice vibrational excitation. The corresponding absorption 
lines are wide. 
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To molecular excitations involving the second electronic transition as 
well as molecular vibrations, there is a corresponding free exciton excita- 
tion in a crystal. This results in a corresponding splitting of crystalline 
absorption bands. ‘Table 18 indicates the split frequencies (in reciprocal 
centimeters) observed in an absorption spectrum corresponding to the 
second transition in a molecule. 


TaBLy 18. FREQUENCIES or LINES oF THE SECOND ELECTRONIC TRANSITION 
IN NAPHTHALENE IN THE VAPOR AND THE CRYSTAL 


Crystal, em7! 
Vapor, 
om b component | a component Magnitudes 
of splitting 
35,905 31,063 31,050 13 
36,398 31,620 31,474 146 
37 ,302 32,255 32,227 28 


Figure 18 represents (according to A. F. Prikhotko’s work) a portion 
of the absorption spectrum (at a temperature of 20.4°K), of a naphthalene 
monocrystal, of light with two polarizations (a is the direction of the 
electric vector parallel to axis a of the crystal, and 6 is the direction of 
the electric vector parallel to axis b of the crystal). The splitting of the 
two absorption bands and the sharp polarization of split components are 


31474 31955 32227 cm! 


Il a-axis 


|| b—axis 
31620 31955 32255| 32673 32948 cm”! 
32411 


Fia. 18. Polarized light absorption spectra of monocrystals of naphthalene at 20.4°K 
for two orientations of the elective vector of the light wave. (After Prikhotko.}3:32) 


clearly visible on the photograph. The bands, whose positions coincide 
in both components, correspond to the excitation of “localized” excitons 
in the crystal. 

It would be interesting to investigate the absorption spectra of the 
naphthalene monocrystal with the ray directed along the monoclinic 
axis and with only one of the two split components present in the spec- 
trum. The intensity of the component should vary from zero to a certain 
maximum magnitude for a Nicol turn of 90°. 
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All the most intensive absorption bands in the naphthalene crystal 
are linked to the second electronic transition in the molecule. For such 
excitations, the interaction of the molecule with its neighbors in the 
crystal varies considerably. Asa result of this, the molecular equilibrium 
positions are displaced, and the excitation energy passes over to lattice 
vibrations. These excited states will possess a short lifetime, shorter 
than the time necessary for radiation. Luminescence from states which 
correspond to absorption with formation of ‘‘localized’”’ excitons of the 
second electronic transition in a molecule is therefore not possible. 

Since the free exciton states have, in this case, less energy than the 
“localized” exciton states, a certain portion of “‘localized” excitons will 
pass to a free exciton state (see Chapter VIII). The latter free excitons 
will give off their energy in the form of radiation. This process is 
apparently responsible for the intense luminescence of naphthalene 
crystals. The luminescence begins with frequencies of 31,062 and 
31,044 cm! (weaker), and these frequencies correspond to the excitation 
energy for free excitons of the second transition in the molecule. Accord- 
ing to a detailed investigation conducted by I. V. Obreimov and K. G. 
Shabaldas!® on the luminescence spectrum for the naphthalene crystal 
at a temperature of 20.4°K, the most intense luminescence spectrum 
begins with the above-indicated lines and extends further toward lower 
frequencies. This is because the transition to a normal state is accom- 
panied by the excitation of intramolecular vibrations and thus lower 
frequency lines. 

Light absorption bands corresponding to the excitation of the first 
electronic transition in the molecule and to its combinations with molecu- 
lar vibrations are of low intensity. The frequencies of this band are not 
much displaced (500 cm-) in relation to the frequencies of corresponding 
excitations in a free molecule. This small displacement indicates a small 
change in interaction forces (this is also corroborated by the comparatively 
small width of absorption bands); it also bears witness to a lifetime suf- 
ficient for luminescence to occur by transition to the ground state. 

To an absorption band of 31,960 cm (this band is dependent on 
“localized”? exciton excitation), there corresponds in the luminescence 
spectrum a doublet of 31,960 and 31,965 cm—. Connected with this 
doublet is a series of lines located at periodic distances from it. The 
distances are equal to the periods of progressions starting with the line 
of 31,062 cm~-'. We also observe in the luminescence spectrum a fre- 
quency of 29,941 cm—!, which appears upon transition to the ground 
state from an excited state corresponding to free excitons of the first 
electronic transition in the molecule. 

b. Anthracene Crystal. The anthracene crystal belongs to the mono- 
clinic-prismatic group of symmetry C$, with two molecules per unit cell; 
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the disposition of these molecules recalls that of the molecules in a 
naphthalene crystal. The first electronic transition in the molecule has 
an electric moment direction parallel to the intermediate axis of the 
molecule. 

The absorption spectrum of anthracene crystals at a temperature of 
20.4°K was investigated by I. V. Obreimov and A. F. Prikhotko.!7 
Figure 19 is a diagram borrowed from the work cited and representing 
the absorption for components a and 6 of the spectrum. From this 
diagram, one can see immediately that the portion of absorption bands 
starting with an electronic transition of 25,370 cm- is simultaneously 
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Fic. 19. Schematic representation of a and b components of the absorption spectrum 
of anthracene at 20.4°K. (After Obretmov and Prikhotko.) 


present in both spectral components (absorption bands in component d are 
more intense). All these absorption bands correspond to the excitation 
of “localized” excitons in the crystal. The general displacement of 
absorption bands in comparison with the gas spectrum is very great 
(~2246 em). This indicates a significant change in interaction forces 
upon molecular excitation. 

Besides absorption bands coinciding in a and b components, there is 
in the b component a series of absorption bands (as a rule, considerably 
narrower), which are situated primarily on the short-wavelength side. 
These are absorption bands of 24,700, 24,962, 25,046, 25,367 cm—, etc. 
(Ref. 17, Table 4). These absorption bands are explained by the excita- 
tion of free excitons in a crystal. 

As was shown in Section 22, in crystals containing two molecules per 
unit cell, two types of free excitons (u = 1, 2), with a wave number equal 
to the wave number of the light wave, may be excited under the action 
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of light. The movement of these excitons is of various velocities. The 
velocity of movement for excitons of one type may be so small that the 
formation of free excitons of this type becomes impossible. The fact that, 
in an anthracene crystal, one observes absorption in one component only 
indicates that this absorption is dependent on the excitation of free exci- 
tons of one type only. 

The luminescence spectrum of an anthracene crystal at a temperature 
of 20.4°K was studied by I. V. Obreimov, A. F. Prikhotko, and K. G. 
Shabaldas.18 The entire luminescence spectrum begins with a fre- 
quency of 24,909 cm—'; therefore, it corresponds to a transition from the 
excited state formed by a free exciton to the ground state of the crystal. 

c. Naphthacene Crystal. The naphthacene molecule belongs to the 
point symmetry group Da». The molecule does not possess degenerate 
energy levels. According to Bethe, it is therefore impossible to expect 
level splittings in its crystalline state. Scheibe and Kandler® investi- 
gated the spectrum of crystalline naphthacene and showed that the 
absorption spectrum differs for two mutually perpendicular polarizations. 
In one component, an absorption band of 525 my is observed; in the other 
component, bands of 475 and 445 my, and others, are observed. In the 
polycrystal reflection spectrum, absorption bands of 484, 462, 495, 
and 526 my are found.®? The same absorption bands are observed upon 
absorption in a fine layer of naphthacene sublimed on glass (see solid 
curve in Figure 20). The dashed curve in Figure 20 shows the optical 
absorption spectrum for a solution of naphthacene in a-bromnaphthalene. 
It is evident from this figure that the spectrum of the solution does not 
contain two extremely long-wave ‘“‘crystalline” absorption bands. 

The presence of dichroism was also confirmed in Jelley’s more recent 
experiments (Ref. 90, p. 338) with naphthacene crystals obtained upon 
crystallization from a-bromnaphthalene. 

The most detailed studies of naphthacene-monocrystal absorption 
spectra in polarized light and their comparison with vapor spectra were 
carried out by A. Y. Hichis in her candidate’s dissertation. The results 
obtained by her are listed in Table 19. 


Tas_e 19. Frequenciss or LINEs oF THE First ELECTRONIC TRANSITION 
in NAPHTHACENE IN THE VAPOR AND THE CRYSTAL 


Absorption spectrum Frequency, em™ 


Naphthacene vapor (15-cm path), 20°C......... 22 ,200 23,625 25,000 
Monocrystal of thickness of order of tenths of 
micron, —190°C: 
Vibration of electric vector || b.............. 19,235 20,715 22,320 
Vibration of electric vector Lb............. 19,810 21,145 22,545 
Magnitude of splitting..................... 575 430 225 


136 DAVYDOV: THEORY OF MOLECULAR EXCITONS 


| J * 
400 450 500 550 $00 
Wavelength, mu 


Fie. 20. Schematic representation of the first absorption band of naphthacene. 
Continuous curve, absorption of a solid layer; Dashed curve, absorption of a solution in 
a~bromnaphthalene. 


Figure 21 shows a photograph of absorption spectra for vapors and for 
a monocrystal in polarized light, according to the data of A. Y. Eichis 
(a is the absorption spectrum for vapors; 6 is the absorption spectrum 
of a naphthacene monocrystal for two polarization components). A. Y. 
Eichis also determined the dispersion curves of monocrystals in polarized 
light. From these curves, she was able to estimate the oscillator strength 
of the first electronic transition in a crystal. The oscillator strength 
proved to be equal to 0.2 for the spectrum with an electric vector of light 
vibration parallel to axis b, and equal to 0.1 for a vibration perpendicular 
to axis 5. 

It follows from Table 19 that, relative to the spectrum of vapors, the 
crystalline spectrum is displaced by approximately 2600 cm—'. Besides, 
each band is split into two sharply polarized bands, each of which is found 
in only one component of the spectrum. The magnitude of the splitting 
is indicated in Table 19. The sharp polarization of absorption bands in 
a naphthacene monocrystal indicates that, upon light absorption, free 
excitons are formed. The formation of free excitons is facilitated by the 
relatively great oscillator strength for this transition, as well as by the 
large molecular size, which causes a large time of deformation in the 
crystal. 
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Fic. 21. Absorption spectra of naphthacene in the vapor state and as a monocrystal 
in polarized light for two polarization components. (After Eichis.) 

d. Phenanthrene Crystal. The phenanthrene molecule belongs to the 
point symmetry group C's. Consequently, it should not have degenerate 
energy levels. 

The spectrum of polarized light absorption at low temperatures by the 
phenanthrene monocrystal was investigated by I. V. Obreimov and A. F. 
Prikhotko.!2!6 The most detailed comparisons of phenanthrene spectra 
in various states of aggregation and the measurement of oscillator 
strengths for normal electronic transitions in the molecule were carried 
out in the candidate dissertation of A. Y. Eichis. In this work, it was 
shown that the first absorption region, which is usually observed in thick 
crystals and which starts with a frequency of 25,967 cm™!, owes its 
existence not to phenanthrene itself, but to an admixture of anthracene 
in the phenanthrene. However, the absorption of phenanthrene itself 
at a temperature of —253°C begins with a frequency of 28,649 cm—! for 
the case of light wave electric vector vibrations parallel to the b axis of 
the crystal and with a frequency of 28,604 cm~ for vibrations perpen- 
dicular to the 6 axis (the monocrystal was several microns thick). 

The phenanthrene monocrystal spectrum is displaced by approxi- 
mately 440 cm~! with respect to the spectrum of the vapor. To each 
electronic-vibrational absorption band in a vapor, there correspond two 
bands in a crystal. The distance between the first two bands, which 
correspond to a purely electronic transition and which are observed in 
various spectral components (a and 6), is equal to 45 cm-!. For elec- 
tronic-vibrational bands, these distances are still smaller. 

The small magnitude of the splitting (smaller by an order of magnitude 
than in naphthacene) indicates a small oscillator strength for correspond- 
ing transitions in the molecule. This deduction is confirmed by the 
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direct measurements carried out by A. Y. Eichis for oscillator strengths 
of electronic transitions in a phenanthrene molecule. IT*or the first 
electronic transition, the oscillator strength was approximately equal to 
0.01. (For naphthacene, as we previously noted, the oscillator strength 
of the first electronic transition was greater by an order of magnitude.) 

e. Benzene and Its Afonoderivatives. The region of 36,000 to 45,000 
em! is the most thoroughly investigated region of benzene absorption. 
The absorption spectrum of benzene vapors in this region consists of a 
complicated system of bands beginning with a frequency of 37,482 em. 

The benzene molecule belongs to the point symmetry group De. 
The above-indicated absorption region corresponds to the forbidden 
transition Ai, > Bo, which becomes permitted owing to an excitation 
with e-type vibrations. The frequency of the electronic transition 
itself is assumed to be equal to 38,089 cm—}.*! 

The absorption spectrum of solid benzene was studied by Pringsheim 
and Kronenberger at temperatures of 93 and 20.4°K.":9% In comparison 
with the absorption spectrum of vapors, the absorption spectrum of the 
crystal is displaced by 253 cm~'. There are several differences in the 
distribution of intensity and in structural details of the crystal spectrum. 
Low frequencies belonging to lattice vibrations occur. The luminescence 
spectrum®*?.* is also very similar to the luminescence spectrum of benzene 
vapors. The displacement relative to the luminescence spectrum of 
vapors is 240 cm7?. 

The small change in the absorption and fluorescence spectra upon 
transition to a crystalline state is in agreement with the forbidden 
character of the corresponding electronic transition. Since the transi- 
tion is forbidden, the resonance interaction is small; as a result, the 
corresponding changes are also small. A more detailed comparison of 
theory with new experimental data was conducted in Section 14d. 

If, for one hydrogen atom in benzene, we substitute an atom of Br, Cl, 
etc., symmetry Ds, becomes symmetry C2, and forbidden transition in 
benzene Ai, —> Bo, will correspond to transition A, — B, in the mono- 
derivative benzene. This transition is permitted in group C2, and the 
electric moment for the transition lies in the molecular plane in a direction 
perpendicular to the bond of the carbon substituent.*! Therefore, in 
the case of benzene derivatives, one should expect considerable changes 
upon transition from a gaseous state to a crystalline one. 

Absorption spectra of benzene derivatives (bromobenzene, iodobenzene, 
and others) at a temperature of 20.4°K were studied in the candidate 
dissertation of E. M. Bronstein under the direction of A. F. Prikhotko. 

The absorption spectrum of gaseous bromobenzene in the ultraviolet 
consists of two regions:*>-*6 (1) absorption bands (36,334 to 40,181 em) 
and (2) continuous absorption regions starting with 41,781 em71. 
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The absorption spectrum for a solution of bromobenzene in hexane%’ 
is displaced into the long-wave region by 250 cm~ with respect to the gas 
spectrum and, as was also to be expected, differs but little from the latter. 

In the spectrum of solid bromobenzene at 20.4°K, Bronstein observed 
absorption in the regions 28,700 to 31,900 em~! and 34,900 to 42,200 cm-!. 
Absorption bands contained in the range of 28,700 to 31,900 cm— do not 
exist in the spectra of the vapor and of solutions. Figure 22 is a diagram 
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Fic. 22. Schematic representation of the near ultraviolet absorption of bromobenzene. 
(After Bronstein.) 
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Fra. 23. Schematic representation of the near ultraviolet absorption of iodobenzene. 
(After Bronstein.) 


of vapor and crystalline bromobenzene absorption according to Bron- 
stein’s work. 

The presence of a new absorption band in the crystal in comparison 
with the vapor was also discovered by Bronstein for iodobenzene, chloro- 
benzene, and several other benzene monoderivatives. Figure 23 is a 
diagram of the absorption of iodobenzene vapors and crystal according to 
Bronstein’s work. 

The appearance of new absorption bands in all benzene monoderivatives 
upon transition from a gaseous state to a crystalline one is connected with 
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the action of intermolecular forces in a crystal. The considerable 
absorption-band displacement, as well as the width of the absorption 
bands, points to the presence of great interaction forces which, apparently, 
possess a resonance character. 

For a fuller clarification of this problem, additional studies of gases in 
large paths and of crystals at low temperature in polarized light are 
needed. 


CHAPTER X 


APPLICATION OF EXCITON CONCEPTS TO THE 
CALCULATION OF THE ENERGY STATES OF 
SEVERAL POLYATOMIC MOLECULESt 


38. Dependence of the Frequency of Light Absorption in 
Para-polyphenyls on the Number of Phenyl Groups 


At the present time, the methods of calculating energy states of 
polyatomic molecules are definitely unsatisfactory. Usually, three 
approximation methods are applied: (1) the valence bond method, (2) 
the molecular orbital method, and (3) the antisymmetrized molecular 
orbital method. In connection with the study of complicated molecules, 
the first method proves to be so cumbersome as to be practically inap- 
plicable. For example, in the calculation on the naphthalene molecule, 
one must solve an equation of the 42d degree; on anthracene, a 429th- 
degree equation, etc. The method of molecular orbitals is significantly 
simpler and is therefore applied more frequently. However, in this 
method, when reducing the polyelectronic problem to the problem of the 
motion of one electron in the coordinated field of the remaining electrons, 
an essential assumption is made concerning the smallness of the change in 
interaction magnitude of one electron with all the remaining ones when 
the state of the given electron changes. Furthermore, in the molecular 
orbital method, the exchange interaction of electrons is ignored. An 
effort to eliminate the latter simplification was made in the method of 
antisymmetrized orbitals, in which antisymmetrized wave functions are 
utilized. However, it seems to us that one must not believe that the 
same relationship exists between the molecular orbital methods with 
antisymmetrized and nonantisymmetrized functions as exists between 
the methods of self-consistent fields with and without quantal exchange, 
worked out for atomic systems by Academician V. A. Fock. Anti- 
symmetrization of wave functions in an atomic system always leads to 
improved results. However, the introduction of antisymmetrized molecu- 
lar orbitals for calculations on complicated molecules is usually connected 
with the separate study of large electronic and ionic interactions in lieu 


t The results of this chapter have been previously published.%:99 
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of study of small interactions of electrons of a given atom with the remain- 
ing neutral atoms. The mistakes introduced thereby are so large (an 
incorrect use of the perturbation method) that at times they exceed the 
improvement obtained by virtue of the antisymmetrization of functions. 
Furthermore, the result of the calculation is expressed in the form of a 
combination of a large number of integrals, each of which can only he 
approximately computed (the integrals are extremely complicated, and 
the functions entering into them are known with insufficient accuracy). 
This may also lead to large errors and, besides, an extreme complication 
of the calculations. 

The calculation on the relatively simple molecule benzene may serve 
as an illustration of what we have just said. As A. London has shown 
(Ref. 48, p. 418), the refinement of the calculation by means of introducing 
integrals referring to nonneighbor atoms leads not to better but to 
decisively worse results (the levels rise by 3 ev). 

On the basis of this method, M. A. Kovner and L. M. Sverdlovt 
calculated the energy levels of the fulvene molecule and established the 
fact that, by applying perturbation theory, the first approximation 
changes even the sequence of levels of the zeroth approximation of the 
theory. This clearly shows that the method is mathematically not well 
founded. 

We can apply the calculation methods for energy states, used by us for 
molecular crystals, to molecules consisting of a certain number of identical 
groups of atoms if the bonds between these groups are weaker than the 
bonds between atoms making up the group. 

We must note that, when calculating the optical behavior of com- 
plicated molecules by means of this method, it is sufficient merely to 
require that the optical electrons should be little involved in the binding 
between the atomic groups. In this connection, the length of the bond 
involving the other electrons, for example, the o electrons in the case of 
aromatic compounds, is of no concern whatsoever. 

If the energy levels of one such group of atoms are known, then one 
can show how they change upon formation of a molecule consisting of a 
certain number of such groups. 

As is well known, a large class of organic dyes includes compounds 
whose absorption spectrum depends upon the presence of several identical 
groups of atoms (chromophoric groups) occurring in a state of mutual 
conjugation. In this regard, the length of the conjugation chain is one of 
the important factors determining the color. Depending upon the type 
of binding between chromophoric groups, we observe dependence, in one 


{ This paper was presented at the Sixth Spectroscopy Colloquium, Kiev, 1948. A 
summary of the paper may be found in the proceedings. 
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or another way, of the wavelength of the absorbed light on the number of 
conjugated groups. The study of this dependence is of great interest. 

In this section, we shall give the quantum-mechanical deduction of the 
dependence of the wavelength and intensity of light absorption by para- 
polyphenyls on the number of phenyl groups in the chain. The mole- 
cules of para-polyphenyls represent chainlike linear formations with the 
phenyl group as a structural element; therefore, they may be treated as 
small unidimensional crystals of finite length. 

Let H; be the energy operator for the kth phenyl group in the para- 
polyphenyl formed by N phenyl groups, and let Vim be the interaction 
operator between the Ath and the mth groups. Then the energy operator 
for the molecule may be written in the form 


N 
H= ¥ H+ ¥ Vin. (10.1) 
k=1 


k<m 


Let us further assume that the normal state of the m electrons of the 
carbon atoms in the kth phenyl group is characterized by the wave 
function ¢,, depending on the coordinates and spins of the 7 electrons; 
the energy of the normal state is EZ; the energy of the lowest excited 
state is Ei, and the corresponding wave function is g,. The positions of 
the phenyl! groups in the molecules will be considered fixed. The normal 
state of the whole molecule will be determined by the wave functions 


&° = [(@N)!]-* } (—1) PY, (10.2) 
where VY = gigs -+- gy. The excited state of the molecule will be 
determined by the function 

& = (NYY axe (10.8) 
where 
xn = [NI S (-1)P Ye, Va = pn OD 


Using (10.1) and the wave functions [(10.2) and (10.3)], we can calculate 
the energies of the normal and excited states of the molecule; their dif- 
ference determines the energy of excitation of the molecule: 


AE =AEntDt+e (10.4) 


where AE,, is the excitation energy of one phenyl group; the magnitude 
D characterizes the difference in interaction energy of the excited and 
normal phenyl group with all remaining (NV — 1) groups in the molecule. 
In this interaction, an essential part will be played only by the interaction 
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between neighboring groups, since the interaction between nonneighbor 
groups will be significantly less. Therefore, the magnitude D will not in 
the first approximation depend on the number of phenyl groups in the 
molecule, and ein (10.4) is determined by the system of linear equations 


YY a.Mnm — m= 0, m=i,2,...,N, (10.5) 
where Af,,m is the matrix element for interaction between the zth and the 
mth molecules. If we limit consideration only to the interaction between 
neighboring groups, and if we introduce the designation Mnn»vi = M, the 
system of equation (10.5) will attain the following form: 


—ae + aM =0 
all — ase + aslf = 0 


(10.6) 


The system of equations (10.6) possesses nontrivial solutions only in 
the case where the determinant of the Nth order, Ay, composed of the 
coefficient of these equations, is equal to zero. One can easily see that 


_ sin (VN + 1)6 


Ay a (10.7) 
€ 
where 2cos 6 = Ww (10.8) 
The equation Ay = 0 will have N roots: 
al 
= y = 2 - 
A; N + 1 l 1, e] ’ N 
Using (10.8), we obtain 
i 
=2 Liceee 
e¢ = 2M cos Vad (10.9) 


For the energy of excitation of the molecule, determined by (10.4), we 
obtain 

al 
N+ 
where A = AE, + Dz 


AE, = A + 2M cos (10.10) 


From (10.9), it follows that the dependence of the lowest [excited] 
energy level of the molecule on the number of phenyl groups in the mole- 
cule will be expressed by the relation (since M < 0) 


_=A—-2 ea eee 
AE ain = A |M| cos Wad (10.11) 
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Consequently, as the number of phenyl groups increases, the red absorp- 
tion limit is displaced toward longer wavelengths. As the length of the 
chain increases, the displacement of the absorption band regularly 
decreases. 

Substituting « [(10.9)] in the system of equations (10.6), we can deter- 
mine the nth coefficient of the linear combination [(10.3)] and, con- 
sequently, the wave function of the excited state as well. After the 
corresponding normalization, we obtain 


N 
-14 
es _ 9 mln . aln ; 
© ee sin? = | y Xn SiN a5 =e (10.12) 
1 n 


In order to determine the oscillator strength for the transitions to the 
ith excited level in the molecule, the matrix element of the transition 
must be calculated: 

Qi = fH * (8° dr. (10.138) 


Using (10.12) and (10.2), we obtain 


N 
~ 
_ y : an ) . ant 
Q} Qa sin ap | ae: ; 


0 cot [rl/2(N + 1)] 
"  /2/(N + 1) 


0, if J is even, 


, if lis odd, 


where Q,, is the matrix element of the transitions in one phenyl group. 
Consequently, the oscillator strength f; of transitions to the /th level will 
be 


2 1 Prerae 
fi= (N + 1) |Qmn|? (ae vy cot? Beeaai if [is odd, 


0, if lis even. 


In the case of large N, the oscillator strength of the transition to the first 
excited level will change in direct proportion to N: 


Si = const »|Qm|2N, (10.14) 


where v is the frequency of the transitions reciprocal centimeters. 
Experimental investigations of the absorption of para-~polyphenyls were 
carried out by Gillam and Hey.!°** They determined the wavelength of 
the maximum of absorption in solutions of five para-polyphenyls (NV = 2, 
3, 4, 5, 6) in chloroform, and of four para-polyphenyls (N = 2, 3, 4, 6) in 
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hexane. Utilizing for biphenyl (NV = 2) and for hexaphenyl (N = 6) the 
values for the absorption of frequencies found by them, we can obtain A 
and .f (in reciprocal centimeters); for the solution in chloroform, 
A = 50,060, 1 = 10,800; for the solution in hexane, A = 50,900, 
AT = 10,200. 
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Fic. 24. Dependence of frequency of absorption on number of pheny! groups in linear 
para-polyphenyls. 


Figure 24 shows, according to formula (10.11) for the values of A and 
MM obtained, the curved or nonlinear dependence of the frequency for the 
absorption of light by para-polyphenyls on the number of phenyl groups. 
Curve A is for solutions of para-polyphenyls in chloroform; curve B is for 
solutions in hexane. The open circles correspond to the experimental 
values obtained by Gillam and Hey for solutions in chloroform and the 
solid circles to those for solutions in hexane. 

The results obtained in this section can be used when studying the 
optical behavior of large molecular chains that are formed in the poly- 
merization of several molecules. 


39. Spectrum of Light Absorption by the Biphenyl Molecule 


a. Energy Levels of the Molecule. The energy states of the biphenyl 
molecule were recently calculated by A. London.§* London used the 
method of antisymmetrized molecular orbitals. London’s calculation is 
exceptionally complicated and is not without basic flaws, which we 
already have mentioned in Section 38. In that section we used (just as 
London did) Goeppert-Mayer and Sklar’s®’ data for calculations on the 
benzene molecule, and we computed the change of the energy spectrum of 
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two benzene molecules (more precisely, of pheny] groups) when they 
form a biphenyl molecule. 

The lowest energy levels of the aromatic molecules are determined by 
the z electrons of carbon atoms. Therefore, we shall consider the phenyl- 
group spectrum to coincide with the benzene spectrum. 

According to data from X rays! and from geometrical analysis, the 
biphenyl molecule is planar and belongs to the point symmetry group 
Dx». Figure 25 shows the schematic location of carbon atoms in the 
biphenyl molecule. 


Fig. 25. Coordinate axes for the biphenyl molecule. 


Let us designate by H, the energy operator for the phenyl group a, and 
by Hy the energy operator for the second group; then 


A. = Ee!, 
Hb = Hee, 
where the #! are the wave functions of the six 7 electrons of each phenyl 
group. 
If V is the operator for the interaction between two groups, then the 
equation for determining the stationary energy states of a biphenyl 
molecule will have the form 


(H. + Hy + V)D = ED, (10.16) 


where D is the wave function of the biphenyl molecule. 

In the zeroth approximation, the product of wave functions of the 
normal state of the benzene molecule (of the phenyl groups) 24? may be 
taken as the wave function of the normal state of the biphenyl molecule. 
The spin function for the normal singlet state of the biphenyl molecule 
may also be represented in the form of the product of the corresponding 
spin states for the combining molecules a and b: 


(10.15) 


Xo = KeXo- 


The complete antisymmetrical wave functions of the normal singlet 
state of the biphenyl molecule may now be written in the form 


Do = (12)-# Y (—L) P8260 x0, (10.17) 


148 DAVYDOV: THEORY OF MOLECULAR EXCITONS 


where the summation Is carried out over all possible permutations of 12 
electrons, and where P, is one of the possible permutations. 

The energy of the normal state of the biphenyl molecule in the first 
approximation will be equal to 


E= EU+ B+ f |®3(I) |?V\@s (1D) |? dr 
— Si f #o*(1)ee(I1) Veg*(11) a9 (1) dr. (10.18) 


Here ®9(I) and # (II) represent the wave functions of the molecules a 
and b for a certain position 1,2, . . . , 6 of the electron in the molecule a 
and 7, 8, ..., 12 of the electron in the molecule 6; @2(1I) and 4?(I) 
correspond to the states, differing from the preceding ones in the permuta-~- 
tion between the molecules a and 6 of any pair of electrons; the summation 


¥ is carried out over all possible permutations of the pairs of electrons 


Z 
between the molecules a and 6 (the number of terms in this sum is equal 
to 36). 

The lowest /th excited state of the biphenyl molecule will correspond to 
the excitation of only one of the ‘‘benzene” molecules in the /th state. 
Therefore, the wave function of the zeroth approximations, corresponding 
to the Ith excited state of the molecule, will be written in the form 


D! = (2)~“leaki + cokgl, (10.19) 
where 
kh = (12!) )) (—1) PBs xtxe, 
a (10.20) 
kk = (12!)-# Y (—1)"P,o38hxoxd. 


Using (10.19), we shall, with the aid of (10.16), obtain two homogeneous 
linear algebraic equations for determining ¢c, and ¢ in (10.19): 


—eée, + M'e, = 0, M'e, — éce, = 0, (10.21) 
where 
—é = BR + ER + f |e) /(?V/as(1D)|? dr 
— ¥ f at ee(Il) Veg (11) a3(1) dr — E, (10.22) 


M' = f #*(1)e9(1) vao*(I1)a4(I1) dr 
~—} f Des) Vap* (11) 44(1) dr. (10.28) 


In order to obtain the nontrivial solutions of the system of equations 
(10.21), it is necessary to have 


= +M1 (10.24) 
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Thus to each energy term of the phenyl group there will correspond 
two terms in the biphenyl molecule, depending on whether the plus or 
minus sign is taken in (10.24). 

Henceforth, only the excitation energy of the biphenyl molecule 
AE' = #' ~— E° willinterest us. If we designate the excitation energy of 
the Ith level in the phenyl group by AE!(= ~AK}), then we arrive at the 
following two values for the excitation energy corresponding to AE! of 
the biphenyl molecule: 


AEt 
AF} 


ARE + B+ Mi, 
AEL + Bi — Ml}, 


(10.25) 


where 


Bi = fe (Tet (T) Vep* (II) 6) (1) dr 
— fee* (1) 62(1) Vb)* (IT)? (II) dr — exchange terms. (10.26) 


B' is characterized by a general displacement of the excitation energy of 
the bipheny] molecule with respect to the excitation energy of the phenyl 
group (benzene). For the first excited states of the molecule, this term 
is small, since the difference in the energy for the interaction of an excited 
and normal “benzene” molecule with another normal molecule is, in this 
case, small compared with M’, which characterizes the resonance inter- 
action between benzene molecules. Henceforth, we shall assume that 


B' = 0. 


Substituting (10.24) in (10.21), we obtain two values for the ratio of the 
coefficients ¢. and ¢; accordingly, the wave functions of the excited states 
[(10.25)] of the biphenyl molecule will have the form: 


Di = 2-4 +h), = Dh = 2-4 — Bb). (10.27) 


In order to determine the magnitude of the splitting of corresponding 
levels, the matrix elements [(10.23)] must be calculated. The function 
V entering into these integrals and forming the classical interaction energy 
between two neutral molecules may be represented in the form of a sum of 
terms which takes into account the dipole-dipole, dipole-quadrupole, 
and other types of interaction: 

6 


12 
2 
Vor ai > (QYaiYr; — Laitns — ZaiZoj) Sper at (10.28) 
i=1j=7 

where RF is the distance between the centers of the molecules; 2a:, Yai, 
. , 2, . . . are the coordinates for the electrons of the ‘‘benzene” 
molecules a and } relative to the systems of coordinates located on the 
centers of the molecules, with the y axis directed along a line which links 

the centers of the molecules; and e is the electronic charge. 
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If we substitute (10.28) in (10.23), we obtain for the dipole-dipole 
interaction 


12 
H 


MW=s — ai J @'*(r) p° sal (2 cos 6% cos 6% — cos 6 cos OF 


— cos 6 cos 67) — > if &!* (1) 9(I1) Vee* (II) Gi (I) dr, (10.29) 


6 
where r= > Lat, > yi = 7 cos OY, etc. 
i=l 


For forbidden dipole transitions, 
f@'*(r)b° dr = 0. 


It is then necessary to take into account in (10.28) the interaction of a 
higher order. For the first two forbidden transitions in the benzene 
molecule, transition matrix elements different from zero correspond to the 
octupole transitions only. Consequently, the terms of the octupole- 
octupole interaction must be taken into account in the matrix element. 

According to Goeppert-Mayer and Sklar,5’ the lowest singlet excited 
state of the benzene molecule belongs to the irreducible representation 
By, of the Dg, point group.'°* The transition from the normal state A1, 
to the level Bo, is forbidden. In the biphenyl! molecule, two wave func- 
tions will correspond to this level: 


Dr= 2+), Di=2*Ga—H), (10.80) 


where k} and k} are determined through ! and ®°, with the aid of (10.20). 
Accordingly, we obtain for the excitation energy 


AE} = AEL+ M1, AE} = AE} — Mm. (10.31) 


By utilizing the table of irreducible-representation characters for the 
symmetry group Ds, (Table 2), to which the biphenyl molecule belongs, 
we can, by comparing the characters, determine the irreducible repre- 
sentations of the functions given by (10.30) if we consider the symmetry 
properties of benzene wave functions. 

We will obtain: 

D} transforms as Biz 
Di transforms as Ao,. 


Knowing the symmetry of wave functions, one can determine the selec- 
tion rules and the polarization of the corresponding quantum transitions 
by using the methods of group theory. Thus we find that the first 
forbidden transition of the benzene is split into two in the bipheny] 
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molecule. The transition to one of the split levels becomes permitted 
and will have the polarization of an electric vector along the intermediate 
axis of the molecule (axis x in Figure 25). 

The second forbidden transition (to B,,) in the benzene also splits into 
two in the biphenyl molecule: 


AER =A\E2+M? (Bu), 


7) 
AE} = AE? — M* = (Ai), (10.32) 


with wave functions 


Di= 2k +) (Bu), 
Di = 28(K — 1) (An), 


respectively. The transition to the state characterized by the function 
D? is permitted and possesses the polarization of an electric vector along 
the long axis of the molecule (axis y in Figure 25). 

The transition to the third, doubly degenerate level £1, in the benzene 
molecule is permitted. In the benzene molecule, two wave functions 
correspond to it: 

WV; and W, = ws (Fi.). 


In lieu of these functions, it is more convenient to take the linear 
combinations 
pm 2h OS) (Bi), 
2 


B} = (5 ‘= ve) (Boy). (10.33) 


In the biphenyl molecule, two wave functions will correspond to the 
benzene function 3 ((10.33)]: 


Di = 2% (ki + hi) (Bau), 


Di = 2-M(K — Rf) (Au). ak, 
The functions given by (10.34) correspond to the energies 
das 3 3 
AB} = AF3 + M (Bu), (10.35) 


AE} = AR} — M3 (Aig). 


The transition to the level AZ? is permitted and possesses a polarization 
of an electric vector along the long axis of the molecule. The transition 
to the level of AE} is forbidden. 

In the biphenyl molecule, two functions will correspond to a benzene 
function &* [(10.33)]: 


Di 
D} 


27 8(ke + hf) (Baw), 


ok — hf) (Ag) at0i88) 
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these have energy values 


AEt = ARS + M4 (Bi); 
AFA = AE‘ —M* (Ay). 0087) 
The transition to the level AZ} is permitted and its polarization will 
correspond to an electric vector along the x axis of the biphenyl molecule. 
The transition to the level AF} is forbidden. 

Although direct transitions from the normal state to states with the 
energy AEi, AH}, AH}, and AE% are forbidden in the biphenyl molecule, 
they may occur in the absorption spectrum in combination with vibra- 
tional frequencies of a certain symmetry. These absorption bands must 
be weak, since they will disappear if the eigenfunction of the electronic- 
vibrational state is represented in the form of the product of the 
electronic and vibrational functions. 

b. Calculation of Matrix Elements. The matrix elements M', which 
determine the magnitude of the splitting of energy levels in the biphenyl 
molecule according to (10.29), are expressed by two terms. The first 
term of (10.29) characterizes the exchange of excitation (excitons) between 
the phenyl groups aand 6. It depends essentially on the character of the 
corresponding transition in the phenyl group. The second term deter- 
mines the interaction between phenyl groups, in connection with which 
the exchange of pairs of electrons takes place between aand b. In order 
of magnitude this term will be equal to the interaction with an exchange 
of electrons between the two closest carbon atoms belonging to different 
phenyl groups. This term depends little on the character of the transi- 
tion; therefore, we may write 


» [ emsanvesrcanaa dr = = S?*+ A —2SB =~ —0.18 ev, 
ap 
“u 
(10.38) 
where 


se i Us (1) Us(1) dr = 0.224, 

i | Uz (1) Us(1) = U-(2)U$(2) dri dre ~ 0.68 ev, 
12 

B= | U* (1) Ua(1) = de, = 3.01 ev. 
1p 


U,(1) and Ug(1) are the wave functions of the first 7 electron in the a and 8 
carbon atoms, respectively (see Figure 25), with a screening constant 
z = 3.18; rog = 1.48 X 10-° cm is the distance’ between the atoms a 
and £8. The value of the exchange integral A is taken from the tables 
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given by Sklar and Lyddane.'"? The integral B has been investigated by 
H. H. Hellman (Ref. 108, p. 510). 

If we now use the values 6! and &*, corresponding to the forbidden 
transitions in the benzene molecule, 


1 
Qh = Vi (Wi — ¥f), 


res 


G2? = y Ws 
/2 ( 1 + TY), 
where 
Wi = go0(1) 0(2)91(8) ¢1(4) ¢-1(5) g2(6), 
we. = 67% Y e2rki6 UI, 
B 
and B° = go(1) g0(2)¢1(3) 1 (4) y-1(5) e-1(6), 
we obtain 
@ iu flare _ {[ 9, if 2 = 1, 
R Je Ne tool ev, ifl=2 
Thus 
M = 0.18 ev, M? = 0.17 ev. (10.39) 


In order to calculate M$ according to formula (10.29), we must first of all 
determine 


Jfb2*(r) 6° dr. 


Substituting the value $3 from (10.33) into this integral and using the 
wave functions of the benzene molecule, 


Ws = VT = go(l)go(2)e1(3)¢1(4) ¢-1(5) e_2(6), (10.40) 
we obtain 
ak og, — 4P_ 
- BCE) de = ST Oy, (10.41) 


where p = 1.39 X 10-8 cm and ry, is the unit vector in the direction of the 
y axis. 
After substituting (10.41) and (10.38) into (10.29) we obtain 


M? => —0.46 ev. (10.42) 
M* is expressed through the transition integral {6**(r)®° dr. Sub- 
stituting 4 from (10.33) and using (10.40), we obtain 


i b*(r)° dr = jee, (10.43) 
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where M1, is the unit vector in the direction of the x axis of the bipheny] 
molecule. Using (10.43) and (10.38), we obtain from (10.29) 


ALl* = 0.42 ev. (10.44) 


c. Comparison with Experiment. There are few experimental studies on 
the absorption spectra of biphenyl pairs. Carr and Stiicklen!®*.10 dis- 
covered a strong absorption in the region of 6.4 ev, as well as a wide 
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Fie. 26. Schematic representation of the ultraviolet absorption bands of biphenyl (I) 
and benzene (III) vapors correlated with calculated energy terms. 


absorption maximum from 7.3 up to 7.6 ev. Because of the conditions of 
the experiment, the higher frequency portion of the spectrum was not 
investigated. Besides, A. London, shows in his work*®® that Beck at the 
University of California discovered in biphenyl pairs a wide absorption 
band in the region 2200 to 2800 A (5.6 to 4.4 ev). Dickson"! also gave 
indications of an absorption band in the vicinity of 4.4 ev. 

Figure 26 (curve I) represents the experimental absorption spectrum 
according to the work of Carr and Stiicklen.‘°* The crosshatched region 
corresponds to Beck’s data. We do not have data on the distribution of 
absorption intensity in this portion of the spectrum. 

Curve III represents the benzene absorption spectrum corresponding 
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to transitions to the levels Bi, and Ey, according to the works of Carr 
and Stiicklen.!9%1!0 Consequently, for benzene levels, we may adopt the 
values 


6.1 ev (Bix), 
AEt=6.9er (Ew). 


AE? 
ABS 


The frequency 38,089 cm~! corresponds to the first transition in benzene 
(according to A. Y. Hichis’ data); consequently, 


AF} = 46 ev (Boy). 


Using these values for the differences (above the normal level) of energy 
levels in the benzene molecule, we obtain, with the aid of (10.31), (10.32), 
(10.35), (10.37), (10.39), (10.42), and (10.44), the energy levels (relative to 
the normal one) of the bipheny! molecule. The position of these levels 
is shown in Figure 26 on line II. 

The terms corresponding to the forbidden transitions from the normal 
level are shown by broken lines, while the permitted ones are indicated by 
solid lines. 

Line IV gives the splitting of the terms for the benzene molecule accord- 
ing to A. London’s theory.*® We must note that, for the benzene levels, 
London adopted the levels calculated in the work of Goeppert-Mayer and 
Sklar®? (5, 5.6, and 8 ev). This seems inconsistent to us, since the 
improved calculation of the benzene molecule, as London himself notes, 
leads to energy levels of 7.3, 8.8, and 1].1 ev. By virtue of the better 
agreement with the experiment, London justifies the use of the values 5, 
5.6, and 8 ev from an older work. However, in order to be more con- 
sistent, one should then simply take the experimental values, as we do 
here. Besides, London’s selection does not improve the agreement with 
the experiment, since all terms obtained from the splitting of the per- 
mitted 8 ev then go to the far ultraviolet. The strong absorption in the 
6.4-ev region is, according to London, explained by the splitting of the 
forbidden benzene level (5.6 ev). 

According to our theory, this absorption band mainly originates from 
the splitting of the permitted term of the benzene molecule. As a result 
of this splitting, a displacement of the permitted transition to the red side 
by 0.46 ev takes place. The second split component of this same transi- 
tion (displacement of 0.42 ev) apparently corresponds to the absorption 
band of biphenyl in the 7.2-ev region. 


CONCLUSION 


The purpose of the present work was to investigate the spectra of light 
absorption by molecular crystals. It was assumed here that the corre- 
sponding spectra of light absorption by free molecules (spectra of rarefied 
gases) were completely known. The solution of the problem had to link 
absorption spectral changes upon transition from a gaseous to a crystalline 
state with the crystalline structure of a solid and with the forces of inter- 
action between molecules. 

Initially, the problem was solved under the assumption that all mole- 
cules of a crystal are rigidly fixed in the points of the crystal lattice. 
Classical mechanics was used in studying the influence of forces of inter- 
action between the molecules in a molecular lattice upon intramolecular 
vibrations (Chapter II). Then we developed the quantum theory of these 
changes, as well as the theory for the change of energy states which 
include electronic excitation. The classical and quantum theories yield 
identical qualitative results, which essentially amount to the following: 

1. In a crystal, each excited state of a molecule is transformed into a 
series of quasi-continuous zones of excited states. The number of zones 
is equal to the number of molecules in the crystalline unit cell. The 
number of sublevels in each zone is equal to the number of unit cells in the 
crystal. The sublevels inside each zone differ in the values of the wave 
vectors for the excitation waves (excitons) that spread in the crystal 
(Sections 7, 9, and 12). 

2. Upon interaction with light only excitons with a wave vector equal 
to the wave vector of a light wave are excited. With the aid of group 
theory we established those quasi-continuous zones of excited states to 
which transitions take place under the action of light; the latter is 
propagated in a given direction and has a given polarization (Sections 
8, 13, and 14). 

3. As a result of the emergence of several zones of excited states corre- 
sponding to one excited state in the isolated molecule, there occurs a dis- 
placement and splitting of the crystalline absorption band in comparison 
with that of a gas. The number of split components and their polariza- 
tion are determined by the number of molecules per unit cell and by their 
mutual orientation. 

Until recently, the opinion generally held was that the splitting of 
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absorption bands in a crystal could only be caused by (a) polymerization 
or association of molecules or (b) removal of degeneracy if the molecule 
possessed degenerate levels. We have shown, however, that the splitting 
of absorption bands may take place in a crystal even when these two 
reasons are lacking. The magnitude of splitting depends mainly on the 
resonance interaction of the molecules in a crystal. It has the greatest 
value for dipole transitions (proportional to the oscillator strengths) and 
decreases rapidly with the increase of the multipole index for the transi- 
tions (Sections 9 and 12). 

4, When the theory is applied to combinational scattering (Raman 
scattering) of light, it leads to the conclusion that, with regard to the 
influence of a crystalline state on intramolecular vibrations, the latter may 
be divided into two groups: (a) vibrations noticeably influenced by 
crystalline structure (these vibrations must be simultaneously active in 
the infrared spectrum) and (6) vibrations not noticeably influenced by 
crystalline structure. Both these groups of vibrations may be simul- 
taneously present in one molecule (Section 10). 

5. We developed the general absorption spectral theory, which takes 
into consideration the possibility of displacements of molecules from their 
equilibrium positions. We pointed out the inadequacy of generally 
accepted conditions with respect to the accuracy of the adiabatic approxi- 
mation when the latter is applied to a system with many degrees of 
freedom. In this case, an essential role is played by the temperature of 
the body and by the correlation between time of molecular displacement 
from equilibrium positions and time of excitation transfer from one 
molecule to another. Excited states of a crystal (with one excited mole- 
cule) prove to be nonstationary. Only in the case in which excitation 
spreads over the crystal so rapidly that there is no time for a local crystal- 
line deformation to take place may we speak about quasi-stationary 
excited states (with respect to transformation of molecular excitation 
energy into lattice vibrational energy). If, however, excitation spreads 
slowly, a local crystal deformation occurs (the exciton is “localized”’). 
In this case, the probability of molecular excitation energy transformation 
into lattice vibrational energy is great (Sections 15 and 16). 

6. We analyzed the possibility for excitation in a crystal of two limiting 
types of excitons: (a) free excitons, which spread rapidly through the 
crystal and do not cause crystalline deformation, and (0) “localized’’ 
excitons, which move slowly because their movement is linked to the 
movement of the local crystal deformation (Section 17). 

7. Equations were obtained which determined the movement of 
molecules in a crystal upon free exciton formation. Their solutions in the 
zeroth approximation result in a crystal excitation energy value, which 
coincides with the case of fixed molecules in lattice points. This estab- 
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lishes the limits for the applicability of the solutions which were obtained 
in Chapter III when studying fixed molecules (Sections 18 and 19). 

8. The group velocities of excitation waves for different excited-state 
zones, but for identical wave numbers, were shown to be not identical. 
For certain excited-state zones, these velocities (for k = Q) may prove so 
small that excitation by light of this type of free exciton becomes impossi- 
ble (localized excitons will be excited) (Section 22). 

9. The proposed theory explains the pleochroism of molecular 
crystals and constitutes the first attempt to form a molecular theory of 
pleochroism in contrast to the usual phenomenological description of 
this phenomenon. Pleochroism of a molecular crystal is caused by a 
property (arising upon interaction in a specific manner of oriented mole- 
cules) of the crystal as a whole, and not by a property of separate mole- 
cules (Sections 22 to 24). 

10. The solution of equations, determining the movement of molecules 
in a crystal upon formation of free excitons, in the first approximation 
permits us to describe the processes of free exciton excitation simul- 
taneously with lattice vibrational excitation. Formulas determining 
absorption band structure and its dependence on temperature are 
deduced. It is shown that the absorption-band width increases with the 
temperature. As a result of interaction with lattice vibrations, absorp- 
tion may attain several maximum values in the absorption band, which 
corresponds to one excited state of a molecule. However, the distances 
from these maxima to the primary band (excitation of a crystal without 
lattice vibrations) may not coincide with the lattice vibrational fre- 
quencies. In certain cases, upon lattice vibrational excitation, maxima 
may appear on the long-wave side (for T’ = 0) of the primary absorption 
band. On the contrary, in luminescence spectra of crystals at low tem- 
perature, one may directly observe lattice vibrational] frequencies, which 
will occur in the form of satellites near the basic lines (Sections 25 to 27). 

11. We have pointed out possible changes in absorption spectral 
structure upon temperature rise and upon translational symmetry viola- 
tion in a crystal. Temperature increase leads to a general intensity 
increase for absorption bands, which are caused by excitation of molecules 
in a crystal upon simultaneous excitation of lattice vibrations. It also 
leads to the appearance of additional absorption bands, conditioned by 
the excitation of a crystal from already-excited vibrational states of a 
lattice. The violation of ideal periodicity of a lattice by the presence of 
admixtures and isotopes and, finally, by thermal fluctuations of density 
will lead to the following: A portion of the ‘‘quasi impulse’’ will be trans- 
mitted to the corresponding heterogeneity. The absorption-band width 
will increase asa result. This broadening effect of absorption bands must 
play an essential role in absorption of light by a liquid (Section 27). 
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12. Upon formation of “localized” excitons, the corresponding excited 
states of a crystal are not stationary. Owing toa local deformation of the 
crystal, the connection between intramolecular excitation and lattice 
vibrations increases. The excitation of molecules is transmitted to lattice 
vibrations. On the one hand, this explains the wide absorption bands 
upon excitation of localized excitons; on the other hand, it results in the 
transformation of absorbed energy into heat. In this case, the widening 
of absorption bands also takes place, because in addition to a ‘‘localized” 
exciton, there may be simultaneous excitation by light of a few phonons 
as well. Upon formation of localized states, the distances between 
various quasi-continuous zones of excited states, corresponding to one 
excited state of a molecule, are small and usually overlap owing to the 
large spread of absorption bands. Therefore, polarization effects would 
play a small role in these cases (Sections 29 and 30). 

13. It was shown that, in a crystal, free excitons may spontaneously 
turn into “‘localized” ones if the energy of the free exciton is greater than 
the energy of the ‘‘localized” one; also, upon decrease in the energy of the 
free exciton, the velocity of the latter decreases. 

Conversely, in certain cases the energy of a “localized” exciton is higher 
than the energy of a free exciton, and the energy decrease of a free exciton 
is linked with the increase of its velocity. In this case, “localized” 
excitons turn into freeones. This last case apparently plays an important 
role in a series of molecular crystals. The velocity increase (and that of 
the wave number) for an exciton upon decrease of its energy may (a) 
condition the presence of luminescence of a solid (in a pure state) ; (b) lead 
to a lifetime increase for the excited state of a crystal in comparison with 
the lifetime of a free molecule (for T = 0); (c) be the cause of an incom- 
plete coincidence of the absorption band and of the luminescence band, 
which corresponds to a purely electronic intramolecular excitation 
(Sections 29 and 32 to 34). 

14. The process of transformation of intramolecular excitation energy 
into heat (into lattice vibrational energy) isextremely probable for excited 
states of a crystal that correspond to the excitation of “localized” excitons. 
If, at the moment of light absorption, free excitons are formed, and if the 
excitation energy of the latter is greater than the excitation energy of 
“localized” excitons, the process of transformation of free exciton energy 
into heat takes place through a state which corresponds to a “‘localized”’ 
exciton. 

The lattice vibrational frequency differs less from intramolecular vibra- 
tional frequencies than from electronic excitation frequencies. Owing 
to this, one must assume that the probability for the transfer of intra- 
molecular vibrational energy to heat will be greater than the probability 
for a transformation of electronic excitation energy. Apparently, this 
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explains the experimentally observed rapid transformation of the vibra- 
tional portion of molecular electronic-vibrational excitation energy to 
lattice vibrational energy (Sections 32 and 33). 

15. An analysis was made of conditions conducive to the emission of 
luminescence in molecular crystals. It was shown that luminescence is 
basically determined by transitions from excited states, corresponding to 
free excitons, to the ground state of a crystal (Section 34). 

16. Using certain experimental data by way of example, we have 
qualitatively confirmed the theoretical deductions of this work for 
combinational scattering (Raman) spectra, infrared spectra, and absorp- 
tion and luminescence spectra in the visible and ultraviolet regions. All 
of the experiments described confirm beyond any doubt the presence in a 
crystalline state of sharply polarized absorption bands which are new in 
comparison with the gaseous state (Sections 36 to 38). 

17. We have shown that exciton representations may be utilized to 
calculate the energy states of several (isolated) polyatomic molecules. 
We have found the dependence of the frequency for light absorption by 
para-polyphenyls on the number of phenyl groups in a molecule. We 
have calculated the spectrum of light absorption by the biphenyl molecule 
(Section 39). 
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